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HAUPTAUFSÄTZE 


Zur Stabilität der Karmänschen Wirbelstraße. 
Von A. W. Mauwe in München. 


er ursprüngliche Plan der Arbeit war, die Kärmänschen Stabilitätsuntersuchungen über 

Wirbelstraßen'!) auf neuem Wege zu wiederholen. Dabei sollte, um das etwas schwer- 
fällige Rechnen mit unendlichen Summen zu vermeiden, voller Gebrauch von den funktionen- 
theoretischen Hilfsmitteln gemacht werden, die für die Behandlung des Problems zur Ver- 
fügung stehen, aber bisher nicht ausgenützt worden sind. Im Laufe der Durchführung dieser 
Aufgabe zeigte es sich, daß der von Kärmän betrachtete Fall einer stabilen Wirbelstraße 
nur ein Sonderfall aus einer ganzen Schar ähnlicher Wirbelanordnungen ist, die ebenfalls 
stabil sind. Die vorliegende Arbeit wird daher diese allgemeineren Wirbelstraßen mit in ihre 
Betrachtungen einbeziehen. Um deutlich zu machen, in welcher Richtung die Verallgemeinerung 
liegt, nehmen wir die wichtigsten Angaben über Inhalt und Ergebnisse der Arbeit vorweg. 


Wir betrachten eine ebene Flüssigkeitsströmung, die, abgesehen von geraden Wirbel- 
fäden senkrecht zur Strömungsebene, wirbelfrei ist. Diese Wirbelfäden, kurz Wirbel genannt, 
seien in zwei parallelen unendlich langen Reihen angeordnet. Innerhalb beider Reihen folgen 
die Wirbel in gleichem Abstand aufeinander. Diesen Abstand nennen wir die Teilung der 
Wirbelreihen. Die Wirbel haben alle gleiche Stärke und in beiden Reihen verschiedenen 
Drehsinn. 


Beide Wirbelreihen können in verschiedener Weise gegeneinander angeordnet sein. 
Kärmän betrachtet die beiden in den Bildern 1 und 2 dargestellten Anordnungen. Die erste 
erweist sich als stets labil, die zweite ist stabil für ein bestimmtes Zahlenverhältnis von 


Straßenbreite h zu Teilung ! der Wirbelreihen. Dieses Verhältnis ist 20,38 und ist fest- 


gelegt durch die Beziehung 


9) Th. v. Kärmän und H. Rubach: Phys. Z. Bd. 13 (1912), S. 49. — L. Prandt]: Abriß der Strömungslehre, 
Braunschweig 1931, Abschn. III, $ 12. — Prandtl-Tietjens: Hydro- und Aeromechanik, Berlin 1929 u. 1931, 11. Bd,, 
S. 155. — W. Kaufmann: Hydromechanik, Berlin 1931 u. 1934, I. Bd., S. 191. — Lamb: Hydrodynamies, 4. Aufl. 
Cambridge 1916, S. 218. — Frank-v. Mises: Differentialgleichungen der Physik, 2. Aufl., II. Bd., Braunschweig 1935, 
S. 465. 


9 


— 


130 Maue, Zur Stabilität der Kärmänschen Wirbelstraße 1946 


Hier betrachten wir den allgemeinen Fall von Bild 3, in dem d eine beliebige Strecke 
mit einer Länge zwischen 0 und / ist. Die Stabilitätsuntersuchung liefert für jeden d-Wert 


ein Verhältnis 3 ‚ für das die Wirbelanordnung stabil ist. Es ist gegeben durch die Beziehung 


han . dn 
Siny =siny 
Bild 1. Wirbelstraße, erste von Bild 2. Wirbelstraße, zweite von Bild 3. Wirbelstraße, 
Kärmän betrachtete Anordnung. Kärmän betrachtete Anordnung. allgemeine Anordnung. 


Unsere allgemeine Betrachtung umschließt die beiden vorgenannten Fälle als Sonder- 
fälle. Für d=5 geht Bild 3 in Bild 2 und Gl. (2) in Gl. (1) über. Für d=0 (Bild 1) folgt 


aus Gl. (2) h=0; d.h. die Wirbelstraße artet aus, indem die Wirbel beider Reihen zusammen- 
fallen. Die Anordnung von Bild 1 mit endlicher Straßenbreite h ist also jedenfalls nicht 

stabil, in Übereinstimmung mit der Kärmänschen Feststellung. 
h : Der in Gl. (2) angegebene Zusammenhang zwischen 
d und h ist in Bild 4 graphisch wiedergegeben. Der 


Anordnung d=4 entspricht danach die größte 


Straßenbreite h bei gegebener Teilung 1. 


. Es sei noch bemerkt, daß die beiden Fälle 


4 d=0 und dadurch ausgezeichnet sind, daß 
Bild 4. hische D ie Wj 
sich die Wirbelstraße, d. h. die ganze Wirbelanord- 
nung, in ihrer Längsrichtung bewegt, wenn wir 
voraussetzen, daß die Flüssigkeit im Unendlichen ruht. In allen anderen Fällen erfolgt die 
Fortbewegung der Wirbelstraße in schräger Richtung. 


$ 1. Die ungestörte Wirbelstraße. 


1. Die ebene Potentialströmung. Wir benutzen die Theorie der ebenen Potentialbewegung. 
Die Strömung erfolge in der komplexen z-Ebene mit z=x-+iy. Die komplexe Geschwindig- 
keit sei w=u-+iv, die zugehörige komplexe Strömungsfunktion wobei 
das Geschwindigkeitspotential und y die Stromfunktion ist. Zwischen f(z2) und w besteht 
die Beziehung 
„_4f 


Der Stern bedeutet Übergang zum konjugiert komplexen. 
2. Strömungen mit Wirbeln. Die Strömungsfunktion 


beschreibt eine Strömung, die in z=z, einen Wirbel mit positivem Drehsinn und der Wirbel- 
stärke /' besitzt. 

Ist g(z) eine ganze Funktion mit nur einfachen Nullstellen, so besitzt die Strömung 


T 
in allen diesen Nullstellen Wirbel der Stärke I’ mit positivem Drehsinn. Führen wir zwei 
ganze Funktionen mit nur einfachen Nullstellen g, (z) und g, (2) ein, so können wir die Strömung 
angeben, die in den Nullstellen von g, positive und in denen von g, negative Wirbel, alle 
von der gleichen Stärke /‘, besitzt. Sie lautet: 


(6). 


2 ni 
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3. Fortbewegung der Wirbel. Während die Strömung eines einzelnen Wirbels stationär 
ist, gilt das für Strömungen mit mehreren Wirbeln im allgemeinen nieht mehr. Vielmehr 
bewegt sich jeder Wirbel unter dem Einfluß der anderen Wirbel; er nimmt an der Bewegung 
des Geschwindigkeitsfeldes der übrigen Wirbel teil?). 

Wir wollen die Geschwindigkeit berechnen, mit der sich ein Wirbel der Strömung (6) 
fortbewegt. Ein positiver Wirbel befinde sich etwa in z2=z2,. Dann gilt 


Wir entwickeln f(z) inz=z,. Die Entwicklung lautet mit Rücksicht auf Gl. m 


2,)+1g 9,(2,) 2g',(z,) 9: (2,) (2 - 2,)+- :lg(@—z, )+f,(2) (8). 
Im letzten Ausdruck von Gl. (8) haben wir f(z) in zwei Teile gespalten. Der erste beschreibt 
das vom Wirbel z=z, herrührende Geschwindigkeitsfeld, der zweite f, (z) das Feld der übrigen 
Wirbel. Das letztere ist für die Fortbewegung des Wirbels in 2, maßgebend. Wir erhalten 
die Geschwindigkeit w, dieser Bewegung, wenn wir Gl. (3) auf f,(z) anwenden und 2=z, 


setzen: 

Entsprechend ergibt sich die Geschwindigkeit w, eines in 2, gelegenen negativen Wirbels: 

Gl. (10) geht aus GI. (9) durch Vorzeichenumkehr und Vertauschung der Indizes 1 und 2 hervor. 
” 4. Anwendung auf die Wirbelstraße. Wir 
gehen zur Betrachtung der Wirbelstraße (Bild 5) 
über, die in den Punkten (n reell ganz- 
zahlig) positive, in negative Wirbel 


gleicher Stärke I’ besitzt. Die zugehörige Strö- 
mungsfunktion erhalten wir, wenn wir in Gl. (6) 


Bild 5. Lage der Wirbel in der komplexen Ebene. 


w, 


(10). 


(11) 


setzen. Diese beiden Funktionen besitzen in den 
geforderten Wirbelpunkten Nullstellen. Es ergibt sich 


sin 


(12). 


f@d)= 


Die Strömung (12) ist außer durch ihre Wirbelanordnung dadurch gekennzeichnet, daß die 
Flüssigkeit im Unendlichen ruht. Davon überzeugen wir uns, indem wir 2 gegen +i co gehen 


Um die zeitliche Änderung des Wirbelbildes zu bestimmen, berechnen wir die Ge- 
schwindigkeit des Wirbels 2=2,=z, nach Gl. (9) und die des Wirbels z2=2,=—z, nach 
Gl. (10) unter Verwendung der Funktionen (11). Für w, und w, ergibt sich derselbe Wert, 
den wir w, nennen, und zwar ist 
2 


Mit den Wirbeln in 2, und —z, kommt wegen der Symmetrie der Anordnung auch allen 
andern Wirbeln die gleiche Geschwindigkeit zu. Das ganze Strömungsbild bewegt sich also 
mit der Geschwindigkeit w,. 


$ 2. Ableitung der Stabilitätsbedingung. 


Um die Stabilität unserer Strömung zu untersuchen, haben wir uns die Wirbel um 
beliebige kleine Strecken aus ihrer Lage verschoben zu denken und die so entstehende 
gestörte Strömung zu betrachten. Für diese die Strömungsfunktion zu finden, ist unsere 


2) Vgl. z. B. Prandtl-Tietjens, I. Bd., S. 194. 


lassen. f(z) strebt dann nämlich dem konstanten Grenzwert — 5 (mod. IT) zu. 
| 
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nächste Aufgabe. Während die funktionentheoretische Beschreibung mittels einer komplexen 
Strömungsfunktion für die ungestörte Wirbelstraße (Gl. (12)) allgemein üblich (a. a. O.) ist, 
ist diese Art der Darstellung für die gestörte Strömung nicht gebräuchlich. 

Zur Vereinfachung der folgenden Rechnung setzen wir in $$2 und3 !’=2a und !=1. 


1. Die gestörte Wirbelstraße. Die Störung berücksichtigen wir, indem wir in Gl. (6) für 

9,(2), 9,(z) statt der Ausdrücke (11) jetzt setzen 

g,(2)= sin (2 — —n-[,(2) 

Die analytischen Funktionen £,(z) und £,(z) werden dabei im Gebiete der Wirbelstraße, d. h. 
in einem Streifen längs der reellen Achse, als dem Betrage nach klein gegen 1 vorausgesetzt. 
Auf diese Weise erreichen wir erstens, daß sich die Nullstellen von g,(z) und g,(z) 
und damit die Wirbel nur wenig gegenüber ihrer ungestörten Lage verschieben, und 
zweitens, daß im Gebiete der Wirbelstraße keine neuen Nullstellen hinzukommen. 
Um auch im übrigen Teil der komplexen Ebene beiderseits der Wirbelstraße das Auftreten 
neuer Nullstellen zu verhindern, genügt es, die Beträge von /, und £, als klein gegenüber 
den (mit dem Abstand von der Wirbelstraße wachsenden) Beträgen der Sinusfunktionen in 
Gl. (14) anzunehmen. Da in den folgenden Rechnungen nur die Werte von Z, und £, im 
Bereich der Wirbelstraße eine Rolle spielen, dürfen wir quadratische Glieder in £,,£, ver- 
nachlässigen. Insbesondere entspricht es einer Vernachlässigung quadratischer Glieder, wenn 
wir zwischen den Funktionswerten von £,,{, in den unverschobenen und den verschobenen 
Wirbelpunkten keinen Unterschied machen. 

Wir berechnen die Verschiebungen der Wirbel gegenüber ihrer ungestörten Lage, wie 
sie durch Gl. (14) gegeben sind. Die unverschobenen Wirbel der oberen Reihe liegen in 
den Punkten 2,+n. Nennen wir die Verschiebung des n-ten Wirbels s,“, so ist seine 
gestörte Lage z2,+n+ s,. Wir entnehmen s„“” aus der Gleichung 


Die Ausrechnung liefert unter Vernachlässigung quadratischer Glieder in £, und s„“ 
also 


Entsprechend ergibt sich für die Verschiebung s„‘” des n-ten Wirbels der unteren Reihe, 
dessen unverschobene Lage — 2, + n ist, 


2. Die zeitliche Änderung der Wirbelverteilung. Wir untersuchen die Fortbewegung des 
Wirbels, dessen ungestörte Lage 2, war. Hierzu verwenden wir Gl. (9), in der wir 2,=2, 
+3,” —=2,+{[,(z,) setzen. Zunächst berechnen wir die benötigten Funktionswerte von g,, 9; 
und ihren Ableitungen und vernachlässigen, wie vorgesehen, quadratische Glieder in £,, £.. 
Wir erhalten 


(14). 


(18). 


&,, £, und ihre Ableitungen sind dabei an der Stelle z, zu nehmen, wofür wir, wie oben be- 
gründet, in diesen kleinen Funktionen auch einfach 2, setzen dürfen. Gl]. (9) liefert schließlich 
mit und !=1 unter neuerlicher Vernachlässigung quadratischer Glieder in die 
Geschwindigkeit w, des betrachteten Wirbels. Es ergibt sich 


Das erste Glied stimmt wegen ’’=27 und !=1 mit w,* von Gl. (13) überein, entspricht also 
der Bewegung der ungestörten Wirbelstraße. Die restlichen Glieder zeigen eine zusätzliche 
Bewegung des betrachteten Wirbels an, entsprechen also einer zeitlichen Änderung der Ver- 
schiebung s,”—=[,(z,). Es gilt demnach, wenn wir die zeitliche Differentiation durch einen 
Punkt andeuten, die Beziehung 

n? 
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Eine entsprechende Beziehung muß auch für die übrigen Wirbel der oberen Reihe gelten. 
Da diese in keiner Weise vor dem betrachteten ausgezeichnet sind, gilt für sie dieselbe Be- 
ziehung (20) mit dem einzigen Unterschied, daß jetzt als Argument von £, und £, statt 2, bei 
Betrachtung des n-ten Wirbels z,+n zu verwenden ist. In den vorkommenden trigonometri- 
schen Funktionen sin und ctg ist diese Ersetzung von z, durch 2,+n belanglos. Wir können 
demnach zu Gl. (20) als Ergänzung hinzufügen: 


u=f, (2.+n) 
+9) 


In derselben Weise kann die untere Reihe untersucht werden. Wir haben hierzu statt 
Gl. (9) die Gl. (10) zu verwenden oder, was dasselbe ist, auf der rechten Seite von Gl. (9) 
das Vorzeichen umzukehren und die Indizes 1 und 2 sowohl der Funktionen g,, 9, als auch 
ihrer Argumente z,, 2, zu vertauschen. Diese Indexvertauschung ist gleichbedeutend mit 
einer Vertauschung von Z, mit {, und von 2, mit —2,. Statt der GIn. (20) und (21) ergibt 
sich so für die untere Wirbelreihe 


mit 
- +0 


Damit haben wir die Differentialgleichungen aufgestellt, die die zeitliche Entwicklung 

einer einmal hervorgerufenen Störung beherrschen; denn die Funktionswerte von £,, Z, in den 

Wirbelpunkten sind ja wegen der Gln. (16) und (17) bis auf den Faktor (— 1)” mit den Ver- 
schiebungen s,„‘”, identisch. 


3. Spektrale Zerlegung der Störung. Es besteht nunmehr die Aufgabe, die Funktionen 
,(z) und Z,(z) so zu bestimmen, daß sie beliebig vorgegebenen Verschiebungen s,„“, s,” der 
Wirbel entsprechen. Wir denken uns zu diesem Zwecke die Störung aus Störungswellen auf- 
gebaut und entwickeln in Fourierscher Weise 


(— 1)” . — (2, n) = \ einy dp 
(24). 


Dabei haben wir, um uns in den folgenden Rechnungen von dem Faktor (—- 1)" zu befreien, 
nicht selbst, sondern (— 1)"- entwickelt. Da es sich bei den Größen 
um die Verschiebungen diskreter Punkte handelt, ist der Integrationsbereich für endlich. 

Es lassen sich nun leicht Funktionen £, (2), £,(2) angeben, die in den Wirbelpunkten 
2,+n, —2,+n die durch Gl. (24) vorgeschriebenen Werte annehmen. Wir können nämlich 
einfach schreiben 


Das ist zwar nicht die einzige Möglichkeit, Z, (2) und £, (z) in Einklang mit Gl. (24) zu wählen. 
Doch ist diese Wahl dadurch ausgezeichnet, daß das Spektrum von £, (2) und £, (z) nur Wellen- 
zahlen |p|<x enthält. Wellenzahlen |p|> dürfen wir aber deswegen nicht zulassen, weil 
sie zur Folge haben würden, daß £,, £, für große positiv oder negativ imaginäre 2 größer als 
sin (2 # 2,)r sein würden, in Widerspruch mit unseren Forderungen von Ziffer 1. 


4. Schwingungsgleichungen für die Amplituden der Störungswellen. Wir gehen mit dem 
Ansatz (25) in Gl. (20) ein und berücksichtigen die Gl. (21). Es ergibt sich 


2 sin’2 72, 


Gleichsetzung der Koeffizienten von ei”% auf beiden Seiten von Gl. (26) liefert mit den Ab- 
kürzungen 
n? in’cos2nz, ap 
sin’ 2x2, sin?2n 2, 


3 (27). 
2 sin’272, 


Z. angew. Math. Mech. 
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Die Gleichung 

Ebenso ergibt sich, wenn wir von Gl. (22) ausgehen: 

u . . . 


Die Gl. (28) und (29) sind lineare Differentialgleichungen, die die zeitliche Änderung 
der Wellenamplituden y,, 7, bestimmen. 


5. Stabilität gegenüber Störungswellen mit 9 —=0. Besonders einfach liegen die Verhältnisse 
im Falle g=0, den wir zunächst betrachten wollen. 

Schreiben wir statt y, (0), x, (0) einfach y,,%,, so lauten die Gl. (28) und (29) mit Rück- 
sicht auf C=0 


Durch Addition und Subtraktion beider Gleichungen erhalten wir, wenn wir die Abkürzungen 


einführen und in der ersten der beiden entstehenden Gleichungen sogleich zum konjugierten 
übergehen, 


Setzen wir die Zeitabhängigkeit von o und ı* in der Form e’®t an, so haben wir 
Um © auszurechnen, bilden wir nach Gl. (27) mit Rücksicht auf = 0 die Ausdrücke 
in?’(1+ cos’2nz,) in? 


Einsetzen in Gl. (33) und Ziehen der Wurzel ergibt 


o=+ti 


2|sin 2 n2, 
Die beiden Werte für » unterscheiden sich nur durch das Vorzeichen. Sind sie reell, so 
besagt der Zeitfaktor ei”t, daß die Störungsamplituden o, 7* und damit auch y,, 2, oszillieren, 
d. h. die Wirbelpunkte kleine Schwingungen um ihre ungestörte Lage ausführen. Wir haben 
dann Stabilität. Sind die &-Werte dagegen wesentlich komplex, so entspricht der eine 
einem exponentiellen Abklingen der Amplituden o und ı* mit der Zeit, der andere einem 
exponentiellen Ansteigen. Es liegt Labilität vor. 


Demnach lautet die 
reell. 
Um diese Bedingung auf den Ausdruck (35) anwenden zu können, spalten wir das 
Produkt (1+cos2r2,*) (1 —cos27z,) in Realteil 1 —|cos27z,|’ und Imaginärteil cos 27 2,* 
—eos2rz,. Als Stabilitätsbedingung erhalten wir dann 


Um die Bedeutung dieser Forderung zu erkennen, setzen wir ei Bild as 


Dann erhält Gl. (36) wegen 


die Form 
........ (88). 
Drücken wir in dieser Gleichung mit Hilfe der Formeln 
die Kosinus durch die Sinus aus, so folgt 
und schließlich, da beide Sinus positiv sind, 
sinoda=Sinhn .... 


als endgültige Stabilitätsbedingung. Das ist mit Rücksicht auf Eh die gleiche Be- 
dingung, die wir am Anfang der Arbeit in Gl. (2) angegeben und besprochen haben. 
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6. Stabilität gegenüber beliebigen Störungswellen. Wir haben die Stabilitätsuntersuchung 
bisher nur für Störungswellen mit 9=0 geführt’). Es bleibt zu zeigen, daß Wirbelstraßen 
mit den durch Gl. (41) geforderten geometrischen Abmessungen gegenüber Störungswellen 
von beliebigem 9*) stabil sind. Der Beweis ist ziemlich langwierig und für das allgemeine 
Verständnis des Gegenstandes von zweitrangiger Bedeutung. Er soll aus diesen Gründen in 
einem gesonderten Abschnitt ($ 3) geführt werden. 

Der entsprechende Beweis für den Kärmänschen Sonderfall d -+ von Bild 2 ist bei 
Lamb (a.a.O., S. 222/223) angegeben, während er in der Kärmänschen Originalarbeit (a.a.O.) 
unterdrückt ist. 


$ 3. Störungen beliebiger Wellenlänge. 


1. Schwingungsgleichungen für die Amplituden beliebiger Störungswellen. Wir schreiben 
in den GIn. (28) und (29) zur Abkürzung 


u)=a, | 
und fügen durch Ersetzung von p durch —g und gleichzeitigen Übergang zum konjugiert 


komplexen zwei weitere Gleichungen hinzu. Das entstehende System von vier Differential- 
gleichungen lautet nach geeigneter Anordnung: 


a= 
“= A:a+(B +0 
—(B — C )e-?izp.a 


(42) 


(43). 


Der Zeitansatz ei“?! führt zu einem System gewöhnlicher linearer homogener Gleichungen. 
© wird durch die Säkulargleichung 


A* (B’ — C*) e2izo* y 
0 (B*+C*)e—2i20° 9 | 
| —0 (4) 
| A (B+ e?izor —im 0 
festgelegt. Die Ausrechnung der Determinante liefert mit Rücksicht auf 
(vgl. Gl. (37)) und mit den Abkürzungen 
Sin2hp—|4A|? 
— 2’ — (46) 
die Bestimmungsgleichung 
für Hieraus folgt 


Dabei ist nach den Gl. (46) p reell und g positiv reell. 
Stabilität liegt vor, wenn & reell oder, was dasselbe ist, 
©? positiv reell 


ist. Damit »’ reell ist, muß p? —:q? > 0 gelten; damit &® für beide Vorzeichen vor der Wurzel 
positiv ist, muß außerdem p>0 sein. Wir werden beweisen, daß 


woraus die beiden vorgenannten Ungleichungen folgen. 


3) Nach Gl, (24) entsprieht 9=0 einer Störung der Wellenlänge 21, bei der die Verschiebungen benachbarter 
Wirbel entgegengesetzt gleich sind. Diese Störung wird gelegentlich als „alternative Störung“ bezeichnet; vgl. die 
erste der beiden Arbeiten von Bl. Dolaptschiew, Z. angew. Math. Mech. Bd. 17 (1937), S. 313; Bd. 18 (1938), S. 263. 

%, d. h. beliebiger Wellenlänge. 
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2. Beweis der Ungleichung (49). In die Ausdrücke für A, B und € von Gl. (27) führen 
wir die Größen d und h gemäß Gl]. (37) ein. Dabei benutzen wir die Formeln (39) und die 


Stabilitätsbedingung (41). Wir erhalten zunächst 
| 


50 
cos2a72,=cos da Cosha—iSin’ha | 
und nach einiger Rechnung 
a _ 
2Sin’ha 2 
n? 
(Cosha icosdn) 


Zur Berechnung von p und q bilden wir unter wiederholter Anwendung von (39) und (41) 


die folgenden zusammengesetzten Ausdrücke: 
n? 


| 4 
(a? — 
B—A fi 


Setzen wir diese in die Gln. (46) ein, so bekommt die zu beweisende Ungleichung p - q>0 
die Gestalt 
Um einzusehen, daß (53) erfüllt ist, nehmen wir eine Reihe von Umformungen vor. 
Multiplikation mit 2h*Sin*kr und Einführung der Abkürzungen hn=u,hp=x liefert 


(58). 


4 a 
(5 Sin’u)Cos22=— Sin uCos u (w — x’)Sin’u>0 . (54). 


Unter Benutzung von (39) und der Beziehungen 
Cos 2 x = Sin?’ Cos’ | 
| 
erhalten wir für die einzelnen Bestandteile der linken Seite von (54), und zwar der Reihe 
nach für die Glieder mit u*, «x, u? x’ und «*, 
Sin? u* Sin? x (Cos? — Sin? x), 


(55) 


— 2 u? x Sin u Cos u Sin x Cos x, 
u? Sin? u (Cos’ u +2Sin?’x) und 
— x* Sin* u. 
Durch geeignete Zusammenfassung dieser Bestandteile ergibt sich für die Ungleichung (54) 


a? (u? Cos? x Sin? x — 2 u x Cos Sin x Cos u Sin u + Cos? Sin? 
— (u* Sin* — 2 u? Sin? Sin? u + x* Sin* u) > 0 
oder mit den Abkürzungen 
P=u(xCos u Sin u — u Cos x Sin x) 
57), 


0 Sin? u — u? Sin? x 
kurz 
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Die Ungleichung ist zu beweisen für alle x&-Werte zwischen — « und u, entsprechend dem 
Wertebereich - a<gp<n von @. Da P?’—@? eine gerade Funktion von x ist, genügt es 
jedoch, sich auf den Wertebereich O<x#=<u für x zu beschränken, was im folgenden ge- 
schehen soll. Ä 

Wir überzeugen uns zunächst davon, daß unter dieser Beschränkung P und @ positiv 
sind. Q@>0 ist gleichbedeutend mit 

Sinn „ 
u 
und daher erfüllt. Für P gilt zunächst 
P > uCos u (x Sin u — u Sin x), 

woraus wieder wegen (59) P >0 folgt. 

Wegen P>0, Q>0 genügt es, statt (58) die einfachere Ungleichung 


zu beweisen. Den Beweis führen wir für O<x<u, indem wir zeigen, 
daß R(x) in den Randpunkten <=0 und = verschwindet und da- 
zwischen überall negative Krümmung besitzt (vgl. Bild 6). Wir haben 


(59) 


Bild 6. Hilfsfigur zum R (x) = u (x Cos u Sin u — u Cos x Sin x) — x? Sin? a + u” Sin’ (61) 


bilität 
Stabilitätsbeweis und ersichtlich 


Weiter gilt 

(x) = — 4 u? Cos x Sin x — 2 Sin? u + 2 u? (Cos? -+ Sin? x) 

— 2 u? (Cos x — Sin x)? — 2 Sin? u | RE 

und wegen 

u (Cos x — Sin x) < u<Sin u 
ist schließlich 


womit die negative Krümmung bewiesen ist. 


Es bleibt noch übrig, etwas über die beiden Grenzfälle <=0, =0 und 2@=u,y=n 
zu sagen. Der erste, 9=0, wurde schon in $ 2 behandelt. Eine Störung mit 9—=r bedeutet 
(vgl. Gl. (24)) eine für alle Wirbel einer Reihe gleichmäßige Verschiebung. Einer solchen 
Störung gegenüber ist die Wirbelstraße indifferent; denn durch die Verschiebung entsteht 
einfach eine neue Wirbelstraße mit etwas geänderten Abmessungen. 


$ 4. Zusammenfassung und Schlußbemerkungen. 


In der vorliegenden Arbeit ist erstens gezeigt worden, daß sich der Stabilitätsbeweis 
für die Kärmänsche Wirbelstraße unter Benutzung funktionentheoretischer Methoden führen 
läßt. Zweitens hat sich ergeben, daß die Kärmänsche Wirbelstraße nur einen Sonder- 
fall einer Schar ähnlicher stabiler Wirbelanordnungen darstellt. 

Die andere Frage, ob die eingeführten Wirbelstraßen allgemeinerer Art wirklich auf- 
treten können, wenn ein Körper von einer Flüssigkeit umströmt wird, ließe sich nur experi- 

> mentell entscheiden. Vom theoretischen Standpunkte 


Br _- kann man hierzu sagen, daß aus Symmetriegründen 

= _9°° notwendig ein unsymmetrischer oder schräg 

Eu Fat gestellter Körper erforderlich wäre. Die Wirbel- 

ET hen 9 straße würde sich nicht wie im Kärmänschen Falle 
et gerade, sondern schräg hinter dem Hindernis aus- 


Bild 7. Schematische Darstellung der Ausbildun wilden, 4. werde gegen die Richtung der 
einer Wirbelstraße bei Flüssigkeitsströmung geneigt (vgl. Bild ?). 
suchungen (a. a. O.) besteht neben den Stabilitäts- 
in der Berechnung des Widerstandes, 
den ein durch eine Flüssigkeit bewegter Körper erfährt, in dessen Nachlauf sich eine Wirbel- 
straße ausgebildet hat. Die Rechnungen lassen sich ohne Schwierigkeit auf den von uns 
betrachteten allgemeineren Fall der Wirbelstraße übertragen. Wir gehen hierauf jedoch nicht 
ein, weil die Rechnungen methodisch nichts Neues bieten. Wir erwähnen nur, daß wegen 
der unsymmetrischen Verhältnisse zu der Widerstandskraft eine quer zur Körperbewegung 
gerichtete Auftriebskraft hinzutritt. 130 
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Über den Einfluß der Wandkrümmung auf die Entstehung 
der Turbulenz. 
Von H. Görtler in Göttingen. 
(Aus dem Kaiser-Wilhelm-Institut für Strömungsforschung.) 


1. Einleitung. Die im folgenden darzulegende Untersuchung reiht sich ein unter jene 
Arbeiten, welche die Entstehung der Turbulenz mittels der Methode der kleinen Schwingungen 
untersuchen. Gerade diese Methode hat nach anfänglichen negativen Ergebnissen gute Er- 
folge gezeitigt, als man daranging, die Zähigkeit und die Krümmung des Geschwindigkeits- 
profils der Hauptströmung in richtiger Weise zu berücksichtigen. 


Diese erfolgreichen Untersuchungen setzten im Jahre 1921 ein (Jenaer Vortrag von 
L. Prandtl [3]'), erschienen in demselben Band 1 der ZAMM, in welchem F. Noether [2] 
in einem zusammenfassenden Bericht über das Turbulenzproblem zu dem Schlusse kam, daß 
die Methode der kleinen Schwingungen „zu keinem Ergebnis hinsichtlich der Entstehung der 
Turbulenz geführt“ habe). Die Göttinger Dissertation von O. Tietjens (1922) [4] brachte 
unter Vernachlässigung der Reibung den Nachweis, daß laminare Grenzschichtprofile an ebenen 
Wänden gegenüber kleinen Störungen stabil sind, solange sie konvex sind, dagegen labil, 
wenn sie eine Einbuchtung besitzen. Dabei wurde jedoch nach dem Vorgange von Lord 
Rayleigh [1] das Geschwindigkeitsprofil durch einen Polygonzug ersetzt. Als man in erster 
Näherung die Reibung berücksichtigte, erhielt man das unerwartete und befremdende Resultat, 
daß alle Profile angefachte Schwingungen ergaben. Die Annäherung des stetig gekrümmten 
Profils der Grundströmung durch ein Streckenprofil mußte also als zu starker Eingriff in die 
vorliegenden Verhältnisse betrachtet werden. 


W. Tollmien [5] führte die Rechnungen für die Grenzschichtströmung an der ebenen 
Wand erstmalig unter Berücksichtigung der Krümmung des Geschwindigkeitsprofils durch. 
Er gelangte zu einer kritischen Reynoldsschen Zahl, welche in befriedigender Weise mit 
den Beobachtungen übereinstimmte. Die Störungswellenlängen, welche sich bei der Rechnung 
ergaben, waren jedoch zu groß — sie betrugen etwa die 30fache Verdrängungsdicke der Grenz- 
schicht — als daß man von einer eigentlichen Turbulenz sprechen konnte. Das solche lang- 
welligen Schwingungen aber Anlaß zur Turbulenzentstehung geben konnten, zeigten schon 
die erheblichen Anfachungen, welche Schlichting [8] für das instabile Gebiet berechnete. 
Nun war durch die Tietjensschen Rechnungen der Hinweis gegeben, daß Geschwindigkeits- 
profile mit Wendepunkt — wie sie u.a. etwa durch solche langwelligen Schwingungen ent- 
stehen können — instabil sein werden. Tollmien [10] gelang unter Berücksichtigung der 
Profilkrümmung der Beweis, daß in der Tat solche Wendepunktprofile an ebenen Wänden 
(oder zwischen parallelen ebenen Wänden) instabil sind, und zwar setzt hier eine ganz andere 
wesentlich kräftigere Instabilität als die obige ein mit kurzwelligen Störungen nach Art einer 
echten Turbulenz. Bei diesen letzteren Rechnungen wurde der Reibungseinfluß vernachlässigt. 
An den errechneten kräftigen Instabilitäten dürfte jedoch eine geringe Zähigkeit keine wesent- 
liche Änderung bringen. 


Unter den weiteren Arbeiten — einen zusammenfassenden Bericht gab Schlichting 
1934 [9] — interessiert uns hier noch eine Stabilitätsuntersuchung an gekrümmten Wänden 
(Sehlichting [7]), in welcher für die Strömung im Innern eines rotierenden Kreis- 
zylinders nach denselben Methoden wie an ebenen Wänden kritische Reynoldssche Zahlen 
berechnet wurden. Man hat einen stabilisierenden Einfluß der Wandkrümmung in dem Sinne 
zu erwarten, daß die mit der Grenzschichtdicke gebildete kritische Reynoldssche Zahl an- 
steigt mit zunehmender Wandkrümmung. Diese Erwartung wird begründet durch die Vor- 
stellung, welche man mit dem Einfluß der Zentrifugalkraft verbindet (siehe Prandtl [6)). 
Sie wurde durch die Schlichtingschen Rechnungen bestätigt. 


In der vorliegenden Arbeit sollen Grenzschichtströmungen an schwach gekrümmten 
ruhenden Wänden untersucht werden. Unsere Betrachtungen zielen darauf hin, das dem 
obengenannten Tollmienschen Kriterium der Instabilität einer Grenzschichtströmung an 
einer ebenen Wand entsprechende Kriterium für gekrümmte Wände zu gewinnen. Es werden 
wie bei allen obengenannten Arbeiten nur zweidimensionale Störungen betrachtet. Eine Unter- 
suchung des Einflusses der Wandkrümmung auf die Entstehung der Turbulenz hat praktisches 
Interesse im Hinblick auf die Formgebung von Strömungskörpern. 


1) Die Zahlen in eckigen Klammern verweisen auf das Schriftenverzeichnis am Schluß der Arbeit. 
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2. Die Störungsgleichung. Als ungestörte zweidimensionale Strömung, deren Stabilität 
gegenüber Störungen in Gestalt von kleinen ebenfalls zweidimensionalen Schwingungen unter- 
sucht werden soll, wird eine Grenzschichtströmung an einer schwach gekrümmten Wand vom 
festen Krümmungsradius R angenommen (Bild 1). R wird positiv bei konvexer, negativ bei 

konkaver Wand gesetzt. Mit (y) wird die Geschwindigkeits- 
/ 4 verteilung dieser Strömung dargestellt, wobei y die Koordinate 
senkrecht zur Wand bezeichnet. Veränderungen der Grund- 
strömung U mit der in Strömungsrichtung weisenden Koor- 
dinate « (= Bogenlänge der Wand) werden dabei als vernach- 
lässigbar klein im zugrunde liegenden Bereich von x voraus- 
gesetzt. Ferner wird vorausgesetzt, daß sich U zeitlich sehr 
langsam im Vergleich mit den in Betracht zu ziehenden Stö- 
rungsschwingungen verändern soll. Im übrigen wird das 
Grenzschichtprofil U(y) in üblicher Weise noch insofern ide- 
alisiert, als an Stelle des asymptotischen Übergangs in die 
äußere Strömung ein Anstieg von U=0 an der Wand (y=0) 
bis zu einem Werte U, bei einem bestimmten y= d= „Grenz- 
schichtdicke* angenommen wird. Außerhalb dieser wandnahen 
Schicht nehmen wir — vorwiegend aus später ersichtlichen 
Gründen der rechnerischen Annehmlichkeit (s. Gl. (2,4) und 
(2,5)) und in Übereinstimmung mit Schlichting [7] S. 166 — 


an, daß sich die Strömung wie die starre Rotation U=const a g 


A um den Kreiszylinder verhält; damit ist der Einfluß der Wand- 
PN krümmung auf die äußere Strömung näherungsweise berück- 

sichtigt. 
In nicht zu großen Abständen y> ö ist U wegen des voraussetzungsgemäß großen Wertes 
von R(>y) praktisch konstant. Wir könnten daher die obige Korrektur für eine erste 
Näherung unterlassen und auch im weiteren Verlauf der Rechnung Glieder der Größenordnung 


2 gegen 1 vernachlässigen. Da jedoch der nachzuweisende Einfluß der Wandkrümmung auf 


das Verhalten der zu untersuchenden zweidimensionalen Schwingungen von eben dieser kleinen 
Ordnung ist, wollen wir, um allen Einwänden zu begegnen, grundsätzlich nur Glieder der 
2 
Ordnung streichen. 
Bei der angegebenen Koordinatenwahl lauten die Navier-Stokesschen Gleichungen 
und die Kontinuitätsgleichung in völliger Strenge: 


R du „» R 28» 2R dv 
R du dev 


Darin bezeichnen « und v die Geschwindigkeitskomponenten in x- bzw. y-Richtung, p den 
Druck, o die Dichte und » die kinematische Zähigkeit. Wir haben uns ferner auf y„=0 zu 
beschränken. Schließlich wird in (2,1) durch A der Differentialoperator 


bezeichnet. 
Der Kontinuität wird durch Einführung der Stromfunktion Y vermöge 
genügt. N 
Der oben bezeichneten Grundströmung, welche die Stromfunktion y, 


haben möge, überlagern wir nun eine kleine zweidimensionale Schwingung mit der Strom- 
funktion y(xz,9,t). Da U=U(y) selbst eine Lösung der hydrodynamischen Gleichungen sein 
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soll — wobei wir dann anzunehmen haben, daß die zeitliche Veränderung von U so langsam 
gegen die zu betrachtenden kleinen Schwingungen stattfindet, daß wir sie vernachlässigen 


können —, erhält man nach Einsetzen von ?=y,-+ y in (2,1) und Elimination des Druckes 
die bezüglich y und ihren Ableitungen linearisierte Störungsgleichung 


Man kann jede Störungsschwingung y (x, y,t) aus einer Reihe von Partialstörungen in 
Fourierscher Weise aufbauen. Wegen der Linearität der Störungsgleichung kann man sich 
daher auf die Untersuchung eines Störungsgliedes der Gestalt 


. . 2.2... 3) 


beschränken. Es wird a reell angenommen und 7, wo 4 die Wellenlänge der Störung be- 


zeichnet. ß ist im allgemeinen komplex; der Realteil 5, ist die Kreisfrequenz der Schwin- 
gung, während der Imaginärteil 5; ein Maß für die Anfachung oder Dämpfung der Störung — 
je nachdem es größer oder kleiner als Null ist — darstellt. Schließlich ist der Realteil c, 
von ce die Phasengeschwindigkeit der Störung. 


Durch Einsetzen von (2,3) in (2,2) ergibt sich: 


--715 +57 +2 +5) (- 
Führen wir hier die Funktion F(y) Den 
also F=const für y>ö, mu EN entwickeln R+y' 
(By binomisch, so ergibt sich unter Vernachlässigung von Gliedern mit dem 


Faktor 5 die Störungsgleichung’) 


[73 1 ’ 2y 3 
(2,5). 


Unter Berücksichtigung der Randbedingungen, welche das Verschwinden der Störungs- 
komponenten an der Wand und den glatten Anschluß in y=ö an die nach y=» abklingende 
Lösung außerhalb der Grenzschicht zu fordern haben, liegt ein Eigenwertproblem vor: Bei 
vorgegebener Geschwindigkeitsverteilung U(y) und Zähigkeit » oder auch Reynoldsscher 
Zahl der Grundströmung ist zu einer bestimmten Störungswellenlänge (Parameter a) der Wert 
von e und damit der Anfachung bzw. Dämpfung /; zu ermitteln. 

Für Flüssigkeiten kleiner Reibung kann man die rechte Seite von (2,5) vernachlässigen. 
Im Hinblick auf das Turbulenzproblem interessieren nur höhere Reynoldssche Zahlen. 
Wir werden daher in der vorliegenden Untersuchung — im Einklang mit Tollmien [10], an 
dessen Ausführungen sich unsere Betrachtungen unmittelbar anschließen — die Reibung ver- 
nachlässigen und sie nur dort, wo sie die Eindeutigkeit einer zu treffenden Bestimmung in 
der reibungslosen Rechnung entscheidet, heranziehen (siehe unten, Abschnitt 4). 

Unseren Betrachtungen liegt somit endgültig die Störungsgleichung 


2 


2) Mit etwas anderer Koordinatenwabl siehe auch Schliehting [7] S. 166. Im weiteren gehen aber die 
Schliehtingschen;Vernachlässigungen über die unseren hinaus, da man es dort mit kleinen «Werten, also mit lang- 
welligen Störungen zu tun hat. 


XUM 
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zugrunde. Wir werden später mehrfach die durch 


definierte Hilfsfunktion y heranziehen, für welche nach (2,6) die Differentialgleichung 
F'+-— F' 
21 R 


besteht. 

Zur Störungsdifferentialgleichung (2,6) bzw. (2,8) sei hier noch bemerkt, daß sie im Gegen- 
satz zu (2,5) nur von der 2. Ordnung ist. Es können daher nur zwei Randbedingungen er- 
füllt werden. Wir kommen im Abschnitt 4 darauf zurück. 


3. Das Lösungssystem und die Randbedingungen im reibungslosen Fall. Die Differential- 
gleichung (2,6) besitzt eine singuläre Stelle der Bestimmtheit dort, wo F=c wird, wenn diese 
Stelle innerhalb der Grenzschicht liegt. Daß letzteres in der Tat für neutrale Lösungen f; =0 
an ebenen Wänden der Fall ist, ist ein auf Rayleigh zurückgehender Satz (siehe hierzu 
den von Tollmien [10], S. 9) angegebenen kurzen Beweis). Für schwach gekrümmte Wände 
ändert sich daran nichts. 

Die konvergenten Reihenentwicklungen eines Fundamentalsystems der Differential- 
gleichung (2,6) um diese kritische Stelle, die mit y=4y. bezeichnet werden möge, sind in be- 
kannter Weise herzustellen und werden durch 


9,=(y— Ye) 2 a,(y — Ye), 


(3,1) 


gegeben. In (3,1) ist «,—=b,=1 und der Index c weist auf den Wert der betreffenden Funktion 
an der Stelle „=. hin. Es werde hier ferner festgelegt, daß In (y— 9.) für positive reelle 
Y—Ye reell sei. 

Zur Herstellung der Randbedingung bei y=ö benötigen wir die Lösung der Differential- 
gleichung (2,6) bzw. (2,8) außerhalb der Grenzschicht, von der wir zu fordern haben, daß sie 
nach außen abklingt. Betrachten wir dazu die letztere Differentialgleichung, welche für „= ö 
die Gestalt . 


annimmt, so ist eine Näherungslösung®), welche der Differentialgleichung bis auf Glieder der 


2 
Ordnung genügt, solange noch y<R ist, durch 


ay? 
x = Const e= 1148, +55 


gegeben. Man hat somit folgende Randbedingungen an die Lösung innerhalb der Grenzschicht 
zu stellen: 


(3,2). 
bzw. 9 (0)=0, 


Die jeweils letztere Randbedingung gewährleistet den stetigen Anschluß mit stetiger Tangente 
bei y=Ö. 


4. Der Einfluß der Reibung. Bei der Darlegung des Reibungseinflusses können wir uns 
kurz fassen, da bereits Schlichting [7], S. 171ff. in Übertragung der Tollmienschen Be- 
trachtung [5], S.26ff. (siehe auch [10] S. 83ff.) an ebenen Wänden auf die Strömung in rotie- 


3) Man erhält sie, indem man zunächst unter Vernachlässigung des Gliedes «a? * X löst und dann dureh Iteration 


einmal verbessert. Übrigens kann man auch die exakte Lösung angeben, welche sich mittels einer Hankelschen 
Funktion darstellen läßt. 
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renden Zylindern die notwendigen Feststellungen gemacht hat, deren Gültigkeit sich unmittel- 
bar auf unseren Fall übertragen läßt. 

Man kann sich auf die Berücksichtigung der Reibungswirkung in einem kleinen Streifen 

1 

von der Größenordnung (a Re F’.) ? beschränken (Re= 2), Neben der den Erwartungen 
entsprechend sich ergebenden Abflachung der Störung — 9, in (3,1) bleibt dann in y= y. endlich, 
so daß die Anwendung der Methode der kleinen Schwingungen ihre Rechtmäßigkeit bewahrt — 
ergibt sich ein Phasensprung in den Störungen, der auch im Grenzfall verschwindender Rei- 
bung bestehen bleibt. Es ergibt sich nämlich für neutrale Schwingungen, daß unter den 
Zweigen des Logarithmus in (3,1) allein die Fortsetzung auf dem Wege von positiv reellen 
Y— Ye durch die untere komplexe Halbebene zu negativen y— y. physikalisch möglich ist. 
Für positive y— y. ist 


3 
und somit für negative y— Ye ; 
3 


Dieselbe Übergangssubstitution gilt auch für angefachte Schwingungen (ec; > 0), während für 
gedämpfte Schwingungen — ir durch +ir zu ersetzen ist. 

Es sei noch bemerkt, daß man bei Berücksichtigung der Reibung neben den in dem 
genannten und bei großen Reynoldsschen Zahlen sehr engen Streifen um y= y. der Grenz- 
schicht korrigierten Lösungen (3,1) zwei weitere linear unabhängige Lösungen für das Fun- 
damentalsystem der vollständigen Differentialgleichung 4. Ordnung erhält. Sie sind ebenfalls 
von Tollmien und Schlichting a. a. O. untersucht worden. Wir können uns daher auf 
wenige Bemerkungen beschränken. Während die eine der beiden hinzukommenden Lösungen 
mit wachsenden y-Werten sehr rasch über alle Grenzen wächst und daher ausscheidet, klingt 
die andere sehr rasch ab in einer Grenzschicht zweiter Ordnung in der gegebenen Grenz- 
schicht. Der Einfluß der Wandkrümmung auf diese beiden Lösungen ist sehr gering im Ver- 
gleich zu dem Einfluß auf 9, und 9,. Unter Hinzunahme der rasch abklingenden Lösung 
könnte man noch die weitere Randbedingung, daß auch die Tangentialkomponente der Störung 
an der Wand verschwindet, erfüllen. Die genannte Lösung hat aber ihrem ganzen Charakter 
gemäß nur geringen Einfluß auf den Störungsverlauf. Wir können daher für hinreichend 
große Reynoldssche Zahlen auf eine über die obengenannte Übergangssubstitution hinaus- 
gehende Berücksichtigung der Reibung verziehten. Es entfällt damit natürlich die Möglichkeit 
der Bestimmung einer kritischen Reynoldsschen Zahl; letzteres wird auch nicht die Auf- 
gabe dieser Untersuchung sein. Wir verfahren damit in Übereinstimmung mit den Tollmien- 
schen Untersuchungen zum Beweise des einleitend genannten Stabilitätskriteriums mit dem 
Ziele, den Einfluß der Wandkrümmung auf die Stabilität einer vorgegebenen Grenzschicht- 
strömung festzustellen. 


5. Eine notwendige Bedingung für die Existenz angefachter 
Schwingungen. Bereits Lord Rayleigh [la], S. 486, hat für 
den Fall der Strömung an ebenen Wänden nachgewiesen, daß 
angefachte Störungen nur dann existieren können, wenn das Grenz- 
schichtprofil einen Wendepunkt y=y, hat (Bild 2), d. h. es muß 
notwendig fürein y=y, die Größe U” (y,) verschwinden (siehe 
hierzu auch Tollmien [10], S. 88£.). 

In Verallgemeinerung dieser Bedingung beweisen wir für 
den Fall schwach gekrümmter Wände folgenden Satz: Notwendig 
für die Existenz angefachter Schwingungen ist ein 
Vorzeichenwechsel von 


Zum Beweise ziehen wir die Hilfsfunktion y heran. Mit 


3 


(5,2), 
2y _ 1 - 


AUM 
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worin der Querstrich die konjugiert komplexen Größen bezeichnet, ist für komplexe Werte 


von e — es befindet sich dann keine Singularität auf der reellen Achse der komplexen 
y-Ebene — unter Berücksichtigung der homogenen Randbedingungen 
"IF-eP .... . GB. 


Es sind also in der Tat angefachte (und überhaupt nichtneutrale) Schwingungen (ce; +0) nur 


dann möglich, wenn r+3r das Vorzeichen wechselt, womit wegen F'+2, F' 


= ltr or) für große R der behauptete Satz bewiesen ist. Wir bezeichnen im 
weiteren jene Stelle, an welcher der Ausdruck (5,1) verschwindet, mit y=9, und versehen 
die zugehörigen Werte der auftretenden Funktionen mit dem Index 0. 
An konvexen Wänden (R >0) ist 
Y >Ys (Y = Wendepunkt des Profils), 
an konkaven Wänden (R<0) dagegen 
ist 9, < Y,, wie aus der Diskussion des 
Ausdruckes (5,1) wegen U”>0O für 
0<ysy und U”<O für „<sysö 
folgt (Bild 3). Hierbei nehmen wir an, 
daß das Geschwindigkeitsprofil einen 
Wendepunkt besitzt. Es folgt aber auch, 
daß an konvexen Wänden bereits ein y, 
ohne Auftreten eines 4, vorhanden sein 
kann, an konkaven Wänden dagegen 
nicht. 
Bild 3. Ur (yo) + Grenzschichtprofle mit Wende- 
punkt sind nach der Grenzschichttheorie 
dadurch gekennzeichnet, daß in der betreffenden Strömung ein Anstieg des von der äußeren 
Strömung der Grenzschicht eingeprägten Druckes p in Strömungsrichtung stattfindet. Es ist 
Auf Grund unserer obigen Betrachtungen können also angefachte Schwingungen an 
konvexen Wänden eventuell schon vor dem Druckminimum, an konkaven jedoch auf jeden 
Fall erst hinter dem Druckminimum auftreten. Der Nachweis, daß sie tatsächlich in den 
genannten Situationen auftreten, d. h. daß unsere obige Bedingung des Verschwindens von 


U” +3 U’ auch hinreichend für die Existenz angefachter Eigenschwingungen ist, wird die 


Aufgabe der weiteren Untersuchung sein. 

Wir erinnern hier nochmals daran, daß bei unseren Rechnungen die Reibung vernach- 
lässigt wird. Wie man von den Untersuchungen spezieller Geschwindigkeitsprofile an ebenen 
Wänden weiß, können dort manche Profile ohne Wendepunkt durch Anbringung einer Korrektur 
für die innere Reibung instabil werden. Es handelt sich aber bei den resultierenden An- 
fachungen auf Grund dieser kleinen Reibungskorrektur für große Reynoldssche Zahlen um 
eine wesentlich kleinere Größenordnung, als bei den für Wendepunktprofile in „reibungsloser“ 
Rechnung von Tollmien aufgefundenen kräftigen Anfachungen. Wenn wir also hier von 
einer notwendigen und hinreichenden Bedingung für die Existenz angefachter Schwingungen 
sprechen, so nur im Sinne der Existenz dieser kräftigen Instabilitäten und unter Außeracht- 
lassung der obengenannten möglichen kleinen Anfachungen, die durch Berücksichtigung der 
Reibung für große Reynoldssche Zahlen bei speziellen unserer Bedingung nicht genügenden 
Profilen zu erhalten sind. 


6. Nachweis, daß ein Vorzeichenwechsel von U’ + H U’ auch hinreichend für die Existenz 


angefachter Eigenschwingungen ist. Wir nehmen an, es liege ein Geschwindigkeitsprofil U (4) 


vor, für das in einem Wandabstand „=, ein Vorzeichenwechsel von U” + R U’ stattfindet: 


U’, + 4 U’,=0. Wir bestimmen für dieses Profil zunächst eine neutrale Eigenlösung unseres 
Eigenwertproblems =0). 
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Dazu gehen wir aus von der Lösung 


unserer Differentialgleichung, welche wir hier am einfachsten in der Gestalt der nullgesetzten 
linken Seite von (2,4) heranziehen; bei unserer Näherung bezüglich R und bei den Parameter- 
werten e=F,, lautet sie: 


Die Lösung (6,1) genügt jedoch nicht den Randbedingungen. 


Wesentlich ist nun folgende Feststellung: Die zu (6,1) gehörige Hilfsfunktion z: 


hat wegen 


x = eul2R (v" + vr) 
2 


(bei Vernachlässigung von Gliedern = U 5 wie vereinbart) in y=,, einen Wendepunkt: 


Die von der neutralen Eigenlösung x, zu erfüllende Differentialgleichung lautet 
. 3 
mit den Randbedingungen 
0) 1 
F” 
Im Falle e=F, ist der in der Differentialgleichung auftretende Ausdruck u ee a für alle 


o 
y-Werte in 0<y=ö eines Vorzeichens, nämlich negativ. Eine Singularität liegt nicht vor. 
Gehen wir daher von der Lösung (6,2) mit a=0,c=F, aus und lassen in der Differential- 
gleichung «a anwachsen, so ändert sich die Lösung mit «a stetig und zwar so, daß bei Fest- 
halten des Wertes y(ö) das Abfallen der Lösung nach dem Inneren der Grenzschicht bis zur 
Wand immer langsamer erfolgt (Bild 4). Hinzu kommt der im gleichen Sinne wirkende Einfluß 
des Parameters «a in der Randbedingung bei y=d. Zwar genügt die Ausgangslösung (6,2) bei 
y=ö nicht der vorgegebenen Randbedingung mit a=0, jedoch ist die Abweichung wegen 


3 
des voraussetzungsgemäß großen Wertes von R gering — es ist für (6,2) 2’ (6) — SR 16) =. 


Rückt man die Lösung zunächst so zurecht, daß sie mit 
a=0 die Randbedingung bei y=ö erfüllt, und läßt dann 
a wachsen unter Einhaltung der Randbedingung bei y=Ö, 
so muß der Nullpunkt der zugehörigen Lösung immer 


A weiter nach der Wand zu wandern. Jene Lösung, für 

X welche der Knotenpunkt in y=0 fällt, erfüllt dann beide 
F Randbedingungen, ist also eine Eigenlösung y=y,„, und 

zwar mit c=c„—=F,. Dieser Existenzbeweis ist jeden- 


falls richtig, wenn > ist wegen des obengenannten 


Fehlers in der Randbedingung, da sonst das Zurecht- 
y rücken der Lösung zur Erfüllung der Randbedingung bei 


Ö a=0 bereits die Wanderung des Knotenpunktes durch 
14 y=0 erwirken könnte. Nun geht aber unsere neutrale 


7 Eigenlösung für R— © in die von Tollmien [10], S. 104 f. 


angegebene neutrale Eigenlösung an der ebenen Wand 
über. Bei dem von Tollmien angegebenen Beispiel [10], 
S. 110f. ist «,ö=1,12, und auch für endliche aber große 


Bild 4. R wird es sich demnach um ay-Werte >; handeln. 
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Wir wollen nun prüfen, ob in der Nachbarschaft der nachgewiesenen neutralen Eigen- 
lösung angefachte Eigenschwingungen existieren. Wir schlagen damit auch hier den von 
Tollmien für R=» beschrittenen Weg ein und werden uns eng an seine Rechnungen an- 
zuschließen haben. Mit der Feststellung angefachter Schwingungen für ein Profil mit Vor- 


zeichenwechsel von U” + 5 U’ ist dann der Beweis erbracht, daß dieser Vorzeichenwechsel für 
die Existenz angefachter Schwingungen hinreichend ist. 


Der Nachweis selbst ist etwas langwierig. Da es sich vielfach um eine wörtliche 
Wiederholung der Tollmienschen Rechnungen [10] handelt, werden wir uns kürzer fassen 
können, indem wir an gegebener Stelle einfach auf seine Ausführungen verweisen und uns 
somit lediglich auf eine Darlegung des Gedankenganges unter Berücksichtigung der Wand- 
krümmung beschränken. 


Es sei x, eine zu unserer oben nachgewiesenen neutralen Eigenlösung y„ mit c„=F, 
benachbarte Lösung, welche die Randbedingung y;(0)=0 erfüllt. Dann ist die Bedingung 


mit (6,4) 


2y 1 1 
notwendig und hinreichend dafür, daß y, eine Eigenlösung ist, d. h., daß y; auch der Rand- 
bedingung bei y=ö genügt. Dabei gehören die Parameter «a„, c„—=F, zur neutralen Eigen- 


lösung x„, die Parameter a, ce zur Lösung x,. Der Beweis ist wie bei Tollmien [10] zu 
erbringen. 


Setzt man 
so ist nach (6,4) 
d 


de e)(F— Andy+ Aa dy= -5) (0) (6) (6,5). 


Entwickelt man zunächst y; und y„ um die singuläre Stelle y=y.(F=e), wobei so normiert 
werde, daß %„(9,)=1, so wird wegen y.=Y, 
c 
27, 3 
c 
(6,6) 
de 


(analog Tollmien [10], S. 102), wo k und k, durch die Bedingungen an der Wand bestimmt 
werden und uns nicht weiter interessieren. Der Index ce verweist auf den Wert der Funktion 


bei y= Ye. 
Man bestimmt nun zunächst den größten Beitrag des ersten Integrales von (6,5) in einem 
|4de| 
Intervall —e<y—y,=<e, wo &® Fr | aber e> also etwa unter Ver- 


wendung der Reihenentwicklungen (6,6) und unter Berücksichtigung der Ü bergangssub- 
stitution (4,1). Wegen In(—&)=1In (es) — ri wird dieser Beitrag gleich 


3 
(F +7) 


Nun ist, wenn man um 9, entwickelt, wegen 
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Unsere Integralbedingung liefert also: 


mit . (6,8). 
0 +: 
Durch Spaltung in Real- und Imaginärteil ergibt sich für E+0 


75) 4er (6,9) 


und 


da 2y 
610) 
0 


Für hinreichend kleine |4«a?’| und |4e| kann man hierin y; in erster Näherung durch y„ er- 
setzen. Wegen F"+ er =0 für yZy, ist der aus der Entwicklung (6,7) hervorgegangene 


Ausdruck WR negativ. Ist also etwa E<[0 (wie in einem von Tollmien für 4=0 

angegebenen Beispiel), so ist c; positiv, wenn was durch Aa? widerspruchslos 

erfüllbar ist. Ebenfalls ist für E >0 ein positives ec; zu erzielen, wenn Der>0, also wieder 

da®<0. Auch für E=0, wo Ac,=0 wird, erhält man aus (6,8) für ec; einen Ausdruck, der 
für da’ <O positiv wird. 

Damit ist nachgewiesen, daß es für Geschwindigkeitsprofile U(y) mit einem 


Vorzeichenwechsel von +5 U’stets angefachte Eigenschwingungen gibt. 


Nach den Ausführungen am Ende des Abschnitts 5 werden also Geschwindigkeitsprofile an 
konvexen Wänden vor dem Minimum des von der äußeren Strömung eingeprägten Druckes, 
an konkaven Wänden erst hinter diesem instabil. Jedoch ist der Einfluß der Wandkrümmung 


1 
außerordentlich gering, denn wegen U-U (145). Die Stelle y,, wo 


U’+ 5 U’ verschwindet, wird demnach nur wenig von dem Wendepunkt , des Profils ab- 


weichen. 

Die Tatsache jedoch, daß konvexe Wände einen im obigen Sinne anfachenden, konkave 
dagegen einen stabilisierenden Einfluß ausüben, dürfte zunächst im Hinblick auf die Vor- 
stellungen, welche man mit der Wirkung der Zentrifugalkräfte verbindet, überraschen; 
siebe hierüber L. Prandtl [6]. Die Zentrifugalkräfte spielen offenbar keine ins Gewicht 
fallende Rolle. 

Unsere obigen Betrachtungen gelten auch unverändert für Grenzschichten an rotierenden 
Zylindern. Hier wäre die Wirkung der Zentrifugalkraft eine ganz andere, wenn sie in Er- 
scheinung treten würde. Ruht die Flüssigkeit außerhalb der Grenzschicht, so nimmt U von 
einem größten Wert U, an der Wand auf U=0 für y=dab. Die Diskussion der Gleichung 


+ U’ (y)=0 (siehe Abschnitt 5) gibt also wiederum das Resultat: 9, >y, an kon- 


vexen, 4,< 4%; an konkaven Wänden. Erst nach Auftreten eines Wendepunkts bei Grenz- 
schichten von Flüssigkeiten innerhalb eines rotierenden Zylinders tritt die berechnete zwei- 
dimensionale Instabilität ein, und entsprechend bereits vor Auftreten eines Wendepunkts bei 
Grenzschichten an der äußeren Wand eines rotierenden Zylinders. 

Um nochmals zu den obigen Rechnungen zurückzukehren, sei ergänzend folgendes 
bemerkt. Tollmien hat auch in der Nachbarschaft einer weiteren neutralen Eigenlösung, 
nämlich der ausgearteten Eigenschwingung = U, a=c=0 im Falle von Wendepunktprofilen 
angefachte Eigenschwingungen erhalten. Wir verzichten hier darauf, die analogen Betrach- 
tungen für gekrümmte Wände durchzuführen. Einmal interessiert der Fall so kleiner a-Werte 


physikalisch nicht (aö< o(7) ) da man es dann nicht mit einer eigentlichen Turbulenz zu 
tun hat; ferner versagt die Betrachtung des Abschnitts 4, da die Größe (a Re F’.)- "s für 
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sehr kleine a-Werte nicht mehr klein bleibt. Schließlich ist zu bemerken, daß im Falle 
a R= O(1) die Randbedingung (3,2) bei y= ö zu verbessern wäre, da in der Differentialgleichung 


für y außerhalb der Grenzschicht ein vernachlässigtes Glied = einen Einfluß von einer von 


uns vereinbarungsgemäß zu berücksichtigenden Ordnung auf die Lösung und damit auf die 
genannte Randbedingung gewinnt. Diese müßte lauten: 


1 


Wir wollen uns damit begnügen hier darauf hinzuweisen, daß e; zwischen einer in der Nach- 
barschaft der ebengenannten (@a=0) und unserer oben nachgewiesenen neutralen Eigenlösung 
(mit endlichem a) positiver Werte fähig ist. 


7. Zusammenfassung. Es ist von praktisch wesentlicher Bedeutung, etwa im Hinblick 
auf die Gestaltung von Tragflügelprofilen, den Einfluß der Flächenkrümmung auf die Ent- 
stehung der Turbulenz zu erforschen. In der vorliegenden Arbeit wird das zweidimensionale 
Problem der Stabilität laminarer Grenzschichtprofile an schwach gekrümmten Wänden vom 
konstanten Krümmungsradius R gegenüber zweidimensionalen Störungen untersucht. 


Das Tollmiensche Kriterium der Instabilität an ebenen Wänden: „Geschwindigkeits- 
profile mit Wendepunkt sind instabil“ wird für gekrümmte Wände abgewandelt in: „Ge- 


schwindigkeitsprofile U (y) mit Vorzeichenwechsel von U’ sind instabil.“ Demnach 


tritt die Instabilität gegenüber zweidimensionalen Störungen an konkaven Wänden erst 
eine gewisse Strecke hinter dem Minimum des von außen der Grenzschicht eingeprägten 
Druckes, bei konvexen Wänden bereits vor demselben auf. Der Einfluß der Wandkrümmung 
ist jedoch außerordentlich gering. Die gefährlichen Störungen halten sich bei schwacher 
Wandkrümmung in der engen Nachbarschaft der von Tollmien für ebene Wände berechneten. 
(Es handelt sich also um Wellenlängen, die durchaus einer echten Turbulenz entsprechen.) 


Zum Schluß weisen wir darauf hin, daß in einer demnächst erscheinenden zweiten Arbeit 
der Nachweis einer dreidimensionalen Instabilität an konkaven Wänden von der Art 
der Taylorschen Instabilität zwischen rotierenden Zylindern (G. I. Taylor, Phil. Trans. 
Roy. Soc. London, Series A, Bd. 223 (1923), S. 289) erbracht werden wird; gegenüber dieser ist 
der oben nachgewiesene schwach stabilisierende Einfluß konkaver Wände auf zweidimensionale 
Störungen von untergeordneter Bedeutung. 


Die Anregung zu der vorliegenden Untersuchung verdanke ich meinem hochverehrten 
Chef und Lehrer, Herrn Prof. L. Prandtl. 
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Beitrag zur Theorie des Tragflügels in schwach inhomogener 
Parallelströmung. 
Von F. Vandrey in Göttingen. 


\ A Tir betrachten im folgenden einen Tragflügel in einer Parallelströmung, deren Richtung 

zur Flügelachse senkrecht und deren Geschwindigkeit längs jeder Stromlinie konstant 
ist, von einer Stromlinie zur anderen aber verschieden sein darf. Die Grundströmung soll 
also quellenfrei, nicht aber wirbelfrei sein. Zunächst nehmen wir an, daß die Wirbel der 
Grundströmung nur einzelne Schichten bilden, die Gebiete konstanter Geschwindigkeit der 
Grundströmung voneinander trennen. Von diesen Gebieten seien im ganzen n vorhanden, die 
wir uns irgendwie numeriert denken. Wir ersetzen in üblicher Weise den Flügel durch 
einen tragenden Wirbelfaden und das System ‚er abgehenden Wirbel, die wir in Richtung 
der Grundströmung annehmen. Das Koordinatensystem legen wir so, daß die x-Achse der 
Strömung entgegenzeigt, die y-Achse mit dem tragenden Wirbel zusammenfällt und in Flug- 
richtung gesehen nach rechts positiv ist und die z-Achse senkrecht zu den beiden anderen 
und nach unten gerichtet ist. Wir betrachten einen Schnitt durch die Strömung sehr weit 
hinter dem Flügel. In diesem Schnitt ist außer der Grundströmung nur noch eine ebene 
Querströmung vorhanden, die in jedem der » Gebiete wirbelfrei sein und daher ein Potential 
besitzen wird. Die Berechnung dieser ebenen Querströmung soll unsere erste Aufgabe sein. 
Aus ihr ergibt sich dann die Zirkulation des Tragflügels in bekannter Weise als der Potential- 
sprung am abgehenden Wirbelband. 

Wir stellen zunächst die Bedingungen zusammen, denen die » Potentialfunktionen an 
den Wirbelflächen der Grundströmung und am abgehenden Wirbelband des Flügels genügen 
müssen. 

An einer Unstetigkeitsfläche der Grundströmung, in der etwa die beiden Gebiete @, und @. 
aneinander stoßen, muß zunächst der Druck auf beiden Seiten derselbe sein, da von der Fläche 
keine Kräfte aufgenommen werden können, ferner muß die Fläche eine Stromfläche sein. Be- 
zeichnet man die Geschwindigkeit der Grundströmung in @, mit V,, das Potential der Quer- 
strömung mit ®,, und entsprechend in @., so läßt sich die Druckbedingung, wie C. Koning') 
gezeigt hat, bis auf quadratische Glieder in den vom Flügel herrührenden Zusatzgeschwindig- 
keiten erfüllen durch die Forderung 


Bu ander Grenze . . (1). 
Die zweite Bedingung läßt sich durch die Forderung erfüllen, daß sich die Normal- 
geschwindigkeiten der Querströmung wie die Geschwindigkeiten der Hauptstiömung verhalten 


sollen 


Zu diesen Bedingungen kommen noch zwei Bedingungen am Wirbelband hinzu. Erstens muß 
die Normalgeschwindigkeit dort stetig sein: 


Zweitens muß dort die Prandtlsche Tragflügelgleichung gelten, die wir in der Form schreiben 
2 


Wir haben jetzt also die folgende Aufgabe der ebenen Potentialtheorie zu lösen: 


In jedem der » Gebiete @, ist eine Potentialfunktion ®, so zu bestimmen, daß an den 
Grenzen die Bedingungen (1) und (2) erfüllt sind und in den Gebieten, die von dem abgehenden 
Wirbelband getroffen werden, außerdem die Bedingungen (3) und (4) gelten. 


Diese Aufgabe ist natürlich in voller Allgemeinheit nicht zu lösen, wir wollen uns 
daher auf den Fall beschränken, daß die Geschwindigkeiten der Grund- 
strömung nur wenig von einem Mittelwert abweichen. Wir setzen 


1) C. Koning: Influence of the Propeller on other Parts of the Airplane Structure. In: Durand, Aerodynamie 
Theory. Berlin 1935, Bd. 4 S. 380 bis 384 und 402 bis 406. 
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Die Zu der Druckbedingung (1) legt es nahe, für das Potential ®, dann den Ansatz zu 
machen 


wo ® eine einheitliche Potentialfunktion in der ganzen Ebene bedeutet. Dieser Ansatz 
erfüllt die Druckbedingung bis auf quadratische, die Kontinuitätsbedingung bis auf lineare 
Glieder in den Zusatzgeschwindigkeiten. Der Einfluß dieser geringeren Berücksichtigung der 
Kontinuitätsbedingung soll unten noch näher untersucht werden. Zur Berechnung von ® kann 
dann die Gleichung (4) dienen. Wir dividieren sie noch durch (1 —v,) und erhalten bis auf 
Glieder höherer Ordnung 


> 


Dies ist aber die Prandtlsche Gleichung für einen Ersatzflügel mit dem geometrischen 
Anstellwinkel a,(1+2v,) in einer gleichförmigen Parallelströmung von der Geschwindigkeit V. 
Sie kann daher nach den hierfür bekannten Verfahren gelöst werden. Ist ® (oder auch nur 
T=(®,,— ®P,):= 0) berechnet, so wird die Zirkulationsverteilung des Flügels in der inhomogenen 
Strömung auf Grund des Ansatzes 


T,=(1—ev,)T. 


Die Zirkulationsverteilung wird also an den Stellen unstetig, wo der Flügel die Unstetigkeits- 
flächen der Grundströmung schneidet. Die Auftriebsverteilung des Flügels in der inhomogenen 
Grundströmung wird dann nach dem Satze von Kutta und Joukowski 


dA 

dy 

bis auf quadratische Glieder. Sie ist also identisch mit der Auftriebsverteilung des Ersatzflügels. 
Die erweiterte Tragflügelgleichung (5) gestattet auch noch den Grenzübergang zu einer 

Grundströmung mit stetig veränderlicher Geschwindigkeit. Ist die Geschwindigkeit am Orte 

des Flügels etwa gegeben durch 


so erhält man die Auftriebsverteilung des Flügels als Auftriebsverteilung eines Ersatzflügels 
mit der c, t-Verteilung des gegebenen und der Anstellwinkelverteilung a,(1+2v(y)) in 
gleichförmiger Parallelströmung von der Geschwindigkeit V. 

Wir wollen jetzt noch auf die Fehler eingehen, die der Rechnung durch die Vernach- 
lässigung der linearen Glieder in der Kontinuitätsbedingung anhaften. Diese Vernachlässigung 
hat u. a. zur Folge, daß nur die Ungleichförmigkeiten der Grundströmung am Orte des Flügels 
in die Rechnung eingehen, während alle den Flügel nicht schneidenden Wirbelflächen der 
Grundströmung ohne Einfluß bleiben. Wir denken uns nun bei einer gegebenen Verteilung 
der Grundgeschwindigkeit am Orte des Flügels die Berechnung der zum Ersatzflügel gehörenden 
Potentialfunktion ® durchgeführt und fragen danach, ob sich die Unstetigkeitsflächen der 
Grundströmung von ihren Schnittpunkten mit dem Flügel aus nicht so fortsetzen lassen, daß 
für die dann entstehende neue Grundströmung die Kontinuitätsbedingung an den Gebietsgrenzen 
von selbst erfüllt ist. Man sieht leicht ein, daß dies möglich ist, denn die Kontinuitäts- 


bedingung ist exakt erfüllt, falls auf allen Grenzen #0 ist, d. h. falls die Grenzen der 


Gebiete Stromlinien der durch ® gegebenen ebenen Querströmung hinter dem Ersatzflügel 
sind. Der Verlauf dieser Stromlinien läßt sich aber bestimmen, man vermag daher zu be- 
urteilen, ob die wirklich vorliegende Grundströmung mit der durch ® gegebenen in der Nähe 
des Flügels so weit übereinstimmt, daß die Berechnung der Auftriebsverteilung nach der oben 
abgeleiteten Tragflügelgleichung zulässig erscheint. 

Bei der Bestimmung derjenigen Grundströmung, für die die Berechnung der Auftriebs- 
verteilung eigentlich gilt, kann noch eine grundsätzliche Schwierigkeit auftreten, falls der 
Abwind des Ersatzflügels ee i längs der Spannweite sein Vorzeichen wechselt. In 
diesem Falle gibt es nämlich Stromlinien der ebenen Querströmung, die oberhalb des Flügels 
von den Aufwindgebieten zu den Abwindgebieten und unterhalb in entgegengesetzter Richtung 
verlaufen. Da diese den Flügel zweimal schneidenden Stromlinien in der inhomogenen Grund- 
strömung Linien konstanter Geschwindigkeit sein sollten, müßte die Grundgeschwindigkeit 
an ihren beiden Schnittpunkten mit dem Flügel dieselbe sein. Dies wird aber im allgemeinen 
nicht der Fall sein, und dann gibt es genau genommen überhaupt keine Grundströmung, für 
die die angewendete Berechnungsmethode richtig ist. Ein Beispiel für diesen Fall läßt sich 
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leicht angeben, wenn man etwa für einen unendlich langen unverwundenen Flügel konstanter 
Tiefe sich das Potential ® der ebenen Querströmung hinter dem Ersatzflügel vorgibt und 
dann nach der Tragflügelgleichung die zugehörige Geschwindigkeitsverteilung der inhomogenen 
Grundströmung am Flügel ausrechnet. Wir betrachten als einfaches Beispiel dieser Art das 
Potential 


2k y’+l 
9,=-7 Va 1—  (aretg -aretg! z* 
Hierin ist 


und k ein als klein vorausgesetzter Proportionalitätsfaktor für die Übergeschwindigkeiten der 
zu berechnenden inhomogenen Grundströmung. In Bild 1 ist ein Stromlinienbild der durch ® 
gegebenen Querströmung und darunter die Abwindverteilung des Ersatzflügels gezeichnet. 
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Bild 1. Beispiel für das Stromlinienbild der ebenen Quer- Bild 2. Verteilung der zusätzlichen Geschwindigkeit der 
strömung hinter einem unendlich langen Flügel konstanter inhomogenen Grundströmung, der zusätzlichen Zirkulation 
Tiefe und Verteilung des Abwindes längs der Spannweite. und des zusätzlichen Auftriebes für das Beispiel des Bildes 1. 


Nach diesem Stromlinienbild müßte in den beiden mit A und B bezeichneten Punkten des 
Flügels in der inhomogenen Grundströmung die gleiche Geschwindigkeit herrschen. Die nach 
der Tragflügelgleichung berechnete Geschwindigkeitsverteilung der inhomogenen Grundströmung 
ist in Bild 2 oben gezeichnet; man erkennt, daß die Geschwindigkeiten in A und B 
nicht übereinstimmen. Es gibt hier also eigentlich überhaupt keine Grundströmung, für die 
die berechneten Zirkulations- und Auftriebsverteilungen des Flügels in der inhomogenen 
Grundströmung richtig. sind. Dies hindert natürlich nicht, daß sie trotzdem ein qualitativ 
brauchbares Bild der in inhomogenen Strömungen vorliegenden Verhältnisse geben können. 
Ähnlich wie in diesem Falle erhalten übrigens bei dem unendlich langen unverwundenen 
Flügel konstanter Tiefe auch in allen anderen Fällen die Gebiete mit kleinerer Grund- 
geschwindigkeit Aufwind, da für den verwundenen Ersatzflügel die Ungleichmäßigkeiten des 
Anstellwinkels in der Zirkulationsverteilung weitgehend ausgeglichen sind. Für den unver- 
wundenen unendlich langen Flügel konstanter Tiefe gibt es daher genau genommen niemals 
eine inhomogene Grundströmung, für die seine Auftriebsverteilung nach der hier dargelegten 
Methode berechnet werden darf. Da aber der Fall eines solchen Flügels nur ein vereinfachtes 
Bild ist für die Verhältnisse im Mittelstück eines endlichen Flügels großer Spannweite, so 
braucht man dieser grundsätzlichen Schwierigkeit keine allzu große Bedeutung beizumessen. 

Man kann weiterhin noch die Frage stellen, ob nicht für gewisse Grundströmungen auf 
die Voraussetzung der Kleinheit der Übergeschwindigkeiten v, verzichtet werden darf. (Aus 
physikalischen Gründen muß hierbei v, > —1 gefordert werden, da der Fall v„<—1 einer 
Anströmung des Flügels von hinten entspricht, bei der andere Verhältnisse vorliegen. Der 
Fall v„=—1 (d. h. Anströmgeschwindigkeit Null in einem der Gebiete) erfordert eine be- 
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sondere Untersuchung, er soll deshalb hier ausgeschlossen werden.) Dies ist möglich, wenn 
man dem Ansatz für die Potentiale ®, die Form gibt 


1+v, 
‘ Für die Potentialfunktion ®* erhält man damit als Bedingung 
1 (a 2 
Pk 


Sind dann die Gebietsgrenzen (bzw. bei stetiger Geschwindigkeitsverteilung die Linien 
konstanter Grundgeschwindigkeit) Stromlinien der durch ®* gegebenen Strömung, so ist auch 
die Kontinuitätsbedingung exakt erfüllt, und man erhält für die Auftriebsverteilung eines 
Flügels in einer solchen Grundströmung wiederum das einfache Ergebnis 

dia 

Zum Schluß unserer Fehlerbetrachtungen sei noch auf einen Umstand hingewiesen, der 
bei der Anwendung auf praktische Fälle beachtet werden muß. Unserer Rechnung liegt die 
Annahme zugrunde, daß die Grundströmung durch den Einfluß des Flügels nicht wesentlich 
geändert wird. Bei größeren Auftriebskräften trifft aber diese Annahme nicht mehr zu, und 
es können Verformungen der Unstetigkeitsflächen der Grundströmung auftreten, die zu Ab- 
weichungen der gemessenen Auftriebsverteilungen von den berechneten führen. Es sei hierzu 
verwiesen auf die Messungen von J. Stüper zur Prüfung der Koningschen Theorie, wo 
derartige Einflüsse vermutet wurden. 

Als Beispiel zu den obigen Entwicklungen 
ist in Bild 3 die Auftriebsverteilung eines dca 
Rechteckflügels vom Seitenverhältnis b/t =5 aa 
gezeichnet, der einen Schraubenstrahl von der | 


Breite 0,2 b mit 10° und 20 °/ Übergeschwin- 
digkeit schneidet. Darunter ist zum Vergleich 
noch die Auftriebsverteilung eingezeichnet, 
wie sie sich ohne Beachtung der Druckbe- 
dingungen durch formales Rechnen mit der 
gewöhnlichen Prandtlschen Tragflügelglei- 
chung ergeben würde. In Bild 4 ist dann 
noch die eigentliche Gestalt des „Schrauben- 
strahles“ gezeichnet, wie sie sich aus den 
Fehlerbetrachtungen ergibt. Man sieht, daß 
diese sehr erheblich von der gewünschten 
Form eines Kreisstrahles abweicht, und man 
wird daher einige Bedenken dagegen haben, 
die Rechnung auf den Fall eines Flügels im 
Kreisstrahle anzuwenden. 
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Bild 3 (oben). Auftriebsverteilung eines Yr=(trv)V 
benstrahl von der Breite 0,2b mit 100%, 
und 200, Übergeschwindigkeit schneidet. 
Darunter zum Vergleich die ohue Berück- 
siehtigung der Druckbedingungen an den 
Strahlrändern berechnete A uftriebsver- f Öebiel I 6ebiet I \ 
teilung des gleichen Flügels. H N # \ 2 
\ \ / 
Bild 4 (rechts). Form der Strahlränder, N 
für die die in Bild 3 gezeichneten Auf- N 
triebsverteilungen in inhomogener Strö- 1 


Ein Vergleich der berechneten Verteilung mit Messungsergebnissen ist z. Zt. leider noch 
nicht möglich. Messungen zur Prüfung der Koningschen Theorie, die in der hier darge- 
legten enthalten ist, sind von J. Stüper?) ausgeführt worden. Sie ergaben beträchtlich 
kleineren Einfluß als die Rechnung (30 °/o der berechneten Werte). Stüper führt dies neben 


J. Stüper: Luftf.-Forschg. Bd. 15 (1938) ,S. 181. 
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Einflüssen der Grenzschicht darauf zurück, daß bei den Messungen die Flügeltiefe größer 
war als der Strahldurchmesser und daß deshalb die Strahlform durch den Flügel erheblich 
gestört wurde, so daß die Voraussetzung der Rechnung nicht mehr zutrifft, daß der Einfluß 
des Flügels auf die Grundströmung zu vernachlässigen ist. Nach unseren Ergebnissen könnte 
man noch an eine andere Möglichkeit zur Erklärung des Unterschiedes zwischen Rechnung 
und Messung denken. Da sich nämlich die Koningsche Rechnung genau genommen gar 
nicht auf einen Kreisstrahl bezieht, sondern in dem bei den Messungen vorliegenden Fall 
(Rechteckflügel mit Endscheiben) auf eine Strahlform, die der in Bild 4 gezeichneten ähnlich 
sein dürfte (falls es überhaupt eine alle Bedingungen erfüllende Grundströmung in diesem 
Falle gibt), so könnte es sein, daß der Kreisstrahl nicht mehr eine zulässige Grundströmung 
ist. Es ist hierbei zu erwarten, daß für den Kreisstrahl der Abwind größer und die Zirkulation 
kleiner wird als für die der Rechnung eigentlich zugrunde liegende Strömung, so daß die 
Unterschiede auch hierdurch erklärt werden könnten. Nach beiden Erklärungsversuchen 
müßte man jedoch erwarten, daß die Übereinstimmung von Rechnung und Messung sich 
bessert, wenn die Flügeltiefe klein gegen den Strahldurchmesser gemacht wird. 


Zusammenfassung. Die Berechnung der Auftriebsverteilung eines Tragflügels in einer 
schwach inhomogenen Parallelströmung mit von Ort zu Ort verschiedener Geschwindigkeit 
wird zurückgeführt auf die Berechnung der Auftriebsverteilung eines Flügels gleicher Tiefen- 
verteilung mit einer zusätzlichen Verwindung in homogener Parallelströmung. Die Verwindung 
wird dabei bestimmt aus der Bedingung der Druckgleichheit auf den beiden Seiten einer 
Fläche konstanter Grundgeschwindigkeit. Voraussetzung für die Anwendbarkeit der Methode 
ist, daß sich die Geschwindigkeit der Grundströmung in der zum Flügel senkrechten Richtung 
nicht viel ändert, genau genommen müssen die Linien konstanter Grundgeschwindigkeit in 
einem vertikalen Schnitt weit hinter dem Flügel Stromlinien der ebenen Querströmung hinter 
dem Ersatzflügel sein. Als Beispiel wird ein Rechteckflügel im Schraubenstrahl behandelt. 
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Ebene Bewegung eines Wirbelkranzes am rotierenden 
radialen Schaufelstern von endlichen Schaufellängen’. 


Von Fr. Florin in Berlin-Tegel. 
(Aus dem Institut für Mechanik an der Technischen Hochschule Berlin.) 


1. Einleitung. Die Bewegung eines Wirbels innerhalb eines von einer idealen Flüssig- 
keit erfüllten rotierenden Kreissektors, bei überlagerter Quellströmung, wurde bereits von 
Kucharski?) untersucht. Die radialen Begrenzungen dieses Sektors erstrecken sich vom Dreh- 
punkt bis ins Unendliche und begrenzen eine nach außen gleichfalls unbegrenzte Flüssigkeit. 

Man erhält nun eine bessere Annäherung an die beispielsweise bei einer Pumpe vor- 
liegenden Verhältnisse, wenn die Wirbelbahnen unter Zugrundelegung des Betriebszustandes 
einer Pumpe gleich für einen ganzen Stern radialer, sich nicht mehr bis ins Unendliche er- 
streckender, aber von einer unbegrenzten Flüssigkeit eingeschlossener Schaufeln untersucht 
werden. Ein solcher radialer Schaufelstern ist das einfachste und deswegen einer exakten 
Berechnung der Wirbelbahnen am besten zugängliche Schema eines zweidimensionalen 
Schaufelrades. Dasselbe kann dem Aussehen eines Schaufelrades dadurch noch weiter an- 
gepaßt werden, daß man die radialen Schaufelspuren nicht nur auf einem Kreis vom endlichen 
Außenradius o„ endigen, sondern auch auf einem zweiten, konzentrischen Kreis vom Innen- 
radius o; ihren Anfang nehmen läßt. Ein dem wirklichen Aussehen eines Schaufelrades noch 
besser sich anpassendes Schema erhält man dadurch, daß man den Schaufelspuren eine nach 
logarithmischen Spiralen gekrümmte Form gibt. Auch in diesen Fällen lassen sich die 
Wirbelbahnen noch mit exakten Methoden berechnen. 

Bislang sind aber in diesen Fällen lediglich die Strömungserscheinungen ohne Einzel- 
wirbel untersucht worden. Als erste Arbeit dieser Art erschien 1918 eine Monographie von 
Kucharski?’), in welcher u.a. für einen radialen Schaufelstern, dessen Schaufeln vom Dreh- 
punkt bis zu einem Kreis vom endlichen Radius o„ reichen, das Bild der relativen Strom- 
linien und das Förderhöhenverhältnis in Abhängigkeit von der Schaufelzahl » ermittelt wurde. 


1) Diese Arbeit wurde von der Fakultät für Maschinenwesen der Technischen Hochschule Berlin als Doktor- 
arbeit angenommen. Referenten: Prof. W. Kucharski und Prof. G. Hamel. Vorsitzender: Prof. H. Föttinger. 
Herrn Prof. Kucharski schuldet der Verfasser für wertvolle Ratschläge und Belehrungen besonderen Dank. 

2) W. Kucharski: Bewegung eines Wirbels in einem nach außen offenen Kreissektor. Forsch. Ing.-Wes. Bd.8 
(1937), 8.14. Ferner eine bisher unveröffentlichte Arbeit über den gleichen Gegenstand, in welcher die Integration der 
Wirbelbahnen geschlossen durchgeführt wurde. 

3) W. Kueharski: Strömungen einer reibungsfreien Flüssigkeit bei Rotation fester Körper. München und 
Berlin 1918. 
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Eine Erweiterung auf den radialen Schaufelstern, dessen Schaufeln von den konzentrischen 
Kreisen mit den Radien o; und o„ begrenzt werden, brachte die 1925 erschienene Arbeit von 
Spannhake*). In derselben wurde das Förderhöhenverhältnis in Abhängigkeit von der 


Schaufelzahl » und dem Radienverhältnis a“ berechnet. Dieser Arbeit folgten nun in kür- 


za 
zeren Zeitabständen die Arbeiten von Sörensen’°) (1927), Schulz*) (1928) und Busemann'’) 
(1928), die eine Krümmung der Schaufelspuren nach logarithmischen Spiralen in Betracht 
ziehen. In ihnen wurde gleichfalls das Förderhöhenverhältnis berechnet. 

In der vorliegenden Arbeit, die auf Anregung von Herrn Prof. Kucharski im 
Institut für Mechanik der Technischen Hochschule Berlin durchgeführt wurde, sollen nun die 
Bewegung und die Ruhelagen eines Wirbelkranzes am radialen Schaufelstern von endlichen 
Schaufellängen untersucht werden, und zwar unter Nachahmung des Betriebszustandes einer 
Kreiselpumpe. Es ist bekannt, daß beim Betrieb einer Pumpe an den Schaufeln Wirbel ent- 
stehen, die sich von den Schaufelenden periodisch ablösen. Auf Veranlassung von Herrn Prof. 
Spannhake sind diese Wirbelbildungen und Wirbelablösungen an umlaufenden radialen 
und logarithmisch-spiraligen Schaufeln, mit und ohne Überlagerung einer Quellströmung, 
zum Gegenstand einer experimentellen Untersuchung gemacht worden, die von Frietsch‘®) 
durchgeführt wurde. 

Da sich die folgenden Untersuchungen auf die Bewegung von Wirbeln in einer reibungs- 
freien Flüssigkeit beziehen, in der die Entstehung von Wirbeln unmöglich ist, müssen wir 
die Wirbel als von vornherein vorhanden ansehen. Die Untersuchungen von Frietsch 
können also nur insoweit zum Vergleich herangezogen werden, als dieselben sich nicht mehr 
auf die eigentliche Entstehung der Wirbel beziehen. Hierzu wird man aber erst noch die 
angekündigte weitere Arbeit abwarten müssen, welche die theoretische Verarbeitung der Ver- 
suchsergebnisse bringen soll. 

Die in der vorliegenden Arbeit in Anlehnung an Spannhake und Schulz benutzte 
konforme Abbildung des Schaufelsterns auf den Umfang eines Kreises bringt es mit sich, daß 
wir gleich die Bewegung eines ganzen Wirbelkranzes studieren müssen, der aus so vielen 
Wirbeln von gleichem Drehsinn und gleicher Stärke besteht, als das Kreiselrad Schaufeln besitzt. 

Die Einführung eines solchen Wirbelkranzes dürfte aber auch vom physikalischen Stand- 
punkt aus durchaus gerechtfertigt sein. Denn die Symmetrie in der Anordnung der Schaufeln 
spricht doch dafür, daß in den einzelnen Schaufelräumen einander entsprechende Vorgänge 
sich abspielen, daß also an den einander entsprechenden Stellen gleich drehende und gleich 
starke Wirbel vorhanden sind, und daß dieselben kongruente Bahnen durchlaufen. Die von 
Frietsch veröffentlichen relativen Stromlinienbilder geben für diese Auffassung eine aus- 
gezeichnete Bestätigung. 

In der Literatur findet man bereits eine große Zahl von Abhandlungen, in denen 
Wirbelbahnen für einfachere Fälle berechnet wurden. Die Verfasser dieser Arbeiten bedienen 
sich entweder des Helmholtz-Kirchhoffschen Grundsatzes, so z.B. die uns besonders 
interessierenden Arbeiten von Hamel®), Kucharski'°) und Bartlh''), oder sie bedienen sich 
der im folgenden ebenfalls benutzten und im wesentlichen als bekannt vorausgesetzten 
Routhschen Methode, wie z. B. die Arbeiten von Paul’, Kneschke') und Kneschke- 
Matthes'*) Zwischen diesen beiden Verfahren besteht natürlich ein analytischer Zusammen- 
hang, über den Herr Prof. Kucharski im Winter-Semester 1937/38 in einem Vortrag im 
Seminar für Mechanik an der Technischen Hochschule Berlin gesprochen hat. 

Über die Theorie der Wirbelbewegung findet man in der Literatur gleichfalls eine 
größere Zahl von Abhandlungen, auf deren Aufzählung hier verzichtet werden kann. Es mag 
genügen, auf eine neuere Arbeit von Kneschke') über diesen Gegenstand hinzuweisen, 
zumal derselben weitere Literaturangaben entnommen werden können. 


4) W. Spannuhake: Die Leistungsaufnahme einer parallelkränzigen Zentrifugalpumpe mit radialen Schaufeln. 
Festschrift zur Hundertjahrfeier der T. H. Karlsruhe. Karlsruhe 1925. 

5) E. Sörensen: Potentialströmungen durch rotierende Kreiselräder. Z. angew. Math. Mech. Bd.7 (1927), 8. 80. 

6) W. Schulz: Das Förderhöhenverhältnis radialer Kreiselpumpen mit logarithmisch -spiraligen Schaufeln. 
Z. aungew. Math. Mech. Bd. 8 (1928), S. 10. 

A. Busemann: Das Förderhöhenverhältnis radialer Kreiselpumpen mit logarithmisch -spiraligen Schaufeln. 
Z. angew. Math. Mech. Bd. 8 (1928), S. 372. 

°), E.Frietsch: Wirbelbildung u. Kräftewirkung an umlaufenden Kreiselradschaufeln. VDI-Forseh.-Heft 384 (1937). 

#») G. Hamel: Bewegung eines geradlinigen Wirbels um eine Buhne. Z. angew. Matlı. Mech. Bd. 13 (19939), S. 98. 

10), W. Kucharski: a.a. 0.?). 

11) W. Barth: Wirbelbahnen um Wände und Platten von unendlich kleiner Wandstärke. Z. augew. Math. Mech. 
Bd. 10 (1930), S. 247. 

12) E.Paul: Bewegung eines Wirbels in geradlinig begrenzten Gebieten. Z. angew. Math. Mech. Bd. 14 (1934), 8. 105. 

” Fe A. Kneschke: Über die Bewegung von Wirbeln in einem einseitig begrenzten Kanal. Ann. Phys. Bd. 14 

(1932), S. 655. 

14) A. Kneschke u. 8. Matthes: Wirbelbewegung um einen Kreiszylinder. Ann, Phys. Bd. 9 (1931), 8. 916. 

15) A. Kneschke: Zur Theorie der Wirbelbewegung. Z. angew. Math. Mech. Bd. 18 (1938), S. 343, 
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2. Die Abbildungsfunktion. Die von Spannhake, Sörensen, Schulz und Busemann 
benutzten Abbildungsfunktionen, durch die der endliche Schaufelstern entweder auf den Um- 
fang eines Kreises oder auf die reelle Achse abgebildet wird, bleiben auch bei der Berechnung 
von Wirbelbahnen ohne weiteres anwendbar, da der Wirbel am Kreis und an der unend- 
lichen Geraden gespiegelt werden kann. Wir werden uns im folgenden der von Schulz be- 
nutzten Abbildung des nach logarithmischen Spiralen gekrümmten Schaufelsterns auf den 
Einheitskreis bedienen, indem wir dieselbe auf den zu untersuchenden Spezialfall des radialen 
Schaufelsterns anwenden. 

Der Einheitskreis sei in der Ebene mit der komplexen Veränderlichen z=e x +iy=r-eiv 
gelegen. Dann überführen ihn die folgenden beiden Transformationen in einen in der Ebene 
mit der komplexen Veränderlichen w=u+iv=o-ei” gelegenen, aus n radialen Schaufeln 
bestehenden Schaufelstern: 


Durch die erste Transformation wird der Einheitskreis in der z-Ebene, dessen Mittelpunkt 
mit dem Punkt z=0 zusammenfällt, auf das von p —q bis p-+gq reichende, auf der reellen 
Achse der 4-Ebene gelegene Geradenstück abgebildet. Durch die zweite Transformation wird 
eine Verkleinerung aller den Koordinatenanfangspunkt als Scheitelpunkt besitzenden Winkel 
auf den n-ten Teil, sowie eine radiale Verzerrung bewirkt. Dadurch entsteht der gewünschte 
radiale Schaufelstern, bestehend aus » vom Innenradius 9; =yYp—q bis zum Außenradius 
0a=Yp-+gq reichenden Schaufeln. In Bild 1 ist dieser Schaufelstern für den Fall n=6 
dargestellt worden. 

Fassen wir die beiden Transformationen (1) und (2) zu einer einzigen zusammen, so 


erhalten wir die Beziehung 


durch welche der Einheitskreis der z-Ebene in den radialen 
Schaufelstern der w-Ebene überführt wird. Dabei wird der 
Außenbereich des Kreises bereits auf die gesamte w-Ebene 
abgebildet. 

Setzt man in Gl. (3) p=g, so erhält man den Spezial- 
fail, daß die radialen Schaufeln bis zum Drehpunkt des 
Schaufelrades reichen. 

Singularitäten treten auf durch die z—4- Abbildung 
in den Endpunkten der Schaufeln, die den beiden Kreis- 
punkten z= + 1 entsprechen, und durch die A— w- Ab- 
bildung im Punkt w,=0, der dem für p>g außerhalb 
des Einheitskreises, und für p=g auf dem Einheitskreis 
gelegenen reellen Punkt 


1 


Bild 1. entspricht. 


3. Das komplexe Potential des Wirbelkranzes.. An einer außerhalb des Einheitskreises 
der z-Ebene gelegenen Stelle 2,=r,: eifı möge sich ein linksdrehender Wirbel von der Zirku- 


lation /', also von der Wirbelstärke "=;; befinden. Durch die Abbildung (3) werden aus 


diesem einen Wirbel » Wirbel von gleichem Drehsinn und gleicher Stärke, die miteinander 
den in Bild 1 eingezeichneten Wirbelkranz bilden. 

Es soll nun das komplexe Potential des durch diesen Wirbelkranz am radialen Schaufel- 
stern in der unbegrenzten Flüssigkeit erzeugten momentanen Stromlinienfeldes ermittelt werden. 
Zu dem Zweck brauchen wir nur das komplexe Potential desjenigen momentanen Strom- 
linienfeldes anzuschreiben, das durch den einen an der Stelle z, der 2-Ebene gelegenen Wirbel 
um die Kontur des Einheitskreises in der unbegrenzten Flüssigkeit erzeugt wird. In diesem 
Stromlinienfeld muß der Einheitskreis eine Stromlinie sein. Das komplexe Potential dieses 
Feldes läßt sich in bekannter Weise dadurch gewinnen, daß der an der Stelle z, gelegene 


AUM 
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Wirbel am Einheitskreis gespiegelt wird, wodurch an der Stelle = ein spiegelbildlicher Wirbel 
1 

von gleicher Stärke, aber von entgegengesetztem Drehsinn entsteht. Dies hat aber zur Folge, 
daß wir in dem durch diese beiden Wirbel erzeugten Stromlinienfeld eine rechtsdrehende 
Zirkulation von der Stärke m um den Einheitskreis erhalten. Da wir jedoch später zwecks 
Berücksichtigung der Zirkulationsströmung um die Schaufeln ohnehin eine Zirkulation um den 
Einheitskreis in Ansatz bringen müssen, und da wir über ihre Stärke noch in ganz bestimmter 
Weise verfügen wollen, ist es zweckmäßig, diese hier unerwünschte Zirkulation unter Wahrung 
der Randbedingung, daß der Einheitskreis Stromlinie sein muß, dadurch zum Verschwinden 
zu bringen, daß wir im Punkt z=0 noch einen dritten, gleich starken, linksdrehenden Wirbel 
anbringen. Das komplexe Potential F(z) dieser drei Wirbel lautet dann 


Dieser Ausdruck stellt zugleich das in den Koordinaten der 2-Ebene ausgedrückte kom- 
plexe Potential des Wirbelkranzes am Schaufelstern der w-Ebene dar. 


4. Die Routhsche Stromfunktion für die Wirbelbahnen am Einheitskreis. Um die Routhsche 
Stromfunktion'*) für die von unserem Wirbel in der z2-Ebene des Einheitskreises innerhalb 
der unbegrenzten Flüssigkeit durchlaufenen Wirbelbahnen zu erhalten, berechnen wir zunächst 
nach dem Helmholtz-Kirchhoffschen Grundsatz die Geschwindigkeitskomponenten r, und 
r,ö, des an der Stelle 2, der z-Ebene befindlichen Wirbels. Hierzu ist eine im vorliegenden, 
sich nur auf die z-Ebene beziehenden Fall einfache Limesbildung vorzunehmen, als deren 
Ergebnis wir zunächst die konjugierte Geschwindigkeit 2,—=&,—ij, des Wirbels erhalten. 
Vor der Grenzwertbildung z>z, ist nach 2 zu differenzieren; diese Differentiation ist aber 
nicht an dem durch Gl. (5) dargestellten komplexen Potential F(z), sondern an dem um das 


Potential —-In (2—z,) des betrachteten Wirbels, der sich ja selbst keine Geschwindigkeit 
induziert, verminderten Potential durchzuführen. Wir erhalten so 


d 1 


2, 

Hieraus ermitteln wir die Radialgeschwindigkeit r, und die Umfangsgeschwindigkeit r,g, des 

Wirbels zu 

m 


Diese Geschwindigkeitskomponenten erfüllen zwar die Kontinuitätsgleichung, nicht aber 
die Bedingung der Wirbelfreiheit. Wir können also die Geschwindigkeitskomponenten des 
Wirbels als die Geschwindigkeitskomponenten einer zwar quellenfreien, aber nicht wirbel- 
freien Hilfsströmung deuten, für die zwar eine Stromfunktion, aber kein Geschwindigkeits- 
potential, und somit auch kein komplexes Potential existiert. 

Wir wollen nun die zu diesen Geschwindigkeitskomponenten gehörende Routhsche 
Stromfunktion y. für die Wirbelbahnen in der z-Ebene berechnen. Dies hat mittels der 
bekannten Beziehungen 


z . 


zu geschehen. Die gesuchte Routhsche Stromfunktion y. ergibt sich hieraus nach Einsetzen 
der Geschwindigkeitskomponenten (7) bis auf eine additive Konstante zu 
m 
Die Gleichung der in der z-Ebene durchlaufenen Wirbelbahnen lautet dann y,=konst. Diese 
Bahnen sind Kreise, die zum Einheitskreis konzentrisch sind. 


5. Die Routhsche Stromfunktion für die Wirbelbahnen am Schaufelstern. Der Routh- 
schen Stromfunktion für die Wirbelbahnen am Schaufelstern muß ein besonderer Abschnitt 
gewidmet werden, da die im vorigen Abschnitt ermittelte Routhsche Stromfunktion y, für 
die Wirbelbahnen am Einheitskreis der z-Ebene nicht zugleich die lediglich in den Koordinaten 
der z-Ebene ausgedrückte Routhsche Stromfunktion zu für die Wirbelbahnen am Schaufel- 


16, Für die in dieser Arbeit untersuchten Wirbelbahnen existiert eine sog. Routhsche Stromfunktion, die zu 
der Schar der Wirbelbahnen in derselben Beziehung steht, wie die gewöhnliche Stromfunktion zu der Sehar der 
zugehörigen Stromlinien. 
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stern der w-Ebene darstellt. Die in der z-Ebene ermittelten Wirbelbahnen gehen nämlich bei 
‘ der konformen Abbildung in die w-Ebene nicht in die gesuchten Wirbelbahnen der w-Ebene 
über. Vielmehr müssen wir uns zu ihrer Ermittlung des Routhschen Verfahrens bedienen, 
das den Zusammenhang zwischen der bekannten Stromfunktion y, und der gesuchten Strom- 
funktion zw klarstellt. Hiernach besteht zwischen diesen beiden Stromfunktionen die wichtige 
Beziehung 


Mit Hilfe dieser Beziehung gestaltet sich die Berechnung von 4. folgendermaßen: Wir 

bilden zunächst aus unserer Abbildungsfunktion (3) den Differentialquotienten I, ‚ indem wir 
1 


ihn nach der Kettenregel zerlegen in 


dh, 
dw, dw, di, da, 
dz, di, de di ' al) 
dw, 
und indem wir die im Zähler und Nenner stehenden Differentialquotienten = und dm, us 
1 1 
den Gl. (1) und (2) errechnen und in Gl. (11) einsetzen. Dadurch erhalten wir 
1 
(12 
dz, n . . . . . . . . . . ). 
Der absolute Betrag dieses Differentialquotienten ergibt sich hieraus zu 


Setzen wir diesen Ausdruck, sowie den für y- erhaltenen Ausdruck (9) in (10) ein, so erhalten wir 
m 


Ze=5 In (14). 


+ 


An diesem Ausdruck für die gesuchte Stromfunktion y,„ stört noch, daß er bis jetzt 
nicht einheitlich entweder in den Koordinaten der z-Ebene, oder in den Koordinaten der w- 
Ebene ausgedrückt ist. Wir erhalten nun den bei weitem einfacheren Ausdruck, wenn wir 
die Stromfunktion y.„ einheitlich in den Koordinaten der z-Ebene ausdrücken. 

Zu dem Zweck müssen wir die in Gl. (14) noch vorkommende Koordinate o der w-Ebene 
durch die Koordinaten r und p der z-Ebene ausdrücken. Dies gelingt am einfachsten dadurch, 
daß wir in unserer Abbildungsfunktion (3) Polarkoordinaten einführen. Dadurch erhalten wir 


Hieraus erhalten wir durch Trennung von Reellem und Imaginärem die beiden Gleichungen 


N (16). 
o" -cos(n -coSp 


Durch Quadrieren und Addieren dieser beiden Gleichungen erhalten wir den gewünschten 


Für spätere Zwecke sei noch vermerkt, daß wir durch Division der beiden Gl. (16) für 9 den 
Ausdruck 


+ 


1 

(r+ 


erhalten, der es uns gestattet, auch 9 in den Polarkoordinaten der 2-Ebene auszudrücken. 


(18) 
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Durch Konstantsetzen des Ausdruckes für y.„. (Gl. 14) erhalten wir dann unter Berück- 
sichtigung von Gl. (17) die in den Polarkoordinaten der z-Ebene ausgedrückte Gleichung der 
Wirbelbahnen in der mw-Ebene des Schaufelsterns. 

Um das erhaltene Resultat zu kontrollieren, wurde der Helmholtz-Kirchhoffsche 
Grundsatz auch noch zur unmittelbaren Berechnung der konjugierten Geschwindigkeit der 
Wirbel am Schaufelstern der w-Ebene benutzt. Hieraus konnte dann nach längeren Zwischen- 
rechnungen die Differentialgleichung für die Wirbelbahnen in der w-Ebene, ausgedrückt in 
den Koordinaten der z-Ebene, ermittelt werden. Da die Integration dieser komplizierten 
Differentialgleichung auf zu große Schwierigkeiten stieß, wurde die Kontrollrechnung in der 
Weise zu Ende geführt, daß die soeben erhaltene Gleichung der Wirbelbahnen implizite 
differenziert wurde. 


6. Die Verdrängungsströmung, Vorbemerkungen. Bisher haben wir die Gleichung der 
Wirbelbahnen am Schaufelstern der w-Ebene für den Fall des Fehlens einer Grundströmung 
abgeleitet. Es fehlt aber noch die Berücksichtigung derjenigen Strömung, die dem Betriebs- 
zustand einer Pumpe entspricht. Diese Strömung kann aus einzelnen Teilströmungen zu- 
sammengesetzt werden, für die sieh in der Literatur die Bezeichnungen Verdrängungs- 
strömung, Durchflußströmung und Zirkulationsströmung eingebürgert haben. Für 
diese Teilströmungen sollen in den folgenden Abschnitten die in den Koordinaten der z-Ebene 
ausgedrückten Stromfunktionen ermittelt werden. Wir erhalten dann die Gleichung der Wirbel- 
bahnen in der w-Ebene beim Vorhandensein dieser Teilströmungen dadurch, daß wir diese 
Stromfunktionen zu der gleichfalls in den Koordinaten der z2-Ebene ausgedrückten Routh- 
schen Stromfunktion 4, (Gl. 14) addieren und ihre Summe gleich einer Konstanten setzen. 

Als erste der bereits von vornherein in der Pumpe vorhandenen Teilströmungen be- 
handeln wir die Verdrängungsströmung, die durch die Rotation des Laufrades in der 
unendlichen Flüssigkeit hervorgerufen wird. Dabei ist zwischen der wirbelfreien, abso- 
luten, nicht stationären, und der nicht wirbelfreien, relativen, stationären 
Verdrängungsströmung zu unterscheiden. Durch Addition der Stromfunktion der relativ- 
stationären Verdrängungsströmung zu der Routhschen Stromfunktion y„ (Gl. 14) erhalten 
wir die Routhsche Stromfunktion für die relativen Wirbelbahnen am Schaufelstern beim 
Vorhandensein der Verdrängungsströmung. Dagegen erhalten wir durch Addition der Strom- 
funktion der absoluten Verdrängungsströmung zu der Routhschen Stromfunktion y,. nicht 
etwa die Routhsche Stromfunktion der absoluten Wirbelbahnen, sondern lediglich die Routh- 
sche Stromfunktion für das momentane Richtungsfeld der Absolutgeschwindigkeiten, mit denen 
sich die Wirbel an den verschiedenen Stellen des Feldes gerade weiterbewegen. 

Die Ermittlung der während der Drehung des Schaufelrades durchlaufenen absoluten 
Wirbelbahnen ist wesentlich schwieriger, als die Ermittlung der relativen Wirbelbahnen ; 
denn sie erfordert auch ein Studium des zeitlichen Ablaufs der Bewegungen auf den Wirbel- 
bahnen. Wir werden uns daher in den folgenden Untersuchungen auf die Ermittlung der 
relativen Wirbelbahnen beschränken. 

Die Ermittlung der absoluten Verdrängungsströmung ist im Sinne der Potentialtheorie 
ein Randwertproblem zweiter Art, dem die Randbedingung entspricht, daß an jeder Stelle 
der Schaufeloberfläche die senkrecht zur Kontur der Schaufel genommene Komponente der 
absoluten Strömungsgeschwindigkeit mit der in der gleichen Richtung genommenen Kompo- 
nente der Umfangsgeschwindigkeit des betreffenden Punktes der Schaufeloberfläche überein- 
stimmt. 


7. Die Stromfunktion für die Verdrängungsströmung. Wir setzen zunächst nach dem 
Vorbild von Spannhake, von dem der Gedankengang der folgenden Rechnung stammt, das 
komplexe Potential F-,,. (2) der absoluten Verdrängungsströmung in der z-Ebene des Einheits- 
kreises für dessen Mittelpunkt als Laurent- Entwicklung an. Da nur das Kreisäußere dar- 
gestellt werden soll, kommen nur die fallenden Potenzen in Betracht. Wir setzen also 


Fr. „(J)=i- > 


Diese Reihe konvergiert unbedingt und gleichmäßig. 

Die in der Reihe auftretenden Koeffizienten a; müssen nun so bestimmt werden, daß 
die Randbedingungen der Aufgabe erfüllt werden, daß also an jeder Stelle der radialen 
Schaufeln, die entgegengesetzt dem Uhrzeiger mit der Winkelgeschwindigkeit » rotieren 
sollen, ihre senkrecht zum Radius gerichtete Umfangsgesehwindigkeit mit der in die gleiche 
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Richtung fallenden Komponenten der dort herrschenden Strömungsgeschwindigkeit überein- 
stimmt. Hieraus ergibt sich nach einigen Zwischenrechnungen die Bedingungsgleichung 
1:0, sinp 
n 
(g:eosp+p) * 
In dieser Gleichung haben wir jetzt die Bestimmungsgleichung für die in der Laurent- 
Reihe (19) auftretenden Koeffizienten a; gewonnen. Ihre Bestimmung erfolgt durch harmo- 
nische Analyse der auf ihrer linken Seite stehenden periodischen Funktion, die sich gemäß 
dieser Gleichung in eine Fouriersche Reihe entwickeln läßt. Da diese Funktion eine un- 
gerade Funktion ist, treten in der Fourierschen Reihe nur Sinusglieder auf. 
Nach Ausführung der harmonischen Analyse und einer für die graphisch-numerische 
Auswertung geeigneten Umformung erhalten wir für die a; 


k-ar-sn(k9) . 2.2. (MO). 


3 q 2 


Damit ist das durch Gl. (19) dargestellte komplexe Potential der absoluten Verdrängungs- 
strömung bekannt. Durch Aufspalten in Real- und Imaginärteil erhalten wir dann hieraus 
die reelle Stromfunktion Y7-,,, dieser Strömung, nämlich 


k=1 


Da wir aber aus dem im vorigen Abschnitt genannten Grunde für die Ermittlung der 
relativen Wirbelbahnen die reelle Stromfunktion Y7-,., der relativ-stationären Verdrängungs- 
strömung benötigen, müssen wir zu der Stromfunktion 77-,,, noch die Stromfunktion 7), der- 
jenigen Drehströmung addieren, welche durch den Übergang von einem raumfesten Koor- 
dinatensystem auf ein mit dem Kreiselrad rotierendes, also körperfestes Koordinatensystem 
zustande kommt. Es ist also 


Da das Kreiselrad entgegen dem Uhrzeigersinn mit der Winkelgeschwindigkeit © rotiert, 
erhalten wir für die Stromfunktion Y), den Ausdruck 


Also lautet die reelle Stromfunktion der relativ-stationären Verdrängungsströmung 
k=1 


8. Die Stromfunktion für die Durchflußströmung. Die Durchflußströmung rührt davon 
her, daß die Flüssigkeit der Pumpe in der Mitte durch ein Saugrohr mit anschließendem 
Leitapparat zugeführt werden muß. Bei dem von uns behandelten ebenen Problem ersetzen 
wir die Zuströmung durch eine im Drehpunkt gelegene Quelle. Die Wirkung des Leitapparates 
soll dadurch ersetzt werden, daß im Drehpunkt auch noch ein Wirbel angebracht wird, dessen 
Drehsinn mit dem Drehsinn des Laufrades übereinstimmt. Man bezeichnet eine solche Quelle, 
der ein Wirbel überlagert ist, als Wirbelquelle. Die in der Wirbelquelle entspringende 
Strömung würde bei ungestörtem Durchfluß nach logarithmischen Spiralen verlaufen. Wir 
wollen in diesem Abschnitt aber die Stromfunktion 7) derjenigen Strömung ermitteln, die 
in der Wirbelquelle entspringt und die ruhend zu denkenden Laufradschaufeln zirkulationslos 
umströmt. Diese Strömung pflegt man als Durchflußströmung zu bezeichnen. 

Die gesuchte Stromfunktion 7), für die Durchflußströmung kann aus der Stromfunktion 
Y, der nur von der Quelle herrührenden Strömung und der Stromfunktion Y', der nur vom 
Wirbel herrührenden Strömung gemäß der Gleichung 


(%6) 


zusammengesetzt werden. Für beide Strömungsanteile müssen die radialen eh 
zum System der Stromlinien gehören. 

Diese Forderung ist für den ersten Strömungsanteil ohnehin erfüllt. Denn die radialen 
Schaufeln vermögen die radiale Quellströmung nicht zu stören. Die reelle Stromfunktion Y%y 
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dieser Strömung kann daher sofort angeschrieben werden. Bezeichnen wir ihre Quellstärke 
mit Q, so lautet ihre reelle Stromfunktion D 


Um nun auch die reelle Stromfunktion für den zweiten Anteil der Durchflußströmung 
zu erhalten, gehen wir von der z-Ebene des Einheitskreises aus. Wir bringen dort in dem 
dem Punkt w,=0 entsprechenden Punkt z, (Gl. 4) einen linksdrehenden Wirbel von der 
Stärke «; an. Wegen der Transformation w = A (Gl. 2), die eine Singularität im Punkt w, 
zur Folge hat, verwandelt sich dieser Wirbel bei der Abbildung in die w-Ebene in einen 
Wirbel von n-facher Stärke. Nach dem Vorbild von Gl. (5) können wir das komplexe Poten- 
tial Fr (2) der durch diesen Wirbel erzeugten Strömung um die Schaufeln, ausgedrückt in 
den Koordinaten der z2-Ebene, sofort anschreiben. Es lautet 


(28). 


Um den zugehörigen Ausdruck für die reelle Stromfunktion zu erhalten, zerlegen wir diesen 
Ausdruck in Real- und Imaginärteil. Nach einigen Zwischenrechnungen ergibt sich unter 
Berücksichtigung von Gl. (4) die reelle Stromfunktion Y,r der vom Wirbel herrührenden 
Strömung zu 
(?) + (r- sin p)? 


r-csp 1)-+(r- sin 
q 


9. Die Stromfunktion für die Zirkulationsströmung um die Schaufeln. Damit trotz des 
idealen Flüssigkeitscharakters Kraftwirkungen zwischen Schaufel und Fördermittel auftreten 
können, muß nach dem Kutta-Joukowskyschen Satz neben den bisher behandelten Teil- 
strömungen auch noch das Bestehen von Zirkulationsströmungen um die Schaufeln angenommen 
werden. Der Drehsinn dieser Zirkulationen muß bei einer Zentrifugalpumpe mit dem Dreh- 
sinn des Laufrades übereinstimmen. 

Das komplexe Potential der Zirkulationsströmung um die Schaufeln wird in der z-Ebene 
einfach dargestellt durch das komplexe Potential eines im Mittelpunkt des Einheitskreises 
angebrachten Wirbels. Bezeichnen wir die Stärke dieses Wirbels mit u,, so hat die Zir- 
kulation um jede der n Schaufeln den Wert u4,-2r, da zum Mittelpunkt des Einheitskreises 
bei der Abbildung kein singulärer Punkt gehört. Der Drehsinn dieses Wirbels muß mit dem 
Drehsinn des Laufrades übereinstimmen, da derselbe bei der Abbildung erhalten bleibt. 

Das komplexe Potential für die Zirkulationsströmung um die Schaufeln kann also ein- 
fach in der z-Ebene als das komplexe Potential F, (z) eines linksdrehenden, im Mittelpunkt des 
Einheitskreises angebrachten Wirbels von der Stärke u, angesetzt werden. Es lautet 


Die zugehörige reelle Stromfunktion lautet 


10. Die Wirbelstärken .; und u, für tangentiales An- und Abströmen. Die Wirbelstärken 
u; und u, sollen für die Berechnung der Wirbelbahnen so bemessen werden, daß die bekannten 
Bedingungen des tangentialen Anströmens an die Eintrittskanten und des tangentialen Ab- 
strömens von den Austrittskanten der Schaufeln ohne Mitberücksichtigung der vom Wirbel- 
kranz herrührenden Strömung erfüllt sind. 

Wenn nämlich diese Bedingungen unter Mitberücksichtigung der vom Wirbelkranz her- 
rührenden Strömung erfüllt sein sollen, müssen die Wirbelstärken «; und «,, und damit auch 
die zugehörigen Zirkulationen von der jeweiligen Lage der Wirbel, und somit auch von der Zeit 
abhängig sein. Dies führt aber auf einen Widerspruch, auf den Hamel'”) aufmerksam gemacht 
hat. Wenn nämlich die Geschwindigkeitspotentiale ® dieser Zirkulationsströmungen, die nicht 
eindeutig sind, auch von der Zeit abhängen, treten in der Bernoullischen Gleichung die ent- 
sprechenden Glieder < auf. Dieselben sind aber vieldeutig, und damit auch der Druck p, 
was nicht zulässig ist. 


17,G. Hamel: a. a. 0.9. 


XUM 
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Hamel hat in der zitierten Arbeit ein anderes Beispiel behandelt, nämlich die Bewegung 
eines Wirbels um eine Buhne, bei dem er, ohne auf einen Widerspruch zu kommen, das Auf- 
treten unendlicher Geschwindigkeiten an der Buhnenkante durch die Überlagerung einer 
parabolischen, von der jeweiligen Lage des Wirbels abhängigen Strömung vermeiden 
konnte. Es läßt sich nun zeigen, daß im Falle der Überlagerung einer solchen, von der Lage 
des Wirbels abhängigen Strömung eine Routhsche Stromfunktion für die zugehörigen Wirbel- 
bahnen nicht mehr existiert. Diese Wirbelbahnen haben also einen anderen Charakter. Es 
bleibt deshalb in einem solchen Fall nichts anderes übrig, als nachı dem Hamelschen Vorbild 
mittels des Helmholtz-Kirchhoffschen Grundsatzes die Differentialgleichung der Wirbel- 
bahnen aufzustellen und zu integrieren. Formal ließe sich diese Rechnung auch bei unserer 
Aufgabe durchführen, wobei die Integration der Differentialgleichung, die eine recht komplizierte 
Form annimmt, allerdings auf graphisch-numerischem Wege erfolgen müßte. 

Aus diesen Gründen sollen die Wirbelstärken «; und «, so bestimmt werden, daß die 
Bedingungen des tangentialen An- und Abströmens wenigstens dann erfüllt sind, wenn die 
vom Wirbelkranz hervorgerufene Strömung nicht vorhanden ist. 

Zu dem Zweck gehen wir aus von dem komplexen Potential F(z) derjenigen Strömung, 
die sich aus der absoluten Verdrängungsströmung, aus der vom Wirbel im Drehpunkt Irer- 
rührenden Strömung und aus der Zirkulationsströmung um die Schaufeln zusammensetzt. Die 
radiale Quellströmung können wir hierbei auslassen, da dieselbe ohnehin bereits tangential 
an- und abströmt. Unter Berücksichtigung der Beziehungen (19), (28) und (30) haben wir also 


2. 


Damit zwecks Erfüllung der Bedingungen des tangentialen An- und Abströmens die absoluten 
Strömungsgeschwindigkeiten an den Schaufelenden endlich sind, muß der Differentialquotient 


ar für die den Schaufelendpunkten entsprechenden, auf dem Einheitskreis gelegenen Punkte 


z2= + 1 verschwinden. Diese beiden Punkte müssen also Staupunkte für die Bildströmung 
werden. Dadurch erhalten wir zwei Bestimmungsgleichungen für die Wirbelstärken «; und 
##s, die sich nach mehreren Umformungen auf die Formen 


dp, 


bringen lassen. 


11. Die Gleichung der Wirbelbahnen. Wir können jetzt die Gleichung der relativen 
Wirbelbahnen aufstellen, die von den Kranzwirbeln am Schaufelstern der w-Ebene beim Vor- 
handensein der Verdrängungsströmung, der Durchflußströmung und der Zirkulationsströmung 
um die Schaufeln durchlaufen werden. Zu dem Zweck brauchen wir nur die Routhsche 
Stromfunktion y„ (Gl. 14) und die gewöhnlichen Stromfunktionen Y7,. (Gl. 25), Po (Gl. 27), 
Py (Gl. 29) und 7, (Gl. 31) zueinander zu addieren und ihre Summe gleich einer Kon- 
stanten K zu setzen. Die Gleichung dieser Wirbelbahnen lautet also 


Dabei muß auch den Koordinaten der überlagerten Teilströmungen der Index 1 gegeben 
werden. Die in dieser Gleichung auftretenden Koeffizienten a; sind durch den Ausdruck (21), 
und die Wirbelstärken «; und 4, durch die beiden Ausdrücke (33) gegeben. Mittels der Be- 
ziehungen (17) und (18) kann diese Gleichung einheitlich in den Koordinaten der z-Ebene 
ausgedrückt werden. 

Da es nicht möglich ist, Gl. (34) nach r, oder , aufzulösen, bedienen wir uns einer 
graphisch-numerischen Methode, um zunächst zu einer Aufzeichnung der Wirbelbahuen in der 
z-Ebene zu gelangen. Hier stellen sie aber nicht etwa die Wirbelbahnen am Einheitskreis 
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dar, sondern lediglich diejenigen Kurven, die durch konforme Abbildung in die w-Ebene in 
die gesuchten Wirbelbahnen am Schaufelstern übergehen. 

Zur Übertragung dieser Kurven aus der z-Ebene in die w-Ebene bedienen wir uns 
zweier, durch konforme Abbildung auseinander hervorgehender Netze, die durch die ortho- 
gonalen Kurvenscharen o=konst. und = konst. gebildet werden. Während das in der 
w-Ebene gelegene Netz einfach aus Radien und konzentrischen Kreisen gebildet wird, er- 
halten wir die entsprechenden Kurvenscharen des Netzes in der z-Ebene durch Auswertung 
derjenigen Gleichungen, die sich aus den Gl. (17) und (18) durch Konstantsetzen von o und 
ergeben. 

In den Bildern 2 und 3 sind die beiden, einander entsprechenden Netze der w-Ebene 
und der z-Ebene für die im folgenden Abschnitt behandelten Beispiele dargestellt worden. 


/ 
\ / 
/ 
\f 
N \ / 
N 
/A| 
| 
/ / 
/ 
Bild 2. Bild 3. 


12. Einige Zahlenbeispiele.e. Wir wollen nun mit Hilfe der abgeleiteten Beziehungen 
einige Zahlenbeispiele durchrechnen und sodann die zugehörigen Wirbelbahnen aufzeichnen. 
Aus der Vielzahl der Möglichkeiten wollen wir hierfür vier besonders charakteristische Fälle 
herausgreifen, die sich, um die umfangreichen Zahlenrechnungen in erträglichen Grenzen zu 
halten, alle auf den in der Mitte geschlossenen radialen Schaufelstern beziehen sollen. Wir 
haben also in den abgeleiteten Beziehungen p=g zu setzen. Ferner sollen sich diese Fälle 
auf einen Schaufelstern mit a—=6 Schaufeln beziehen, deren Länge in der Weise festgelegt 
werden soll, daß p=64 gesetzt wird. Endlich soll in diesen Fällen der Schaufelstern mit der 
Winkelgeschwindigkeit rotieren. 


Hiermit ergaben sich der Innenradius o; und der Außenradius o„ nach den Formeln 
und =Yp+g zu —=0 und 0, =2,245. 

Ferner ergaben sich nach den Formeln (33) die Wirbelstärken «; und u, zu ;=0 und 
u; = 0,5989. 


Endlich ergaben sich nach Formel (21) die Koeffizienten a, der ersten sechs Glieder 
von Reihe (2) zu a4,=—08%9, ,=+0258, ,=— 0,125, a,=+008%0, a,= — 0,055, 
a,—= + 0,040. Für die Stromfunktion der absoluten Verdrängungsströmung erhalten wir unter 
Benutzung der ersten sechs Glieder dieser Reihe bereits einen für unsere Zwecke ausreichend 
genauen Ausdruck. Dabei ist die Abnahme der entsprechenden Glieder der Reihe im all- 
gemeinen noch stärker, als es in der Abnahme der Koeffizienten a, zum Ausdruck kommt, 


da dieselben noch mit multipliziert sind, 
ı 


Aus der in der Gleichung der Wirbelbahnen noch vorkommenden Wirbelstärke m der 
einzelnen Kranzwirbel, sowie aus der Quellstärke @ der radialen Quellströmung sollen nun 
vier Wertenpaare so zusammengestellt werden, daß dadurch die eingangs erwähnten vier Fälle 
unterschieden und an Hand von Kurvenbildern möglichst klar herausgearbeitet werden können. 


Zu dem Zweck treffen wir in den vier Fällen für m und @ die folgende Wahl: 


1 


| 
W028 


N 
>> 
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l.m=+1, 
2. m=—2, 
3. m=+1l, 9=+--; 


4.m=—-2, Q=+—. 


Um ein klares Bild der Wirbelbahnen in der w-Ebene zu erhalten, wählen wir die 
Differenzen AK zwischen den K-Werten je zweier aufeinander folgender Kurven in den vier 
Fällen wie folgt: 

1. AK=0,05; 3. AK=01; 
2. AK=05; 4. AK=06. 


13. Kurvendiskussion. In den Bildern 4, 5, 6 und 7 sind die in den vier Fällen ermittelten 
Wirbelbahnen wiedergegeben worden. Wie man sieht, kommen in diesen Bildern die für die 
einzelnen Fälle charakteristischen Eigenschaften der Wirbelbahnen deutlich zum Ausdruck. 

Zur Beurteilung der erhaltenen Kurvenscharen ist es wertvoll, in den einzelnen Fällen 
auch noch den Verlauf der relativen Stromlinien der zugehörigen Grundströmung zu kennen. 
Die Gleichung dieser Stromlinien erhalten wir ohne weiteres aus Gl. (34), wenn wir darin 
4w=0 setzen. Wir können hier davon Gebrauch machen, daß diese Stromlinienbilder in den 
dann noch verbleibenden beiden Fällen bereits von Kucharski‘*) ermittelt wurden. Diese 
Bilder, die sich auf die Schaufelzahl 20, wie überhaupt auf andere Zahlenwerte beziehen, sind 
in den Bildern 8 und 9 wiedergegeben worden. 

Bild 8 zeigt die bei der Rotation der radialen Schaufeln auftretende relativ-stationäre 
Verdrängungsströmung, welcher noch diejenige Zirkulationsströmung um die Schaufeln über- 
lagert wurde, die die Schaufelendpunkte in dem gezeichneten Relativstrombild zu Stau- 
punkten macht. 

Bild 9 zeigt diejenige Relativströmung, welche aus der in Bild 8 dargestellten Relativ- 
strömung durch Überlagerung einer radialen Quellströmung hervorgeht. 

Wenn nun in die durch die Bilder 8 und 9 dargestellten Strömungen einmal ein links- 
drehender, und einmal ein rechtsdrehender Wirbelkranz hineingesetzt wird, kommen wir 
wieder auf die von uns behandelten vier Fälle. Von den betreffenden Kranzwirbeln werden 
dann Wirbelbahnen von der Art der in den Bildern 4 bis 7 dargestellten durchlaufen. 

Der Vergleich von Bild 8 mit den ihm entsprechenden Bildern 4 und 5, sowie von Bild 9 
mit den ihm entsprechenden Bildern 6 und 7 zeigt einerseits den älınlichen Charakter der 
Kurvenscharen in größerer Entfernung vom Schaufelstern, sowie andererseits den verschieden- 
artigen Charakter in der Nähe des Schaufelsterns. In der näheren Umgebung des Drehpunktes 
zeigen die Wirbelbahnen eine auffallende Ähnlichkeit mit denjenigen bislang noch nicht ver- 
öffentlichten Wirbelbahnen, die Kucharski für den rotierenden radialen Schaufelstern mit 
vom Drehpunkt bis ins Unendliche reichenden Schaufeln für die gleichen Fälle und für die 
gleiche Schaufelzahl, aber für sonst andere Zahlenwerte gefunden hat. 

Die genannten Ähnlichkeiten sind leicht zu verstehen. In einiger Entfernung vom 
Schaufelstern ist sein Einfluß auf das Aussehen der Wirbelbahnen bereits gering, und in 
größerer Entfernung findet so gut wie keine Beeinflussung mehr statt. Da die gegenseitige 
Beeinflussung der alsdann auf einen entsprechend großen Kreis verteilten Kranzwirbel dann 
auch nicht mehr ins Gewicht fällt, müssen die Wirbelbahnen dort praktisch mit den Strom- 
linien der überlagerten Relativströmung zusammenfallen. Die Wirbel werden dort einfach 
von der Strömung mitgenommen. Dagegen wird das Aussehen der Wirbelbahnen in der 
näheren Umgebung des Drehpunktes in erster Linie von den näher gelegenen Konturen be- 
einflußt. Es wird also für dieses Aussehen nicht viel ausmachen, ob die Schaufeln eine 
endliche Länge besitzen, oder ob sie bis ins Unendliche reichen. 

Die Bilder 4 bis 7 lassen ferner erkennen, daß der Einfluß der Schaufeln auf das Aus- 
sehen der Wirbelbahnen in ihrer nächsten Umgebung den Einfluß der überlagerten Teil- 
strömungen bei weitem überwiegt. Die große Dichte der Kurvenschar in nächster Nähe der 
Schaufeln ist kennzeichnend für die dort herrschenden großen Wirbelgeschwindigkeiten, die 
mit zunehmender Annäherung an die Schaufeln bekanntlich sogar dem Grenzwert „unendlich“ 
zustreben. 

Aus früheren Arbeiten ‘®) ist ferner bekannt, daß ein in genügender Nähe einer Schaufel 
befindlicher Wirbel sich längs dieser Schaufel gerade in entgegengesetzter Richtung vorwärts 
bewegt, als ein mit dem Drehsinn des Wirbels an der Schaufel abrollendes Rad sich fort- 

18) W. Kucharski: Strömungen einer reibungsfreien Flüssigkeit bei Rotation fester Körper. München und 
Berlin 1918, Bild 55 und 57. 


19) Z.B. W. Kucharski:a. a. O.?). 
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bewegen würde, und daß die in der Nähe des Drehpunktes eines innen geschlossenen radialen 
Schaufelsterns gelegenen Wirbelbahnen hyperbelartig innerhalb eines jeden Sektors verlaufen. 
Auch ist aus der Kucharskischen Arbeit bekannt, daß diese Wirbelbahnen bei fehlender 
Quellströmung symmetrisch zum Mittelradius des Sektors verlaufen, wogegen dieselben durch 
Überlagerung einer Quellströmung je nach dem Drehsinn der Wirbel nach der einen oder 
anderen Seite verschoben werden. 

Die Bilder 4 und 5 zeigen, daß bei fehlender Quellströmung die Wirbelbahnen in ihrer 
Gesamtheit symmetrisch zum Mittelradius des Sektors verlaufen. Dagegen zeigen die Bilder 6 
und 7, daß bei vorhandener Quellströmung keine Symmetrie mehr vorhanden ist. 

Es ist bemerkenswert, daß nach Überlagerung der radialen Quellströmung die Wirbel- 
bahnen doch nicht im Ursprung der Quelle ihren Anfang zu nehmen scheinen. Wenn man 
aber die nicht geschlossenen Wirbelbahnen, auf denen sich die Wirbel nach außen hin end- 
gültig vom Laufrad entfernen, in der entgegengesetzten Richtung zurückverfolgt, bemerkt 
man, daß dieselben sich immer näher und näher an den Schaufelstern heranschlängeln, so 
daß sie auf den Zeichnungen wegen der dort herrschenden großen Dichte der Wirbelbahnen 
sehr bald nicht weiter verfolgt werden können. Aus den Zeichnungen geht also nicht hervor, 
an welcher Stelle der Schaufelkontur sich die Wirbel ablösen. Die Beantwortung dieser Frage 
auf Grund unserer Zeichnungen würde aber ohnehin keinen Sinn haben, da sich dieselben 
auf eine reibungsfreie Flüssigkeit beziehen, in der die Entsteliung von Wirbeln unmöglich ist. 

Eine besondere Eigentümlichkeit der in den Bildern 4 bis 7 dargestellten Wirbelbahnen 
sind die zahlreichen relativen Ruhelagen der Wirbel. (Die Kreuzungspunkte zweier Wirbel- 
bahnen sowie die isolierten „Mittelpunkte“ der geschlossen umlaufenden Wirbelbahnen.) Die 
Beantwortung der Frage nach der Stabilität dieser Ruhelagen bereitet keine Schwierigkeiten 
mehr, nachdem die Wirbelbahnen in den betrachteten Fällen gezeichnet vorliegen. Einer 
Arbeit von Kneschke”), die sich mit den Ruhelagen von Wirbeln beschäftigt, können wir 
nämlich folgendes entnehmen: 

Es gibt drei Arten von Ruhepunkten, nämlich Kreuzungspunkte, Rückkehr- 
punkte und isolierte Punkte. 

Bezüglich dieser Ruhepunkte gelten die folgenden Sätze: 

1. Kreuzungspunkte sind bedingt stabile Ruhepunkte. 
2. Rückkehrpunkte sind bedingt stabile Ruhepunkte. 
3. Isolierte Punkte sind stabile Ruhepunkte. 

Da die in den Bildern 4 bis 7 vorhandenen Ruhepunkte Kreuzungspunkte und isolierte 
Punkte sind, haben wir es hier mit bedingt stabilen und stabilen Ruhelagen zu tun. Die 
stabilen Ruhelagen rühren zum Teil von den entsprechenden, aus den Bildern 8 und 9 zu 
ersehenden stabilen Ruhelagen der Grundströmung her. 

Daß für den Wirbelkranz stabile Rulıelagen möglich sind, ist — besonders für die dem 
Betrieb als Pumpe entsprechenden Fälle (Bild 6 und 7) — ein wichtiges Resultat unserer 
Untersuchungen, denn die Wirbel können in solchen Lagen nicht von dem Seclıaufelrad los- 
kommen. Da diese Ruhelagen relative Ruhelagen sind, rotiert der eine solche Ruhelage ein- 
nehmende Wirbelkranz mit dem Schaufelrad so, als ob er mit ihm fest verbunden wäre. Die 
sämtlichen Wirbel des Kranzes durchlaufen hierbei als absolute Bahn den gleichen Kreis. 

Bei fehlender Quellströmung haben wir im Innern der Schaufelräume symmetrisch zu ihrer 
Mittellinie gelegene Bereiche geschlossener, relativer Wirbelbahnen um die zugehörigen stabilen 
relativen Ruhelagen (Bild 4 und 5). Diese Bereiche füllen nahezu die ganzen Schaufelräume aus. 

Nach Überlagerung der Quellströmung werden diese Bereiche geschlossener Wirbel- 
bahnen in den einzelnen Schaufelräumen an die Seite gedrängt und dabei verquetscht und ver- 
kleinert, wodurch an der anderen Seite eines jeden Schaufelraumes Platz für solche Wirbel- 
bahnen geschaffen wird, auf denen die sich von den Schaufeln ablösenden Wirbel das Schaufel- 
rad verlassen können (Bild 6 und 7). 


Zusammenfassung. Es wurde die ebene Bewegung eines Wirbelkranzes an einem als 
Kreiselpumpe betriebenen radialen Schaufelstern von endlichen Schaufellängen unter der Vor- 
aussetzung untersucht, daß die Bedingungen des tangentialen An- und Abströmens ohne Mit- 
berücksichtigung der vom Wirbelkranz herrührenden Strömung erfüllt sind. Mittels des 
Routhschen Verfahrens wurde die Gleichung der von den einzelnen Kranzwirbeln durch- 
laufenen relativen Wirbelbahnen aufgestellt. Diese Gleichung wurde dann für vier besonders 
charakteristische, sich auf den innen geschlossenen radialen Schaufelstern beziehende Fälle 
numerisch ausgewertet, und die zugehörigen relativen Wirbelbahnen wurden graphisch dar- 
gestellt. Die Ergebnisse, insbesondere auch die mit den Ruhepunkten zusammenhängenden 
Fragen, wurden sodann besprochen. 140 


20) A. Kneschke: Über die Ruhelagen von Wirbeln. Z. angew. Math. Mech. Bd. 14 (1934), 8. 178. 
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Untersuchungen zur Verminderung störender Rüttel- 
schwingungen an Kreiselgeräten. 
Von K. Magnus in Göttingen. 
Bericht des Institutes für angewandte Mechanik der Universität Göttingen!). 


1. Einleitung. Jeder Kreisel besitzt restliche Unwuchten, die beim Lauf Zentrifugal- 
momente und damit störende Schwingungen hervorrufen. Es soll untersucht werden, wie die 
Trägheitsmomente eines Kreiselgerätes zu wählen sind, damit die durch die Zentrifugalmomente 
erregten Schwingungen möglichst klein bleiben. 


Die Untersuchungen wurden angeregt durch das Verhalten eines Kreiselhorizontes, der 
in der Werkstatt des Institutes für angewandte Mechanik der Göttinger Universität hergestellt 
worden war. Der Horizont bestand aus einem elektrisch angetriebenen Kreisel mit senk- 
rechter Drehachse, der durch Kardanringe Bewegungsfreiheit um die waagerechten Achsen 
bekam. Das System hatte eine geringe Schwerpunktstieferlage. 

Obwohl der Kreisel bei fest eingespannter Kappe ruhig lief und nur die üblichen kritischen 
Drehzahlen zeigte, war der Lauf bei frei beweglich aufgehängter Kappe so unruhig, daß der 
Kreisel nicht auf volle Drehzahl gebracht werden konnte. Die überaus starken Rüttelschwin- 
gungen ließen eine Überbeanspruchung des Materials und damit eine Zerstörung des Apparates 
befürchten. 

Man hat es hier offenbar mit erzwungenen Schwingungen des Kreiselpendels zu tun, 
deren erregende Kraft die durch restliche Wuchtfehler erzeugte umlaufende Zentrifugalkraft 
ist. Nun sind zwar die durch Unwuchten erzeugten Schwingungen sowie die damit zusammen- 
hängende Selbsteinstellung elastischer Wellen in der Literatur sehr ausführlich behandelt 
worden’). Die meisten Verfasser rechnen jedoch, entsprechend ihrer etwas anderen Problem- 
stellung, mit festen Lagern. Lediglich Den Hartog berücksichtigt die Nachgiebigkeit der 
Lager, nimmt jedoch deren Masse so groß an, daß die Kopplung mit den Lagermassen ver- 
nachlässigt werden kann. 

Im folgenden soll nun das Problem für den speziellen Fall des Kreiselpendels behandelt 
werden. Der Kreisel soll mit elastischer Welle in einer Kappe liegen. Kreisel und Kappe 
erhalten durch Lagerung in einem Kardanring Bewegungsfreiheit um zwei zueinander senk- 
rechte horizontale Achsen. Wenn der Kreisel läuft, so entstehen durch die vorhandenen Un- 
wuchten Zentrifugalmomente. In der folgenden Rechnung wird nun ein durch die elastische 
Kreiselwelle gekoppeltes System von Kreiselkörper und den Massen des Kardangehänges be- 
trachtet, das durch die Zentrifugalmomente zu erzwungenen Schwingungen erregt wird. 


Die Untersuchungen werden zeigen, daß das Trägheitsmoment der Kreiselkappe eine 
bestimmte Größe haben muß, damit Resonanzerscheinungen mit der Eigenfrequenz der Nutation 
des Kreisels vermieden werden. Es bleibt jedoch stets die Resonanz mit der durch die 
elastische Welle verursachten Querschwingung des Kreiselkörpers bestehen. 


2. Die erregende Zentrifugalkraft. Die die Schwingungen verursachenden Unwuchten 
können dynamischer oder statischer Natur sein. Statische Unwuchten, die durch einen kleinen 
Fehlerabstand des Schwerpunktes von der Drehachse hervorgerufen werden, sollen nicht be- 
rücksichtigt werden. Sie sind zudem bei technischen Kreiseln durch Auswuchten fast restlos 
beseitigt. Die durch statische Unwuchten bedingten Schwingungen sind in der Literatur aus- 
führlich behandelt worden; sie sind im wesentlichen Translationsschwingungen, die sich den 
durch die dynamischen Unwuchten hervorgerufenen Drehschwingungen überlagern und so das 
durch die letzteren erzeugte Bild der Schwingungserscheinungen etwas verändern. 


Die dynamischen Unwuchten sind durch einen Fehlerwinkel zwischen der Drehachse 
und der Hauptträgheitsachse gekennzeichnet. Die Größe der dadurch hervorgerufenen Zentri- 
fugalmomente soll für einen symmetrischen Kreiselkörper bestimmt werden. Es sei in einem 
kreiselfesten Koordinatensystem A, B, C die Achse C die ausgezeichnete Hauptträgheitsachse. 
Sie bilde mit der Drehachse den Fehlerwinkel e (Bild 1). Die Achse B sei senkrecht auf 
der durch Drehachse und Hauptträgheitsachse bestimmten Ebene. Dann können wir das 
Zentrifugalmoment aus den Eulerschen Gleichungen der Kreisellehre bestimmen. Diese 
lauten im allgemeinen Fall: 

1) Diese Untersuchungen wurden auf Anregung von Herrn Prof. Schuler durchgeführt, dem auch an dieser 
Stelle für viele wertvolle Hinweise gedankt sei. 

2) Z. B. Pöschl: Z. angew. Math. Mech. Bd. 3 (193), S. 297. Karas: Ing.-Arch. Bd. 1 (1930), S. 84 und 158. 


Ferner die Bücher: Stodola: Dampfturbinen, 3. Aufl. Berlin 195. Klein-Sommerfeld: Theorie des Kreisels, 
Heft IV, Leipzig 1910. Den Hartog: Mechanische Schwingungen, deutsch von G. Mesmer, Berlin 1936. 


166 Magnus, Untersuchung störender Rüttelschwingungen an Kreiselgeräten 


9,84 = — %) we — Zu, 
= (I — we — ZB, 
oc=(9ı 9;) Ze: 
Dabei sind ©,, 9, ©c die Hauptträgheitsmomente, 
Ze die Zentrifugalmomente und ®,, ®c die 
Komponenten der Drehgeschwindigkeit um die kreiselfesten 
Achsen A, B, ©. Im vorliegenden Fall ist nun: 
94 Op, 
=wsine, 


Damit ergibt sich: \ 
Zı 


Bild I. Zerlegung der Dreh- 
= (Ic — 9,) sin cose. geschwindigkeit am Kreisel. 

Das Zentrifugalmoment wächst also mit dem Quadrat der Umlauffrequenz, es hat ein 
Maximum bei einem Fehlerwinkel von 45°. Da in unserem Fall jedoch nur kleine & vor- 
kommen, so können wir für das Moment einfach schreiben: 

Ze — 9,) ©? . (1). 

Da ferner das Zentrifugalmoment die Differenz der Trägheitsmomente we Faktor hat, so 
könnte man annehmen, daß für einen Kugelkreisel, d. h. für einen Kreisel mit kugelförmigem 
Trägheitsellipsoid, das Zentrifugalmoment verschwinden müßte. Dem widersprechen jedoch 
die Ergebnisse der Erfahrungen. Die weitere Rechnung wird nun auch die Erfahrungstatsachen 
theoretisch bestätigen. Es soll jedoch hier noch eine ganz allgemeine Bemerkung über den 
Fehlerwinkel e bzw. über das Zentrifugalmoment gemacht werden. 

Bei einem Kugelkreisel ist jede Achse gleichberechtigt, jede Drehachse kann als Haupt- 
trägheitsachse aufgefaßt werden. Die geringste Inhomogenität des Kreiselmaterials oder ganz 
unbedeutende Massenverlagerungen durch Verbeulen können nun bewirken, daß eine Achse 
als ausgezeichnete Hauptträgheitsachse erscheint, die mit der Drehachse einen beliebig großen 
Winkel einschließt. Andererseits wird dieselbe Inhomogenität bei einem stark abgeplatteten 
Kreisel nur eine geringfügige Verlagerung der Hauptträgheitsachse zur Folge haben. Da es 
also offensichtlich auf das Produkt (9. — 9 ,)e ankommt, so soll dieses für eine vorgegebene 
Inhomogenität untersucht werden. 

Das Trägheitsmoment des jetzt als fehlerfrei angenommenen symmetrischen Kreisel- 
körpers in einer beliebigen Richtung «a ist gegeben durch: 

O.sin’a 

Die durch Materialfehler und schlechte Bearbeitung bewirkten Abweichungen von der 
ideaien Massenverteilung im Kreiselkörper denken wir uns zu einem besonderen Körper zu- 
sammengefaßt, dessen Trägheitsellipsoid wir als symmetrisch annehmen wollen. Er habe die 
Hauptträgheitsmomente ©, und ©,, deren Achsen gegenüber dem als fehlerlos angenommenen 
Kreiselkörper um den Winkel @ verdreht sein sollen. Dann stellt sich das Fehlerträgheits- 
moment in der Richtung a dar durch: 

Or = 9, 608? (a—y)+ 9, sin’ (a—p). 
Um das Trägheitsellipsoid des wirklichen Kreiselkörpers zu bekommen, müssen wir $% 
und ©, addieren. 
(2) 
mit 
=9c+ 9, sin’p+9,cos’p, 


D= (8, — 9,)sin’ 


©, und ©- sind dabei die Trägheitsmomente des fehlerhaften Kreiselkörpers, D ist eine 
nur von der Inhomogenität abhängige Konstante. Um nun die Hauptträgheitsachsen des 
wirklichen Kreiselkörpers zu finden, differenzieren wir (2) und bekommen: 


2a+2Dcos2a. 


\ N 

| 
| xum 
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Für eine Hauptträgheitsachse verschwindet dieser Ausdruck. «a gibt in diesem Fall die 
Abweichungen der Hauptträgheitsachsen von idealem und wirklichen Kreiselkörper, also den 
Fehlerwinkel e an. Man bekommt: 


2D 
(94 -%e) 


Da wir e als klein vorausgesetzt hatten, so können wir die Tangente durch den Winkel 
ersetzen. Multipliziert man dann e noch mit dem im Zentrifugalmoment Gl. (1) vorkommenden 
Ausdruck (9c — © ,), so ergibt sich: 


Z 1 
Dez 


Das bedeutet also, daß das Zentrifugalmoment für eine vorgegebene Inhomogenität eine 
von der Form des Kreisels unabhängige Größe ist. 


Die GI. (3) zeigt ferner, weshalb für einen Kugelkreisel die Verhältnisse am ungünstigsten 
liegen. Beim Kugelkreisel ist es seiner Form wegen wahrscheinlicher als beim abgeplatteten 
oder beim verlängerten Kreisel, daß der Winkel x den Wert 45° erreicht, daß also der in 
Gl. (3) auf der rechten Seite stehende Ausdruck ein Maximum wird. Für den abgeplatteten 
Kreisel ist 9>45°, für den verlängerten Kreisel ist 9 <45° am wahrscheinlichsten. 


Das für das körperfeste System abgeleitete Zentrifugalmoment (1) läuft vom raumfesten 
System aus gesehen mit der Drehgeschwindigkeit » um. Wir wollen es in seine Komponenten 
in Richtung der raumfesten x- und y-Achsen, die der folgenden Rechnung zugrunde gelegt 
werden, zerlegen. Dabei können wir statt des Momentenpfeiles selbst dessen Projektion in 
der waagerechten Ebene nehmen, da nur kleine Winkelausschläge des Kreisels betrachtet 
werden sollen. Damit bekommen wir die Momentkomponenten: 


= Zpsinot, 


Z, = — Zp wos ot. 


3. Die Differentialgleichungen der erzwungenen Schwingungen des Kreiselpendels. Wie 
schon erwähnt wurde, sollen nur kleine Auslenkungen des Kreiselpendels betrachtet werden. 
Wir können dann in einem raumfesten Koordinatensystem (Bild 2) 
die Bewegungen des Kreiselpendels durch die Föpplschen 
Gleichungen der Kreisellehre beschreiben®). Als äußere Momente 
haben wir darin zu setzen: 


1. Das Schweremoment des Kreiselkörpers: 5: a; =Ge- 
wicht des Kreiselkörpers, s; — Abstand Schwerpunkt-Auf- 
hängepunkt, Auslenkung aus der Gleichgewichtslage). 

2. Das durch die elastische Welle übertragene Kopplungs- 
moment zwischen Kreiselkörper und Kardangehänge. 

3. Das Zentrifugalmoment. 


Bild 2. Koordinatensystem. 


Zu den beiden Föppl-Gleichungen treten dann noch die 
Schwingungsgleichungen des Gehänges (Kreiselkappe und Kardan- 
ring), in denen als äußere Momente die Schweremomente @,s,a bzw. G,s,ß des Gehänges 
(ohne Kreisel!) in den beiden Ebenen sowie die Kopplungsmomente auftreten. Damit ergibt 
sich das Gleichungssystem: 


Il 9, + ce la —a)=(I9c — 9 ,)Ew*sin wt, 
d’ d 
0.58 - 7077 + (Pr — (Ic — 9) cost, 


(4). 
III. G,,a+cla —a,)=0, 


3) Siehe Schuler: Die Föpplsche Methode zur Lösung von Kreiselproblemen. Z. angew. Math. Mech. Bd. 15 
(1935), S. 88 bis 91. 
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Dabei bedeuten: 
©c = Trägheitsmoment des Kreiselkörpers um seine Drehachse. 
0 ‚= Trägheitsmoment des Kreiselkörpers um eine Querachse. 
9,n = Trägheitsmoment des Gehänges um die x-Achse. 
©,n = Trägheitsmoment des Gehänges um die y-Achse. 
Gk Sk, die obengenannten Schweremomente. 
c = Kopplungskoeffizient der elastischen Kreiselwelle. 
a — Auslenkung des Gehänges um die x-Achse. 
ß = Auslenkung des Gehänges um die y-Achse. 
a; — Auslenkung der Hauptträgheitsachse des Kreiselkörpers um die x-Achse. 
Pr = Auslenkung der Hauptträgheitsachse des Kreiselkörpers um die %y-Achse. 


Im Differentialgleichungssystem (4) tritt uns ein Fall entgegen, in dem das in den Föpp|- 
Gleichungen auftretende Moment nicht nur vom Winkel allein abhängig ist, sondern vom 
Winkel und von der Zeit t. Wir haben ein inhomogenes System vor uns, bei dem I.mit II 
in den Geschwindigkeitsgliedern (durch die Kreiselkräfte), ferner I mit III und II mit IV in 
den Kraftgliedern (durch die elastische Welle) gekoppelt sind. 

Bevor das Gleichungssystem (4) in voller Allgemeinheit behandelt wird, soll noch kurz 
auf den Fall der starren Kreiselwelle eingegangen werden, 


4. Betrachtung der Resonanzfunktionen für eine starre Welle. Im Grenzfall c=» wird 
a«=a und fx=Pß, d. h. der Kreiselkörper macht dieselben Bewegungen wie die Kappe. Wir 
können dann im Gleichungssystem (4) die Gl. I mit III und II mit IV zusammenfassen und 
bekommen mit Einführung der Abkürzungen: 


9. =94 +9 = Trägheitsmoment des Gesamtsystems (d. h. Kreisel und Ge- 
hänge!) um die x-Achse, 
= Trägheitsmoment des Gesamtsystems um die y-Achse, 
+ @, s, = Gesamtes Stabilitätsmoment um die x-Achse, 
G,8,= + @, s, = Gesamtes Stabilitätsmoment um die y-Achse, 
zwei Gleichungen, die mit Hilfe des Ansatzes 
a=4Asinot, 
ß=Beoswt 
für A und B die Lösungen ergeben: 


(94 — 90) 8, — (9, + 96) 
(6, 8, — ©?) (G, 8, — — w* 
(%— 5, — (9, + 80) 
(@, 8, — 9% 0°) (@, 8, — 9,0?) — 90? 


A= 


8). 
B= 


Der Nenner ist den beiden Resonanzfunktionen gemeinsam. Er ist eine quadratische 
Form in @*, deren Nullstellen sich sofort angeben lassen. Es läßt sich zeigen, daß stets zwei 
reelle Lösungen für &® vorhanden sind, da die Diskriminante der quadratischen Form stets 
positiv ist. Das Vorzeichen dieser Lösungen wird durch den Quotienten: 


Q G, G, 8; 
bestimmt. Da aus Stabilitätsgründen @, s, und G, s, entweder beide positiv oder beide negativ 


sind, so ist das Produkt beider stets positiv. Das Vorzeichen von @ wird also nur durch den 
Nenner bestimmt. Wir bekommen: Für 


haben wir stets zwei reelle Lösungen für ®, also zwei Unendlichkeitsstellen der Resonanz- 
funktion. Ist «>1, so gibt es nur eine reelle Nullstelle für . Im Grenzfall a=1 rückt 
die zweite Nullstelle zuo=x». 


Die Nullstellen der Resonanzfunktion lassen sich aus (5) sofort entnehmen. Sie liegen 
bei =0 sowie bei: 
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Der Verlauf der Resonanzfunktionen ist in | per 
Bild 3 für drei typische Werte von u aufgetragen. 
Den Kurven sind die an dem vorgenannten Ho- A 


rizont gemessenen Werte zugrunde gelegt. Die 
Resonanzstellen treten im vorliegenden Fall schon 
bei sehr niedrigen »-Werten auf, so daß sie f 

praktisch vollkommen ungefährlich sind, denn 7 Bi 
bei einer Umdrehung des Kreisels in 2x Se- 
kunden sind die auftretenden Massenkräfte so 
gering, daß eine Beschädigung des Apparates w-7 


ausgeschlossen ist. Jedoch können die Resonanz- 
stellen bei entsprechend stärkerer Fesselung des 4 
Kreisels auch wesentlich höher rutschen, wie 
man aus (5) sofort entnehmen kann. 

Von Bedeutung ist noch der Grenzwert, 
dem die Resonanzkurven für wachsendes & zu- 
streben. Man bekommt z. B. für As: 


(9c + 9) 

Der Grenzwert ist um so kleiner, je weiter | 

wir vom Kugelkreisel entfernt sind, d. h. je mehr A 
das Trägheitsellipsoid des Systemes Kreiselkörper 

+ Kappe von der Kugelgestalt abweicht. Im 
Sonderfall = 9, 9, bzw. u=1 hat die Re- 
sonanzfunktion für große » Parabeleigenschaften, 7 2 ® 

sie geht also mit >» gegen Unendlich. Das Bild 3. Resonanzfunktion bei starrer Kreiselwelle 
ist klar, wenn man bedenkt, daß dann die Um- u nn En 
lauffrequenz stets gleich der Nutationsfrequenz ist. Dieser Fall muß also vermieden werden. 
Es zeigt sich jedoch, daß das Verhalten der Resonanzfunktion für große &-Werte durch 
die Elastizität der Kreiselwelle ganz erheblich beeinflußt wird. 


5. Die Resonanzfunktion für eine elastische Welle. Wir kehren nun zum Gleichungs- 
system (4) zurück und führen dort dimensionslose Größen ein, um der Rechnung größere 
Allgemeinheit zu verleihen. Wir nennen: 


S 


ce 


Dabei ist » die Eigenfrequenz der Drehschwingung des Kreiselkörpers um eine Quer- 
achse bei fest eingespannten Lagern. y ist als dimensionslose Frequenzgröße ein Maß für 
die Umlauffrequenz des Kreisels. Das Trägheitsmoment der Hängung um die x-Achse setzt 
sich zusammen aus dem Trägheitsmoment der Kreiselkappe und dem Trägheitsmoment des 
äußeren Kardanringes. Nimmt man, was meist der Fall ist, die Kappe als symmetrisch an, 
so kann man u, zerlegen in «, + 4,, wobei u, das Verhältnis des Trägheitsmomentes des 
Kardanringes zu ©, bedeutet. 

Mit diesen Beziehungen und dem Ansatz: 


a4=4Asinot, a=(Csinwt, 
Pr = Beoswt, =Deswt 
bekommt man für die gesuchten Amplitudenfunktionen A, B, C, D das Gleichungssystem: 
D=—(1—6)ey? 
—A+Cl(i+a,— =0 
—B+Di+a,— u, y?) 
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Da es hier nur auf das Verhalten der Resonanzfunktionen für große » ankommt, d. h. 
auf Drehzahlen, die in der Größenordnung der Betriebsdrehzahlen liegen, so können wir im 
System (7) unbedenklich die Größen a,, a,, a, vernachlässigen, da sie um etwa zwei Zehner- 
potenzen kleiner als Eins sind. Das bedeutet praktisch eine Vernachlässigung der Schwere- 
momente und der durch diese hervorgerufenen langsamen Präzessionen. 


Die Lösung des Gleichungssystemes (7) ist: 


4 für die Kappe, 

N für den 
Kreisel 


Dabei ist der Nenner N allen Resonanzfunktionen gemeinsam: 

Der Verlauf der vier Resonanzfunktionen in Abhängigkeit von der dimensionslosen Um- 
lauffrequenz y ist in Bild 4 aufgetragen, der Rechnung sind die an dem ausgeführten Horizont 
gemessenen Werte für die Trägheitsmomente zugrunde gelegt. Die ausgezogenen Kurven stellen 
die Resonanzkurven € und D, also die 
Auslenkungen der Kreiselkappe dar, die 
gestrichelten Kurven entsprechen den 
Auslenkungen der Hauptträgheitsachse 
des Kreisels in den beiden Ebenen. Die 
Betriebsdrehzahl lag bei dem unter- 
suchten Gerät etwa bei y=3, also ober- 
halb der Resonanzstellen. Der Kreisel 
muß also beim Anlaufen diese Stellen 
durchlaufen. Dabei ist eine Resonanz- 
stelle um so gefährlicher, je breiter sie 
ist, denn eine um so längere Zeit braucht 
der Kreisel, um sie zu durchlaufen. Es H 
kann sogar der Fall eintreten, daß die ; 
Drehgeschwindigkeit nicht über den 
Wert einer solchen kritischen Drehzahl 
gesteigert werden kann, nämlich dann, 
wenn alle von außen durch den Antrieb in das System hineingesteckte Energie durch das 
Anfachen der Rüttelschwingungen verbraucht wird, so daß dann keine Energie zum Steigern 
der Drehzahl mehr übrigbleibt. Ein solcher Fall hat bei dem beschriebenen Horizontmodell 
tatsächlich vorgelegen, denn es gelang zwar ohne Schwierigkeit, den Kreisel bei festgehaltener 
Kappe auf Betriebsdrehzahl zu bringen, doch bei frei hängendem System war ein Steigern 
der Drehzahl nicht möglich, wollte man nicht das Gerät zerstören. 

An den Resonanzkurven des Bildes 4 interessieren uns 

1. die Nullstellen, 

2. die Unendlichkeitsstellen, 

3. die Grenzwerte, denen die Resonanzfunktionen für große » zustreben, denn diese 
entsprechen etwa den Werten der Resonanzfunktionen für die Betriebsdrehzahlen. 


Die Nullstellen lassen sich aus den Gl. (8) sofort angeben. Wir erhalten: 


Bild 4. Resonanzfunktionen für Kreiselkörper und Kreiselkappe. 


1 1 
für C: y,? 
für die Kappe, 
(9). 
für A: und 
für den Kreisel 
für B: und 9», 


XUM 
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Die Unendlichkeitsstellen ergeben sich aus den Nullstellen des Nenners von (8): 
1 1 1 1 ı2 
— +, - t >. 10). 


Da sich der Radikand als Summe zweier Quadrate darstellt, so kann eine komplexe 
Lösung für y»?” nicht vorkommen; dagegen können imaginäre Lösungen für y» auftreten, 
wenn y»? negativ wird. 

Es sei nun die Frage vorgelegt: Wie kann man die Lage der Nullstellen und Unend- . 
lichkeitsstellen so beeinflussen, daß sich die Eigenschaften des Kreiselverbandes bessern ? 
Man könnte zunächst eine Nullstelle der Resonanzfunktion in die Betriebsdrehzahl legen. 
Das ist jedoch aus zwei Gründen abwegig. Einmal tritt eine solche Nullstelle niemals bei 
allen vier Resonanzfunktionen gleichzeitig auf. Das ersieht man aus (9), da ja wegen des 
endlichen Trägheitsmomentes des Kardanringes stets u,=+ u, ist. Man könnte also nur eine 
Beruhigung jeweils in einer Ebene erzielen. Zweitens aber würde man durch ein solches 
Höherlegen der Nullstellen die Lage der Unendlichkeitsstellen ungünstig beeinflussen. 

Man wird versuchen, die Unendlichkeitsstellen zu kleineren Drehzahlen hin zu ver- 
schieben, um sie auf diese Weise ungefährlicher zu machen. Das kann durch Verändern von 
Trägheitsmomenten geschehen. Man kann entweder am Kardanring Zusatzgewichte anbringen 
und damit u, vergrößern, oder man kann das Trägheitsmoment der Kreiselkappe durch An- 
bringen von Zusatzgewichten vergrößern; damit wird «, und wegen 4,=u,+ u, auch u, ge- 
ändert. Den Kreisel selbst wollen wir zunächst als gegeben betrachten, so daß eine Änderung 
seiner Trägheitsmomente vorläufig ausgeschlossen sein soll. 

Bild 5 zeigt die Verschiebung der Nullstellen und der Unendlichkeitsstellen mit einer 
Vergrößerung von u,. (N, bedeutet Nullstelle von A usw. U bedeutet Unendlichkeitsstelle.) 
Dabei ist zu beachten, daß gemäß Gl. (9) die Nullstellen von € bzw. D gleichzeitig Nullstellen 
von A bzw. B sind. Bild 5 zeigt, daß die Vergrößerung von u, wenig Einfluß auf die Ver- 
schiebung der Nullstellen sowie der oberen Unendlichkeitsstelle hat. Dagegen verschwindet 
die untere Unendlichkeitsstelle für u, größer als etwa 1,2 vollkommen, d.h. der Nenner von 
Gl. (8) hat für größere u, nur noch eine reelle Nullstelle. 


\ 
2 Bild 5. Verschiebung der Nullstellen und 
M Unendlichkeitsstellen der Resonanzfunktion 
mit einer Vergrößerung des Kardanringes. 
25 
4, Bild 6. Verschiebung der Nullstellen und 
Unendlichkeitsstellen der Resonanzfunktion 
mit einer Vergrößerung der Kreiselkappe. 
\ Bild 7. Resonanzfunktionen für Kreisel- 
15 körper und Kreiselkappe für größeres Träg- 
v heitsmoment der Kappe. 
05 20 75 20 25 
Bild 5. 


IN 
A 
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Bild 7. 
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Bringt man die Zusatzgewichte nicht am Kardanring, sondern an der Kreiselkappe selbst 
an, so ändern sich die Verhältnisse nach Bild 6. Man sieht durch Vergleich mit Bild 5, daß 
das Anbringen der Zusatzgewichte an der Kappe sehr viel mehr Einfluß auf die Verschiebung 
der Null- und Unendlichkeitsstellen hat, als das Vergrößern des Kardanringes. Auch hier ver- 
schwindet die untere Unendlichkeitsstelle für genügend großes u,. 

Die durch ihre Breite gefährliche untere Resonanzstelle in Bild 4 kann also durch ge- 
eignete Wahl der Trägheitsmomente fortgeschafft werden. An dem untersuchten Horizont 
wurden die dazu notwendigen rechnerisch ermittelten Zusatzgewichte an der Kappe angebracht. 
Der Versuch zeigte die Richtigkeit der Überlegungen. Der Kreisel ließ sich ohne jede 
Schwierigkeit auf seine Betriebsdrehzahl bringen und zeigte keinerlei anormales Verhalten. 
Der Verlauf der vier Resonanzfunktionen für die nunmehr veränderten Werte der Trägheits- 
momente ist in Bild 7 aufgetragen. Durch Vergleich mit Bild 4 sieht man, welche Ver- 
besserung erzielt ist. Es tritt nur noch die obere sehr schmale Resonanzstelle auf. Aller- 
dings nehmen die Resonanzkurven für kleines y ziemlich große Werte an. Das ist aber 
ungefährlich, da der Kreisel beim Anlaufen dieses Gebiet sehr schnell durchläuft und die 
Drehenergie noch gering ist. Die Versuche haben auch, vor allem beim Auslauf, merkliche 
aber ungefährliche Schwingungen gezeigt. 

Schließlich sollen noch die Grenzwerte der Resonanzfunktionen für große y betrachtet 
werden. An den Gl. (8) kann man diese ablesen. Die Auslenkungen der Kappe (C und D) 
gehen für y>» gegen Null. Die Kappe bleibt also ruhig stehen. Dagegen gehen die beiden 
Resonanzfunktionen A und B gegen den endlichen Grenzwert e, gleichgültig wie groß die 
Verhältniswerte ö, «, und «, sind. Das bedeutet eine Dauerauslenkung um den Wert e, also 
eine Verbiegung der Kreiselachse derart, daß die Hauptträgheitsachse mit der Drehachse zu- 
sammenfällt. Die Selbsteinstellung der elastischen Kreiselwelle wird also durch die Rechnung 
bestätigt. Diese Tatsache beweist noch einmal, weshalb ein Kreiselkörper mit kugelförmigem 
Trägheitsellipsoid ungünstiger ist, als ein solcher mit abgeplattetem Trägheitsellipsoid. Denn 
für einen gegebenen Materialfehler ist der dadurch entstehende Fehlerwinkel e zwischen 
Symmetrieachse und Hauptträgheitsachse für einen Kugelkreisel am größten, es werden also 
auch für ihn die Resonanzfunktionen A und B maximale Werte annehmen. 

Welche Bedingungen zwischen ö, #, und «, ziehen nun das Vorhandensein von nur 
einer reellen Nullstelle des Nenners nach sich? Dazu sei die quadratische Gleichung N =0 
(Gl. (8)) näher betrachtet. Diese hat die Form: 


Die Koeffizienten a und 5 lassen sich sofort aus (8) entnehmen. Aus der Lösung der 
Gl. (11): 


sieht man, daß zwei reelle Wurzelpaare für y nur dann auftreten können, wenn beide Lösungen 
für y? positiv reell werden. Notwendige und hinreichende Bedingung dafür ist 
a<0; b>0. 

Diese Bedingung ist erfüllt für: 

1. ö°<1: das entspricht dem verlängerten Kreisel. 

2.0°>1 +4, + 4, + 4, 4,; das entspricht, je nach den Werten der Verhältnisse u, 

und «,, einem mehr oder weniger stark abgeplattetem Kreisel. 
Um diese Verhältnisse besser überblicken zu können, sollen sie an einem Kurvenbild 

veranschaulicht werden. Bild 8 zeigt die Abszissen der Resonanzstellen yr in Abhängigkeit 
von der Größe ö? für verschiedene Werte von u, aufgetragen. Für ö°< 1 (verlängerter Kreisel) 


3 = 3 
2 2 
7 

— 

15 45 [7 25 

Ir 


Bild 8. Die Resonanzstellen in Abhängigkeit von d? für Bild 9. Die Resonanzstellen in Abhängigkeit von d? für 
verschiedene Werte von verschiedene Werte von 
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treten zwei Wurzeln auf, deren Lage sich mit einer Veränderung von «, nur unwesentlich 
ändert. Die obere Nullstelle liegt jedoch sehr ungünstig, da sie in die Nähe der Betriebs- 
drehzahl kommen kann. Ein verlängerter Kreisel kommt deshalb für die Konstruktion der- 
artiger Kreiselgeräte nicht in Frage. Für den abgeplatteten Kreisel gibt es ein Gebiet, in 
dem nur eine reelle Wurzel auftritt, die ihre Lage sowohl mit ö? als auch mit «, nur wenig 
ändert. Für den stark abgeplatteten Kreisel treten dann wieder zwei reelle Wurzeln auf, 
von denen eine jedoch durch entsprechende Wahl von u, ganz zum Verschwinden gebracht 
werden kann. Der größte Wert, den ö? erreichen kann, ist bekanntlich ö°=4. Das entspricht 
dem ebenen Kreisel. Für den bei den Versuchen benutzten Kreisel (von Anschütz &Co., 
Kiel) war ö6°?=3,43 gemessen worden. 

Das entsprechende Schaubild für verschiedene Werte von u, zeigt Bild 9 (d. h. An- 
bringen der Zusatzgewichte an der Kreiselkappe statt am Kardanring). Wieder haben wir 
für den verlängerten Kreisel zwei sehr ungünstig liegende Wurzeln. Für ö°>1 zunächst nur 
eine reelle Wurzel, die sich mit ö? fast nicht verschiebt, aber durch Vergrößern von u, zu 
kleineren y-Werten verschoben werden kann. Der Fall des sehr stark abgeplatteten Kreisels 
bringt wieder das Auftreten einer weiteren reellen Wurzel, die aber durch Vergrößerung von 
4, fortgeschafft werden kann. 


Aus diesen Betrachtungen kann man ersehen: 
1. Der verlängerte Kreisel ist ungünstig. 
2. Es muß die Bedingung 


erfüllt sein. Da es aus zuvor klargestellten Gründen zweckmäßig ist, 6? möglichst groß zu 
wählen, so muß man durch eine entsprechende Vergrößerung von u die Bedingung (12) er- 
füllen. Bemerkenswert ist, daß die Elastizität der Kreiselwelle in der Bedingung (12) über- 
haupt nicht auftritt. 

Im vorliegenden Fall läßt sich eine sehr anschauliche Deutung der Tatsache geben, 
daß es durch geeignete Wahl der Trägheitsmomente gelingt, eine Resonanzstelle zum Ver- 
schwinden zu bringen. Die untere Resonanzstelle entsteht durch die Übereinstimmung von 
Umlauffrequenz und Nutationsfrequenz des Kreisels. Das soll näher erklärt werden. 

Für die Nutationsfrequenz eines Kreisels mit elastischer Welle gilt für großen Impuls 
die Näherungsformel*): 

—(Oyn+ 9x n) 


Führt man darin die hier benutzten dimensionslosen Bezeichnungen ein, so ergibt sich 
für die bezogene Nutationsfrequenz der Wert: 


n\? y’ 

Trägt man sich über y als Abszisse 
die Frequenz als Ordinate auf, so ergibt 


die Formel (14) für ”- die in Bild 10 ge- 


zeichnete Kurve. Es ist dies im wesent- | “ 
lichen eine Gerade, die für große y-Werte |! 2 13 

etwas abbiegt. Trägt man sich noch dazu 5 v S 7 
die Kurve der Kreiseldrehzahl, also eine & & v 


45 °.Gerade ein, so bekommt man entweder y 
einen Schnittpunkt beider Kurven (Fall a) 
oder gar keinen Schnittpunkt (Fall b). Die 
Abszisse des Schnittpunktes erhält man 


Bild W. Erklärung der unteren Resonanzstelle. 


Yv — 
aus (14), wenn man dort = =y ansetzt: 


(u, + 


Soll kein reeller Schnittpunkt S auftreten, soll also die Umlauffrequenz des Kreisels 
niemals mit der Nutationsfrequenz zusammenfallen, so muß 


4) Siehe Magnus: Schwingungen kraftgekoppelter Kreisel. Ing.-Arch. Bd. 9 (1938), S. 178 bis 204. 


(15). 


| 
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gelten. Das deckt sich aber vollständig mit der zuvor abgeleiteten Bedingung Gl. (12). Im 
Diagramm (Bild 10) bedeutet das, daß die anfängliche Steilheit der Kurve kleiner als 45° 
sein soll, denn nur dann ist bei dem vorliegenden Kurventyp kein Schnittpunkt möglich. 


6. Zusammenfassung. Jeder Kreisel besitzt restliche Unwuchten, die beim Lauf Zentri- 
fugalmomente hervorrufen. Es ist untersucht worden, wie die Trägheitsmomente eines Kreisel- 
gerätes zu wählen sind, damit die durch die Zentrifugalmomente erregten Schwingungen 
möglichst klein bleiben. Damit keine Resonanz von Kreiseldrehung und Nutationsfrequenz 
auftritt, muß zwischen den Trägheitsmomenten des Kreiselkörpers und denen der kardanischen 
Hängung die Beziehung (12) erfüllt sein. Diese Bedingung besagt, daß das Trägheitsmoment 
der Hängung eine bestimmte Mindestgröße haben muß. Ist das nicht der Fall, so kann man 
die Bedingung durch Vergrößerung der Trägheitsmomente der Hängung nachträglich befriedigen. 
Dabei ist ein Vergrößern des Trägheitsmomentes der Kreiselkappe wirksamer, als ein Ver- 
größern des Trägheitsmomentes des Kardanringes. 

Die Resonanzkurve eines Kreiselgerätes mit elastischer Kreiselwelle behält mindestens 
eine Unendlichkeitsstelle, die bei Übereinstimmung der durch die elastische Welle ver- 
ursachten Querschwingung mit der Umlauffrequenz des Kreisels auftritt. Ihre Lage wird 
durch die Wirkung der Kreiselkräfte etwas verschoben. 

Der verlängerte Kreisel ist ungünstig, da für ihn eine meist sehr hoch liegende zweite 
Resonanzstelle auftritt. Aber auch Kreiselkörper mit kugelförmigem Trägheitsellipsoid sind 
unbrauchbar, da für sie die wahrscheinliche Größe des durch Herstellungsfehler hervor- 
gerufenen Fehlerwinkels zwischen Symmetrieachse und Hauptträgheitsachse am größten ist. 
Dieser Fehlerwinkel ist aber für den Lauf oberhalb der kritischen Drehzahl maßgebend. 

Das günstigste Trägheitsmomentenverhältnis für einen symmetrischen Kreiselkörper 
liegt etwa bei Y2 und ist unabhängig von der Elastizität der Kreiselwelle. 148 


Graphische Statik räumlicher Kräftesysteme mit Hilfe der 
dualen Kräfteabbildung. 


Von Robert Sauer in Aachen. 


nter einer projektiven Kräfteabbildung‘) versteht man eine ebene oder räumliche Ab- 

bildung von Kräften, bei der die Wirkungslinien aller Kräfte projektiv transformiert 
werden und jedes Gleichgewichtssystem wieder in ein Gleichgewichtssystem übergeht. In 
der vorliegenden Arbeit wird eine besonders einfache projektive Kräfteabbildung, die duale 
Kräfteabbildung, eingeführt und zur Lösung der Grundaufgaben der graphischen Statik 
räumlicher Kräftesysteme benutzt. Dabei zeigt sich, daß die duale Abbildung die früher von 
R. v. Mises?) und von W. Prager’) angegebenen Abbildungen in einer besonderen Weise 
miteinander verknüpft, indem die Kraftvektoren eines Bündels nach v. Mises und die 
Momentvektoren eines Bündels nach Prager transformiert werden. Durch diese Verkoppelung 
der Misesschen und Pragerschen Abbildungen, die sich zwangsläufig aus dem Begriff der 
dualen Kräftetransformation ergibt, werden die Konstruktionen der graphischen Statik räum- 
licher Kräftesysteme erheblich vereinfacht. 


1. Plückersche Kraftkoordinaten. Wir erinnern zunächst an die bekannte Darstellung 
einer an einem starren Körper angreifenden Kraft (eines „linienflüchtigen“ Vektors) durch 
Plückersche Koordinaten (Bild 1): 

Der in der Wirkungslinie g angreifende Kraftvektor 
f wird bei Annahme eines festen Nullpunktes O in ein- 
deutig bestimmter Weise ersetzt durch den nach O parallel 
verschobenen Kraftvektor f und das Kräftepaar mit dem 
Momentvektor 


den wir ebenfalls von O aus antragen wollen; rt ist der 
Ortsvektor eines beliebigen Punktes A der Wirkungs- 
linie g. Die 2X3 Koordinaten der Vektoren ft, f in einem 
rechtwinkeligen Koordinatensystem mit O als Nullpunkt, 
Bild ı. nämlich k,,k,,k,, k,,k,,k,, nennen wir die Plücker- 


1) R. Sauer: Math. Ann. Bd. 110 (1934), S. 464 bis 472 u. Z. angew. Math. Mech. Bd. 14 (1934), S. 193 bis 198. 
2) Z. Math. u. Phys. Bd. 64 (1917), S. 209 bis 232 u. Z. angew. Math. Mech. Bd. 4 (1924), S. 212. 
3) Z. angew. Math. Mech. Bd. 16 (1926), S. 341 bis 355. 
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schen Koordinaten der gegebenen Kraft und fassen sie im Symbol 


k,,k,,k,, k,, k, 
zusammen. 
Die Plückerschen Koordinaten sind an die aus (1) folgende notwendige Bedingung 


gebunden. Im übrigen aber sind sie willkürlich vorgebbar, d. h. irgend 6 Zahlen, die nicht 
- alle verschwinden und der Nebenbedingungen (2) genügen, bestimmen stets genau eine Kraft 
X als Plückersche Koordinaten. 

Sonderfälle: 

a) B=p,0 (d. h. Momentvektor P=0) stellt eine Kraft dar, deren Wirkungslinie g durch 
den Nullpunkt O geht. 

b) M=0, m (d. h. Kraftvektor m=0) stellt ein Kräftepaar dar mit dem Momentvektor im. 
Als Wirkungslinie betrachten wir in diesem Fall die unendlich ferne Schnittgerade der zu m 
senkrechten Parallelebenenschar. 


2. Definition der dualen Kräfteabbildung. Wir definieren die duale Kräfteabbildung 
durch die Gleichungen 


Sie ordnet jeder Kraft &=f,? mit den Koordinaten k;,k; umkehrbar eindeutig eine Bildkraft 
&=f',f' mit den Koordinaten k’,k’; zu; denn wenn die k,%; Kraftkoordinaten sind und 
daher die Bedingung (2) erfüllen, gilt nach (3) auch 


d.h. auch die k’,%’; sind Kraftkoordinaten. Der Faktor c in (3) ist eine beliebige nichtver- 
schwindende Konstante, die je nach den Maßverhältnissen der vorliegenden Aufgabe zweck- 
mäßig gewählt werden kann. 

Für das Folgende genügt es, die geometrische Bedeutung der dualen Kräfteabbildung 
nur für die beiden Sonderfälle a) und b) der Nr. 1 anzugeben: 

a) Kräfte B=p,0 des Bündels O0. 


Aus $p=0, d.h. p, =p,=Pp,=0 und aus (3) folgt für die Bildkraft P’=p,' pP’: 
Wegen p',=p „=_ liegt der Momentvektor p’ der Bildkraft in der Z-Achse und demnach die 
Wirkungslinie g’ des Bildkraftvektors p’ in der X Y-Ebene. Wir haben damit den Satz: 


Die duale Kräfteabbildung transformiert die Kräfte B=p,0, d.h. 
sämtliche Kräfte des Raums, deren Wirkungslinien g den Punkt O 
enthalten, in die Kräfte ®’=p’,P, deren Wirkungslinien g’ in der 
XY-Ebene liegen (=Kräftefeld der XY-Ebene). 


Außerdem ergibt sich aus (4) folgende Konstruktion der Bildkraft ® (Bild 2): 
Der Grundriß des Kraftvektors p in der X Y-Ebene liefert Größe und Richtung des Bildkraft- 
vektors p’. Die Wirkungslinie g’ des Bildkraftvektors p’ ist parallel zum Grundriß der 
Wirkungslinie 9 der gegebenen Kraft im Abstand 0G@'=ctgp(p=spitzer Neigungswinkel 


von g gegen die X Y-Ebene). Die Vektoren 06, p’, p bilden ein Rechtssystem. 


AN 
Bild 2 
(links). E 

m 
m’ 
Bild 3 
(rechts), 
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Die Wirkungslinie g’ der Bildkraft hat hiernach die Gleichung 


Ausnahmefall: 


Die Kräfte ® mit der Z-Achse als Wirkungslinie (n= =9,=0,9=t 5) bilden sich wegen 


p,ı=p,=p,=0 ab in die Kräftepaare ®’ der X Z-Ebene; die z-Koordinate p, der gegebenen 
Kraft bestimmt das Moment s„=cp, des Bildkräftepaars 

b) Kräftepaare N=0,m. 

Aus m=0, d.h. neu =m,-0 und aus (3) folgt für die Bildkraft M’= m’, m’: 


Wegen m’, = m’,=0 steht der Bildkraftvektor m’ senkrecht auf der X Y-Ebene und wir 
haben damit den Satz: 


Die duale Kräfteabbildung transformiert die Kräfte N=0, m, d.h. 
sämtliche Kräftepaare des Raums, in die Kräfte W=m’,m’, deren 
Wirkungslinien Ah’, senkrecht auf der XY-Ebene stehen (=Parallel- 
kräftebündel senkrecht zur XY-Ebene). 


Außerdem ergibt sich aus (6) folgende Konstruktion der Bildkraft W’ (Bild 3): 
Die z-Komponente des Momentvektors m liefert Größe und Richtung des Bildkraftvektors m’. 
Die Wirkungslinie h’ des Bildkraftvektors m’ ist parallel zur Z-Achse im Abstand OM’=c-ctgg 
(p =spitzer Neigungswinkel von m gegen die X Y-Ebene). O0 M’ ist senkrecht zur Ebene des 


Neigungswinkels und die Vektoren oM, m, m’ bilden ein Rechtssystem. 
Der Spurpunkt M’ der Wirkungslinie h’ der Bildkraft hat hiernach die Koordinaten 


m, m, 
u’ = — 
M m,’ Ym + m, 


M. 
Ausnahmefall: 

Die Kräftepaare M mit den Momentvektoren m in der X Y-Ebene (m, =g9=0) bilden 
sich wegen. m’, =m',=m’,—0 ab in die Kräftepaare M’ mit den zu m parallelen Moment- 
vektoren m’ 

Die Abbildungen eines Vektorbündels O nach a) in ein ebenes Vektorfeld und nach b) 
in ein Parallelvektorbündel wurden schon von R. v. Mises?) und von W. Prager‘) benutzt. 
Wir werden darauf in Nr. 4 zurückkommen. 


3. Eigenschaften der dualen Kräfteabbildung. Ehe wir zeigen, wie sich die duale Kräfte- 
abbildung für die Konstruktionen der graphischen Statik räumlicher Kräftesysteme verwerten 
läßt, stellen wir hier die wichtigsten Eigenschaften der dualen Kräfteabbildung zusammen. 

Die grundlegende mechanische Eigenschaft, auf der die praktische Verwendbarkeit 
der dualen Kräfteabbildung beruht, bringt folgender Satz zum Ausdruck: 


Jedem Gleichgewichtssystem ..., eines starren Körpers 
entspricht bei der dualen Kräfteabbildung wieder ein Gleich- 
gewichtssystem --.,®’y und umgekehrt. 


Beweis: Die Gleichgewichtsbedingungen 
des gegebenen Kräftesystems werden durch (3) in die Gleichgewichtsbedingungen 


des Bildkräftesystems transformiert. 
Die grundlegende geometrische Eigenschaft der dualen Kräfteabbildung läßt sich 
folgendermaßen aussprechen: 


Die Wirkungslinien g der gegebenen Kräfte und die Wirkungs- 
linien g’ der Bildkräfte stehen in einer dualen projektiven Beziehung, 
d.h. den Geraden g eines Bündels entsprechen stets die Geraden y’ 
einer Ebene und den Geraden g einer Ebene entsprechen stets die 
Geraden g’ eines Bündels. Die Punktkoordinaten x,y,2 der Scheitel 
der 9-Bündel und die Ebenenkoordinaten w’, v’, m’, der zugeordneten 
g’-Ebenen sind verknüpft durch die Gleichungen 

1 


c2’ c2’ c2' 
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Da dieser Satz nur theoretisches Interesse hat und im folgenden nicht benutzt wird, 
wollen wir auf die Angabe des Beweises, der einige liniengeometrische Vorkenntnisse erfordert, 
verzichten. Wir beschränken uns darauf, den Satz an den beiden Sonderfällen a) und b) von 
Nr. 2 zu bestätigen: 

a) Das Geradenbündel mit dem Scheitel O wird in die X Y-Ebene abgebildet. 

b) Die Geraden der unendlich fernen Ebene (= Wirkungslinien der Kräftepaare des 
Raums) werden in das Parallelgeradenbündel mit dem unendlich fernen Punkt der 
Z-Achse als Scheitel transformiert. 

Wie wir eingangs erwähnten, bezeichnet man Kräfteabbildungen, bei denen die Wirkungs- 
linien projektiv transformiert werden und die Gleichgewichtsbedingungen erhalten bleiben, 
als projektive Kräfteabbildungen. Aus den beiden vorhergehenden Sätzen folgt: 

Die duale Kräfteabbildung ist eine spezielle projektive Kräfte- 
abbildung. 


4. Verwendung der dualen Kräfteabbildung in der graphischen Statik. Jedes an einem 
starren Körper angreifende Kräftesystem kann bekanntlich ersetzt werden durch eine am 
Nullpunkt O angreifende Einzelkraft B=p,0 und ein Kräftepaar M—0,m. Nach den Sonder- 
fällen a) und b) von Nr. 2 werden bei der dualen Kräfteabbildung die in O angreifenden 
Kräfte ® in Bildkräfte ®P’=p’,p’ der XY-Ebene und die Kräftepaare M in Bildkräfte 
M’ — m’, m’ senkrecht zur X Y-Ebene transformiert. 
Wir erhalten daher in der X Y-Ebene als Bild- 
ebene folgende Abbildung eines beliebigen 
räumlichen Kräftesystems (Bild 4): m’ 
a) einen in einer Geraden g’ der Zeichenebene 

b) einen Vektor m’, der in einer zur Zeichen- 
ebene senkrechten Geraden linienflüchtig 
ist und durch seinen Grundriß M’ sowie 
durch einen Aufriß oder eine Umklappung 
gegeben werden kann. 

Die 4 das Kräftesystem ®, M bestimmenden 
Elemente g’, p’, M’, m’ bezeichnen wir fortan kurz ans 
als Bild des gegebenen Kräftesystems, g’ Bild 4. 
nennen wir gelegentlich das Bild der Geraden g. 

Aus der Erhaltung der Gleichgewichtsbedingungen (vgl. Nr. 3) folgt: 

Alle Zusammensetzungen bzw. Zerlegungen räumlicher Kräfte- 
systeme lassen sich in der Zeichenebene durchführen, indem man 
mit den Hilfsmitteln der ebenen graphischen Statik einerseits die 
Vektoren p’ und andererseits die Vektoren m’ zusammensetzt bzw. 
zerlegt. 

Die Abbildungen eines Vektorbündels O a) in die Vektoren p’ einer Bildebene oder b) in 
die Vektoren m’ senkrecht zu einer Bildebene wurden schon in Arbeiten von R. v. Mises?) 
und W. Prager°) benutzt. Während es sich dort aber zunächst nur um Kräftesysteme eines 
Bündels handelte, können wir hier allgemeine räumliche Kräftesysteme konstruktiv erfassen. 
Dabei führt die duale Kräfteabbildung zwangsläufig dazu, daß wir die Kraftvektoren nach 
v. Mises gemäß a) und die Momentvektoren nach Prager gemäß b) abbilden. Diese sich 
hier naturgemäß ergebende Verkoppelung der Misesschen und Pragerschen Abbildungen 
wird in Nr. 5 und 6 durch einfache und durchsichtige Konstruktionen für die graphische Statik 
allgemeiner räumlicher Kräftesysteme auch ihre praktische Rechtfertigung finden. 


5. Abbildung von Einzelkräften und Kraftschrauben. Wir wenden jetzt das Abbildungs- 
verfahren von Nr. 4 auf Einzelkräfte und Kraftschrauben an. 

I. Abbildung einer Einzelkraft. 

Wenn das Kräftesystem aus einer Einzelkraft &=1,f besteht, sind bei der Zerlegung 
nach Nr. 4 in B=p,0 und NM=0, m die Vektoren p, m mit den Vektoren f£,f von Nr. 1 (vgl. 
Bild 1) identisch. Nach (2) gilt also 


+p,m,+ Mm, 
und weiter mit Hilfe von (7) 


Durch Vergleich mit (5) folgt der Satz: 
Im Bild g’‘, p’, M’, m’ liegt der Punkt M’ dann und nur dann auf der 
Geraden gg’, wenn das gegebene räumliche Kräftesystem sich auf eine 
Einzelkraft reduziert. 


12 
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In Bild 4 ist ein nicht auf eine Einzelkraft reduzierbares Kräftesystem, in Bild 5 da- 
gegen eine Einzelkraft dargestellt. 


Bild 6 (rechts). 


Bestimmung der Kraft $ aus ihrem Bild g’,p', M’,m’: 

a) Der Kraftvektor p und damit die Richtung der Wirkungslinie g ist durch g’ und p’ 
allein festgelegt; die diesbezügliche Konstruktion wurde schon in Nr. 2 (vgl. Bild 2) behandelt. 
P) Der Spurpunkt S der Wirkungslinie g ergibt sich aus g’,p’, M’,m’ nach folgender 
Konstruktionsvorschrift (Bild 6): 


Man ziehe den Vektor q’= OF’ von der Länge ja senkrecht zu g’ und zwar so, daß 


’ 


q,p’, m’ ein Rechtssystem bilden. Der Spurpunkt S der Wirkungslinie g ist dann der Schnitt- 
punkt der Geraden OM’ mit der Parallelen zu g’ durch F’. 


Beweis: Nach (1) ist m—=rxXp, wobei jetzt v den Ortsvektor OS des Spurpunktes S mit 
den Koordinaten #,, 9; (2,0) bedeuten soll. Daraus folgen die 3 Komponentengleichungen 
Ps — Ps, X, Ps — YsPı- 

Wegen 2,—0 erhält man = 
m, m, 
=— 
Ps 
also nach (7) 
Y— 
d. hı. der Spurpunkt S liegt auf der Geraden O M’. Ferner ist die --Komponente des Moment- 
vektors m der gegebenen Kraft einerseits definitionsgemäß gleich q’ X p’, andererseits auf Grund 
von (6) wegen m’, — m’,—0, m’, = m, gleich m’; daher bilden g’, p’, m’ ein Rechtssystem und, 
da p’ und q’ außerdem senkrecht sind, hat man 
Sonderfälle und Ausnahmefälle: 
a) Wenn die Wirkungslinie der Kraft durch den Nullpunkt O geht, wird m’=0 und M’ 
unbestimmt, das Bild der Kraft reduziert sich auf g’, p’. 
b) Wenn die Kraft in ein Kräftepaar entartet, wird p'=0 und g’ unbestimmt; das Bild 
der Kraft reduziert sich auf M’, m’. 
c) Wenn die Wirkungslinie der Kraft die Z-Achse schneidet, entartet nach Nr. 2 m’ in 
ein Kräftepaar; im Bild der Kraft ist M’ der unendlich ferne Punkt der Geraden g’. 


d) Wenn die Wirkungslinie der Kraft auf der Bild- N 
ebene senkrecht steht, entarten nach Nr. 2 m’ und p’ H 
in Kräftepaare; im Bild der Kraft rücken M’ und g’ Im} 


ins Unendliche. 

Um nicht zu weitläufig zu werden, wollen wir von 
nun an die beiden Ausnahmefälle c) und d) von der Be- 
trachtung ausschließen. 

HU. Abbildung zweier Kräfte. 

Zwei Kräfte 8, und St}, seien vorgegeben durch ihre 
Der Fallunterscheidung nach parallelen, sich schneiden- 
den oder windschiefen Wirkungslinien g;, 9 der ge- 
- gebenen Kräfte $,, 8, entspricht folgende Fallunter- 
scheidung im Bild der Kräfte: 

Wenn g,=gyn Sind 9, parallel. 
Wenn g’; verschieden ist von g’y, setze 
man und m’, m’; zum resultieren- 
den Vektor p’5; bzw. m’; zusammen. Bild 7. 


XUM 
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Falls dann der Grundrißpunkt M’, von m’; auf der Angriffslinie g' 
von p’z liegt, schneiden sich die Wirkungslinien 9, 9n und &, 
lassen sich zur resultierenden Kraft 8, mit dem Bild g’r,Pp’r, M’'r, m’r 
zusammensetzen; liegt dagegen M’z nicht auf g’z, dann sind die Wir- 
kungslinien 9, 91 windschief. 


II. Abbildung einer Kraftschraube. Eine Kraft 8&=f,f und ein Kräftepaar 
R=0,T faßt man bekanntlich unter der Bezeichnung Kraftschraube zusammen, wenn 
der Kraftvektor £ von & und der Momentvektor t von R parallel sind. 

Außerdem schicken wir folgende Definition voraus (vgl. 

Bild 8): Ein Punkt R’ der Bildebene heißt Antipol einer Ge- 

raden g’ der Bildebene, wenn die Verlängerung des Lotes 06’ 

e » zu g' über O hinaus den Punkt R’ enthält und die Strecken 0@’, 
OR’ in der Beziehung 0@'- OR’=c? stehen. 

Für das Bild einer Kraftschraube gilt dann folgender Satz 
(Bild 8): 

Eine Kraft & mit dem Bild g, p’, M’, m’ und ein 
Kräftepaar A mit dem Bild X’, r’ bilden eine Kraft- 
schraube, wenn R’ der Antipol von g’ ist. 

Beweis: Nach (7) hat R’ die Koordinaten 


r, 


und nach (5) wird g’ vom Lot aus O im Punkt @’ mit den Koordinaten 


Pı Ps 
Pr 


geschnitten. Da p und T parallel sind, also 
r, 


gilt, folgt aus den vorhergehenden Gleichungen 


Ps Ps + 

Yp = c = — Yır 

Ps 9a 96 + 


womit die Behauptung bewiesen ist. 


6. Grundaufgaben der graphischen Statik räumlicher Kräftesysteme. Zum Schlusse zeigen 
wir, wie man mit Hilfe der dualen Kräfteabbildung zwei Grundaufgaben der graphischen 
Statik räumlicher Kräftesysteme lösen kann. 


Aufgabe I: Vorgegeben ist das Bild g’, p, M’, m’ eines Kräftesystems, das 
sich nicht auf eine Einzelkraft und nicht auf ein Kräftepaar 
reduzieren soll. Gesucht ist k 
das Bild der dasKräftesystem 
ersetzenden Kraftschraube. Ya 
Konstruktion (Bild 9): 
Man fälle das Lot 0@’ aus O auf g’ und ermittle # ee SE 
auf der Verlängerung von 0G’ über 0 den Antipol R’ ,# 
der Geraden g’ (vgl. Nr. 5, II). Dann bestimme man 
den Schnittpunkt M’, der Geraden g’ und R’M’ und 
zerlege den Vektor ın’ in die parallelen, in R’ und M’, 
angreifenden Vektoren r’ und m’,. Damit ist das 
Kräftesystem g’, p’, M’, m’ zerlegt in die Kraft mit 
dem Bild 9’, (= pP’, M’,, m’, und das Kräfte- 
paar mit dem Bild R’, v’; nach Nr. 5, III bildet die 
Kraft g’,, p’,, M’,, m’, mit dem Kräftepaar R’, v’ die 
gesuchte Kraftschraube. Bild 9. 
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Aufgabe Il: Vorgegeben ist das Bild g’, p', M’ m’ eines Kräftesystems, das 
sich nicht auf eine Einzelkraft und nicht auf ein Kräftepaar 
reduzieren soll, sowie der Spurpunkt U und das Bild a’ einer 


Geraden a. 


Gesucht ist das Bild des das Kräftesystem er- 


setzenden Kraftkreuzes, dessen eine Wirkungslinie die Ge- 


rade a ist. 
Konstruktion (Bild 10): 


Wir wählen den Spurpunkt U als Nullpunkt O der 
dualen Kräfteabbildung; 9’, p’, M’, m’ bzw. a’ sei das Bild 
des gegebenen Kräftesystems bzw. der Wirkungslinie a in 


bezug auf diesen Nullpunkt. 
Schnittpunkt P’ von g’ und a’ und den Punkt M’ die Gerade b’ 


Dann ziehe man durch den 


und zerlege den Vektor p’ nach den Geraden a’ und b’ in die 


M= My 
Bild 10. 


Komponenten p’,; und 
g', p’, M’, m’ zerlegt in die Kraft mit dem Bild a’, p’,, 
deren Wirkungslinie a durch O= U geht, und in die Kraft 
mit dem Bild d’, M’„(= M'), m’, (= m’); die Kräfte a’, p’, 
und m’j bilden das gesuchte Kraftkreuz. 

Die Konstruktion versagt, wenn M’ auf a’ liegt; in 
diesem Falle ist a eine Nullinie des durch das gegebene 
Kräftesystem bestimmten Nullsystems. 129 


Damit ist das Kräftesystem 


BUCHBESPRECHUNGEN 


FRITZ WEINIG, Aerodynamik der Luft- 
schraube. XVI-+484S. ın. 374 Abb. Berlin 1940, 
Verlag Julius Springer. Preis geb. 48 M. 


Die eingehende Behandlung eines solchen Spezial- 
zebietes, wie es die Aerodynamik der Luftschraube 
darstellt, kommt einem schon lange gehegten 
Wunsche entgegen. Zwar gab es viele und gute 
Sonderarbeiten, auch zusammenfassende Darstel- 
lungen, aber die ersteren blieben meist ohne Zu- 
sammenhang miteinander, und die letzteren stam- 
men größtenteils aus älterer Zeit, so daß wichtige 
Teile der neueren Entwicklung fehlen. Das vor- 
liegende Buch von Weinig füllt diese Lücke aus 
und erweitert dazu unsere Kenntnis von der Luft- 
schraube nach verschiedenen Richtungen hin sehr 
erheblich. 

Die Aufgaben und Fragen, die ein Lehrbuch 
dieser Art aufzuzeigen hat, sind: 

Erkenntnis der theoretischen Grundlagen. 

Gegebenheiten der praktischen Ausführungen. 

Wie ist eine Schraube zu berechnen? 

Wie sind die besten Verhältnisse zu finden? 

Wie ist eine vorgelegte Schraube zu beurteilen? 

Wie verhält sich ein Propeller unter verschie- 
denen Betriebsbedingungen? 

Wie kann die bisherige Theorie verfeinert 
werden? 

Über alle diese Dinge gibt das Buch erschöp- 
fende Auskunft. Dankenswert ist schon im Eingang 
die sehr große Zahl von Zahlenangaben und von 
Verfahren für die Aufmessung und Ermittlung fest- 
stehender geometrischer und physikalischer Daten 
bei vorgelegten Luftschrauben (Gewichte, Trägheits- 
momente, Abmaße ausgeführter Schrauben usw.). 

Die früher üblichen verschiedenartigen Betrach- 
tungsweisen für die Propeller als Strahltheorie, 
Flügelblattheorie, Wirbeltheorie waren nur zum 
Teil durch sachliche, dem Gegenstand angepaßte 
Gesichtspunkte entstanden, vielmehr jedoch aus 
. der historischen Entwicklung heraus bedingt. Für 
die verschiedensten Zwecke wurde bald die eine, 
bald die andere Lehre benutzt, man kombinierte sie 
auch wohl für bestimmte Rechenaufgaben (Kärmän, 
Bienen, Glauert u. a... In dem Buche wird dies 
ersetzt durch eine Darstellung, die dem Problem 
viel mehr und unmittelbarer angepaßt ist: Zu- 


nächst „die Elementarschraube“: sie ist für 
schwache und große Belastung exakt berechenbar, 
die Reibung läßt sich mitberücksichtigen, und vor 
allem lehrt diese Theorie den Einfluß von Stei- 
gung und Blattbreite wie Leistung auch unter ver- 
schiedenen Betriebsverhältnissen vollkommen exakt 
beurteilen. Man kann diese Lehre von der Ele- 
mentarschraube als die wichtigste und vielfach 
schon genügende Annäherung an die genaue 
Schraubentheorie ansprechen. 

Viele wertvolle Erkenntnisse (auch aus eigenen 
Arbeiten des Verfassers) sind in diesem ersten 
großen Abschnitt untergebracht. Bemerkenswert 
erscheint unter anderen die exakte Darstellung 
von Schub- und Drehmoment in Abhängigkeit vom 
Fortschrittsgrad 4. Die danach für kleine Stei- 
ungsverhältnisse vereinfachten Formeln und 
iagramme bedeuten eine Rechtfertigung für die 
von verschiedenen Autoren als erste Näherung an- 
genommene lineare bzw. quadratische Abhängigkeit 
dieser Größen von 4. 

Interessant ist auch das eingeschobene Kapitel 
über Leitvorrichtungen, weil es die Theorie gegen- 
läufiger Schrauben in sich schließt. Dabei werden 
auch sonst selten vorkommende, aber für die zu- 
künftige Entwicklung (für sehr schnelle Flugzeuge) 
wichtige Erkenntnisse gewonnen und ausdrücklich 
herausgestellt. Auch die schwierigen Fälle un- 
gleichmäßiger Schraubenströmungen (z. B. Schräg- 
anblasen) behandelt das Buch in seltener Ausführ- 
lichkeit. 

Die folgerichtige Weiterentwicklung führt dann 
zur „Gesamtschraube“ (und zwar: 1. schwach- 
belastete Schraube mit sehr großer Flügelzalıl, 
2. Luftschraube mit endlicher Flügelzahl), die durch 
radiale Zusammenziehung des Schraubenstrahles 
und durch Druckänderung in demselben von der 
Elementarschraube unterschieden ist. Hier ist das 
wichtige Froblem der günstigsten Schubverteilung 
zuerst näherungsweise und dann als exakte Rand- 
wertaufgabe durchgeführt. 

In diesem, wie auch in den anderen Kapiteln 
ist die Darstellung so gehalten, daß nach einem 
kurzen Überblick über die zu erwartenden Resultate, 
die jedesmal hinzugenommenen Erweiterungen und 
Verfeinerungen der Theorie (Reibung, endliche 
Flügelzahl usw.) Schritt für Schritt nacheinander 
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entwickelt werden. Zu allen genannten theore- 
tischen Entwicklungen bringt das Buch auch 
Zahlenwerte und stellt die Ergebnisse in Kurven 
fest, die für die praktische Berechnung sofort her- 
angezogen werden können. In einem weiteren be- 
sonderen Abschnitt „Berechnung von Luftschrau- 
ben“ wird gezeigt, wie der Berechnungsgang nor- 
malerweise durchgeführt werden kann, entweder 
auf Grund der benützten Tragflügelprofile oder 
nach der Methode der „gleichwertigen Profil- 
polaren“, 

Im letzten Teil „Zusammenwirken von Luft- 
schraube und Flugzeug“ werden die Grundlagen 
für die Behandlung dieser nicht ganz einfachen Zu- 
sammenhänge dargelegt und durch Zahlenrech- 
nungen und Vergleiche mit Versuchsergebnissen 
erläutert. 

Sehr dankenswert ist die Beigabe einer großen 
Zahl von Versuchsergebnissen mit Luftschrauben 
und ein Schrifttumsverzeichnis von sehr erheb- 
lichem Umfang (296 Arbeiten!). 

Dagegen ist die konstruktive Ausführung 
von Luftschrauben in dem Buche nicht enthalten. 
Zwar rechtfertigt der Titel des Werkes diesen 
Verzicht, aber der Verfasser gibt im Vorwort auch 
noch einen besonderen Grund dafür an, daß näm- 
lich gerade die konstruktive Gestaltung in sehr 
schneller Entwicklung begriffen ist. Eine Darstel- 
lung in dem Buche würde es (wenigstens nach 
dieser Richtung hin) schnell veralten lassen. 

Die hier gegebene kurze Übersicht zeigt schon, 
wie umfangreich das Wissensmaterial ist, das der 
Verfasser in diesem Buche zusammengetragen hat. 
Es ist nicht überall leicht zu lesen; das liegt aber 
an der Vielseitigkeit und den zum Teil sehr ver- 
wickelten Zusammenhängen des Stoffes. Dem 
Verfasser ist zuzugestehen, daß er es gut verstan- 
den hat, auch solche schwierige Stellen anschaulich 
zu gestalten, und daß er überhaupt in allen Teilen 
die Darstellung flüssig und klar wiedergegeben 
hat. Das Werk von Weinig bedeutet jedenfalls 
einen sehr wertvollen Zuwachs in der neueren 
aerodynamischen Literatur. 


Hannover. A. Pröll. 168 


H. W. DROSTE, Die Lösung angewandter 
Differentialgleichungen mittels La- 
placescher Transformation. Mit einem 
Vorwort von Prof. Dr. G. Doetsch. .‚Neuere 
Rechenverfahren der Technik“. Heft 1. 35 S. Berlin 
1939, Verlag E. S. Mittler u. Sohn. Preis 5 M. 


Das algebraische Rechnen mit dem Differential- 
symbol, das bis auf Leibniz zurückgeht, wurde 
von dem Ingenieur Oliver Heaviside mit 
eroßer Virtuosität für die Lösung von Aufgaben 
ausgebaut, die auf partielle Differentialgleichungen 
führen. Da Heaviside seine Entwicklungen rein 
formal durchführt, ohne eine eigentliche mathema- 
tische Begründung zu geben, bedarf ihr Gebrauch 
eroßer Vorsicht, wenn man falsche oder unvoll- 
ständige Ergebnisse vermeiden will. Schon vor 
einer Reihe von Jahren ist die Erkenntnis aufge- 
taucht, daß der Heavisidesche Formalismus auf 
eine Funktionentransformation hinauskommt, die in 
den Arbeiten von G. Giorgi, K.W. Wagner, 
J.R. Carson, G. Doetsch und B. van der 
Pol in verschiedener Form eingeführt wurde. G. 
Doetsch insbesondere verdanken wir eine kon- 
sequent auf die Laplacesche Transformation auf- 
gebaute Darlegung des Öperatorenkalküls vom 
Standpunkt der modernen Mathematik. 

H. W. Droste erläutert in dem vorliegenden Heft 
in aller Kürze die Lösung von Differentialgleichun- 
gen mittels der Funktionentransformation nach La - 
place in engem Anschluß an die Darstellung von 
G. Doetsch. Er leitet die wichtigsten Regeln 
für das Rechnen mit der Laplaceschen Transfor- 


mation her und bringt Beispiele für ihre Anwen- 
dung. Es folgt die Lösung von gewöhnlichen line- 
aren Differentialgleichungen mit konstanten Koef- 
fizienten und der Heavisidesche Entwicklungssatz, 
der hier allerdings nur für rationale Operatoren be- 
wiesen wird, während seine wichtigsten technischen 
Anwendungen doch auf dem Gebiet der mero- 
morphen Funktionen liegen. Im nächsten Abschnitt 
zeigt der Verfasser im Anschluß an Arbeiten von 
B. van der Pol und H. Fischer, wie man 
durch die Laplacesche Transformation der Bessel- 
schen Differentialgleichung auf einfachem Wege 
wichtige Beziehungen zwischen den Besselschen 
Funktionen erhält. Auch bei den linearen partiellen 
Differentialgleichungen bietet die Anwendung der 
Laplaceschen Transformation wesentliche Vorteile, 
was hier am Beispiel der Telegraphengleichung ge- 
zeigt wird. Auf die Lösung durch asymptotische 
Reihen, für die seinerzeit Heaviside zahlreiche 
Beispiele gegeben hat, geht der Verfasser nicht ein, 
weil dies über den von ihm gesteckten mathema- 
tischen Rahmen hinausgeht. Eine Sammlung von 
Laplace-Integralen am Schluß des Heftes rundet 
die Darstellung ab. 


Die vorstehende Inhaltsübersicht zeigt, daß der 
Verfasser auf wenig Seiten ein ausgedehntes Sach- 
gebiet behandelt und die wichtigsten Grundzüge 
einer Methode erläutert, die für den theoretisch ar- 
beitenden Ingenieur und Physiker von großem 
Wert und Nutzen ist. In ihrer sehr konzentrierten 
Form ist diese Darstellung keine leichte Lektüre, 
obgleich an mathematischen Vorkenntnissen nicht 
mehr vorausgesetzt wird, als was ein Ingenieur im 
Diplomexamen nachweisen muß. 


Berlin. Karl Willy Wagner. 149 


Dr. HERMANN DIESSELHORST, Prof. em. d. 
Physik a. d. Techn. Hochsch. in Braunschweig, 
Magnetische Felder und Kräfte mit 
einer Übersicht über die Vektoren- 
rechnung, eine mit Zusätzen versehene Sonder- 
ausgabe des Beitrages „Elektrodynamik“ aus Hand- 
buch der Elektrizität und des Magnetismus, Bd. IV 
(1920). VI+ 216 S. m. 56 Abb. Leipzig 1939, Ver- 
lag Johann Ambrosius Barth. Preis kart. 8 M. 


In einem ersten kurzen Abschnitt werden die 
experimentellen Grundlagen gegeben, wobei der 
Verfasser besonders auf die Gedankengänge von 
Ampere eingeht, für weitere Versuche dagegen 
auf die bekannten ausführlichen Physikbücher ver- 
weist. Den Hauptinhalt bilden auf 140 Seiten die 
Systeme von Magneten und elektrischen Strömen 
unter starker Anwendung der Vektoren- und 
Affinoren-Rechnung: die mechanischen Kräfte und 
Energieverhältnisse werden ausführlich behandelt. 
In einem Geleitwort betont Fritz Emde, daß die 
Vertreter der Mathematik bei der Vektorenrech- 
nung meist nur die Geometrie und Mechanik be- 
rücksichtigen, während ihnen Elektrizitätslehre 
ferner liegt. Und doch lassen sich viele Probleme 
der Elektrodynamik, z. B. die sprunghaften Ände- 
rungen an den Körperoberflächen in einfacher 
Weise nur unter Zuhilfenahme der Vektorenrech- 
nung behandeln. Aus diesem Grunde ist ein län- 
gerer Abschnitt von fast 50 Seiten über diese 
Rechnungsmethode vorausgeschickt. Da die wei- 
teren Fortschritte der Wissenschaft an den hier 
gegebenen Grundlagen nicht viel ändern werden, 
so dürfte dieses Buch für lange Zeit allen denen, 
die sich mit magnetischen Fragen zu befassen 
haben, eine wertvolle Hilfe sein. Das Buch ist für 
denjenigen, der mit der Vektorenrechnung ver- 
traut ist, klar und leicht verständlich geschrieben, 
setzt allerdings eine Beherrschung der mathema- 
tischen Denkweise voraus. 


Freiberg (Sachsen). Brion. 16 
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ERWIN LOHR, Prof. a. d. Deutschen Techn. 
Hochschule Brünn, Vektor- und Dyaden- 
reehnung für Physiker und Tech- 
niker. („Arbeitsmethoden der modernen Natur- 
wissenschaften.“ XV +411 S. m. 34 Fig. Berlin 
1939, Verlag Walter de Gruyter & Co. Preis geb. 
18 M. 


Als erster Band einer neuen Sammlung „Arbeits- 
methoden der modernen Naturwissenschaften“ 
erscheint das vorliegende Buch, das vor allem 
der Förderung der praktischen Verwendung der 
Vektor- und Dyadenrechnung in Physik und 
Technik dienen will. Offenbar der Brünner Tra- 
dition folgend, verwendet Verf. die Gibbs-Jau- 
inannsche Schreibweise. Er ordnet den Stoff in 
drei Abschnitte. In den beiden ersten: Arithmetik 
und Algebra extensiver Größen und Analysis 
dieser Größen entwickelt er zunächst den ganzen 
Formelapparat. Er knüpft dabei an das Rechnen 
mit reellen und komplexen Größen an und sucht 
durch das Verfahren des „Projizierens“ die syste- 
matische Verbindung mit der gewöhnlichen Arith- 
metik, Algebra und Analysis möglichst aufrecht 
zu halten. Der Hauptabschnitt bringt dann An- 
wendungen insbesondere auf verschiedene Gebiete 
der Physik, so auf das der Mechanik, der Elasti- 
zitätstheorie, der Hydrodynamik, -der elektro- 
magnetischen und optischen Erscheinungen, 
ferner einen Ausblick auf die Verwendung in der 
Quantentheorie. Ein kurzer Abschnitt ist den 
Anwendungen in der Geometrie gewidmet. 

Das Buch ist im wesentlichen auf die Bedürf- 
nisse des theoretischen Physikers zugeschnitten. 
Als Einführung in das Gebiet eignet es sich, 
wenigstens für den Techniker, kaum, einmal 
weil zunächst ein Übermaß von Formeln zu be- 
wältigen ist, deren Verwendungsmöglichkeit meist 
erst viel später aufgezeigt wird. dann aber 
auch, weil die Anschauung ziemlich zurücktritt; 
dagegen wird es dem, der schon einigermaßen mit 
dem Gegenstand vertraut ist, zur Weiterarbeit 
und auf Grund des ausführlichen Inhaltsverzeich- 
nisses zum Nachschlagen gute Dienste leisten. 


Dresden. Willers. 115 


A. ADRIAN ALBERT, Associate Professor a. d. 
Universität Chicago, Strueture of Algebras. 
(American Mathematical Society. Colloquium Publi- 
eations, Vol. XXIV.) XI+ 210 S. New York 1939, 
American Mathematical Society. Preis 4 Dollar. 


Man pflegt einen Größenbereich, in dem Addition, 
Subtraktion und Multiplikation unbeschränkt aus- 
führbar sind, d. h. einen Ring, dann als Algebra zu 
bezeichnen, wenn er über einem Körper liegt, mit 
dessen Elementen auch im erweiterten Bereich un- 
beschränkt und in gewöhnlicher Weise gerechnet 
werden kann. Die reiche Entwicklung, die die 
Theorie der Algebren in den letzten Jahren genom- 
men hat, ist eine schöne Frucht internationaler Zu- 
sammenarbeit. Dem Verfasser selbst sind große 
Fortschritte dabei zu verdanken, er ist daher wie 
kein anderer berufen, eine Gesamtdarstellung zu 
geben, die dreiteilig geplant ist. Die Einführung in 
das Gebiet brachte das Buch Modern Higher Al- 
gebra, hier liegt der zweite Teil vor, ein dritter 
über die Arithmetik der Algebren soll folgen. Das 
Buch bringt die Fortschritte, die seit dem Bericht 
Deurings (Ergebnisse der Mathematik und ihrer 
Grenzgebiete IV, 1 1935) zu verzeichnen sind, und 
gibt bei älteren Untersuchungen neue Darstellun- 
gen und wesentliche Vereinfachungen. Vorwort 
und Inhaltsverzeichnis unterrichten ausführlich über 
er ru auf den einzugehen hier zu weit führen 
würde. 


Halle (Saale). H. Brandt. 150 


Dr. WOLFGANG KRULL, o. Prof. a. d. Universi- 
tät Bonn, Elementare Algebra vom 
höheren Standpunkt. (Sammlung Göschen, 
Bd. 930.) 142 S. m. 6 Zeichn. Berlin 1939, Verlag 
Walter de Gruyter & Co. Preis geb. 1,62 M. 


Das reichhaltige Büchlein bringt die elementare 
Algebra in dem üblichen Umfang unter Einschluß 
der genäherten Wurzelberechnung bis zur Einfüh- 
rung in die Galoissche Theorie, die da, wo sie auch 
historisch begonnen hat, nämlich an den Kreis- 
teilungsgleichungen, erläutert wird. In 6 Kapiteln 
werden behandelt: Formales Buchstabenrechnen, 
Nullstellen und Zerlegung von Polynomen, Auflösung 
der Gleichungen ersten bis vierten Grades, höhere 
Gleichungstheorie, Kreisteilungstheorie und zuletzt 
numerische Wurzelberechnung. Dabei wird überall 
von modernen Begriffen und Gesichtspunkten Ge- 
brauch gemacht, und es werden Hinweise auf 
weiterführende Untersuchungen gegeben. So wird 
derjenige, der Bekanntes wiederholen will, durch 
neue Ausblicke angeregt, während der Weiter- 
strebende zugleich eine erste Einführung in neuere 
algebraische Untersuchungen findet. 


Halle (Saale). H. Brandt. 151 


L’Istituto Nazionale per le Applicazioni del 
Caleolo, Nel Quadriennio 28 Ottobre 1933 — — 
27 Ottobre 1937—XV, PublicazionidellIsti- 
tuto perle Applicazioni del Calcolo, 
N 27. VIII+133+X Seiten. Roma 1938—XVI, 
Consiglio nazionale delle Ricerche. 


Dieser Tätigkeitsbericht des römischen Rechen- 
instituts enthält zunächst einen kurzen Überblick 
über die Tätigkeit in den Jahren 1933 bis 1937, 
eine Schilderung der Organisation des Instituts, 
Urteile über das Institut, Personalangaben (S. 1 bis 
38); eine Liste der Veröffentlichungen der Insti- 
tutsmitglieder (S. 39 bis 48), sodann 118 Berichte 
über die einzelnen Arbeiten (S. 51 bis 116), wo bei 
einigen kriegswichtigen Arbeiten allerdings nur 
die Titel angegeben werden; Angaben über die Zu- 
sammenarbeit mit andern Anstalten (S. 119 bis 
129); ein Namenverzeichnis (S. 131 bis 133) und 
eine Liste früherer Veröffentlichungen aus den 
Jahren 1927 bis 1936 (S. I bis X). 

Das Recheninstitut hat 14 Mitarbeiter und 7 Mit- 
arbeiterinnen. Die meisten davon sind Mathe- 
ınatiker. Direktor ist Prof. M. Picone. Das Institut 
ist aufs reichste mit graphischen und mechanischen 
Instrumenten ausgestattet; doch werden Einzel- 
heiten darüber nicht angegeben. Das Institut 
nimmt sowohl von den Behörden wie von der In- 
dustrie Aufträge entgegen. Die Arbeiten sind 
außerordentlich vielseitig. sie betreffen fast alle 
Gebiete der Mechanik, Physik, Technik und Sta- 
tistik. Es wird berichtet, daß auch neue Berech- 
nungsmethoden entwickelt worden sind. 

Wer das Buch kostenlos zugesandt zu erhalten 
wünscht, wende sich unmittelbar an das Istituto 
Nazionale per le applicazioni del calcolo, Roma, 
Piazzale delle Scienze. 


Stuttgart. Fritz Emde. 158 


Prof. Dr. J. PETERS, Sechsstellige Ta- 
feln der trigonometrischen Funk- 
tionen. 2. Aufl. VIL+293 S. Bonn und Berlin 
1939. Verlag Ferdinand Dümmler. Preis geb. 
32,40 M. 

Gut zehn Jahre nach Erscheinen der ersten 
kommt jetzt mit einigen Berichtigungen, im 
großen Ganzen aber unverändert die zweite Auf- 
lage dieser für das Rechnen mit Maschinen außer- 
ordentlich wertvollen Tafeln heraus. Sie ent- 
halten in der Hauptsache mit einem Argument- 
schritt von 10” die sechs trigonometrischen Funk- 
tionen mit sechs geltenden Ziffern, dazu die ersten 
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Differenzen und Proportionaltafeln, ferner um in 
möglichst allen Gebieten linear interpolieren zu 
können w-ctg w und w-cosee w von 0° bis 
1° 20° mit einem Argumentschritt von 1”. An- 
ordnung, Druck und Ausstattung der Tafeln sind 
ausgezeichnet. 

Dresden. Willers. 114 

Dr. phil. F. DRENCKHAHN, Prof. a. d. Hoch- 
schule für Lehrerbildung, Seestadt Rostock, unter 
Mitarbeit von Dr.-Ing. habil. U. Graf, Prof. a. d. 
Techn. Hochschule Danzig, Geländemathe- 
matik. VII, 211 S. mit 283 Abb. Langensalza- 
Berlin-Leipzig 1939, Verlag von Julius Beltz. Preis 
geb. 7 M. 

Nach dem Vorwort stellt die „Geländemathema- 
tik“ eine Zusammenfassung der Geländekunde und 
der Vermessungskunde vor. Das Buch soll eine 
ohne besonderes Fachstudium verständliche Dar- 
stellung sowohl der theoretischen als auch der 
praktischen Grundlagen der Vermessungs- und Ge- 
ländekunde sein; es wendet sich an „Studenten, die 
Angewandte Mathematik nicht gerade als Spezial- 
fach studieren, an Lehrkräfte an den Schulen, an 
die Besucher von Fachschulen, und schließlich an 
alle diejenigen, die in der Ausübung von Wehr- 
und Geländesport mit der Geländemathematik in 
Berührung kommen“, 

Das Buch ist in elf Teile eingeteilt, in denen die 
Vermessungsebene, Hilfsmittel für Absteckungen 
und Messungen in der Horizontalebene, Hilfsmittel 
für Messungen in der Vertikalebene, elementare Be- 
trachtungen über die bei Vermessungen auftreten- 
den Beobachtungsfehler, Horizontalaufnahmen, Ver- 
tikalaufnahmen, vereinigte Längen- und Höhen- 
messungen — besser wohl vereinigte Horizontal- 
und Vertikalaufnahme — Geländemathematik im 
Geländedienst (Wehrsport), der Mensch als das Maß 
der Dinge — gemeint ist das „Messen“ ohne be- 
sondere Hilfsmittel —. die geometrischen Grund- 
lagen des Luftbildwesens — oder der Luftbild- 
messung — und Kartenentwürfe, die beiden letzten 
Teile von U. Graf, behandelt werden. 

Die ersten sieben Teile beschäftigen sich mit den 
Hilfsmitteln und Verfahren der Vermessungskunde 
in einer für Nichtfachleute bestimmten Darstellung. 
Wie weit bzw. wie wenige weit man dabei im ein- 
zelnen Fall gehen soll, ist eine Frage der persön- 
lichen Auffassung. Ich kann dem Verfasser nicht 
überall beistimmen. Z. B. dürfte m. E., nachdem 
einerseits die einfachsten Zielvorrichtungen in brei- 
tester Form behandelt wurden und andererseits 
Theodolite wie der kleine Reisetheodolit der Firma 
M. Hildebrand und sogar der Theodolit II der 
Firma ©. Zeiß vorgeführt werden, eine Behandlung 
des Zielfernrohrs in seiner Bauart nach Rams- 
den und in der Bauart mit innerer Einstellinse 
unter Beifügung von schematischen Abbildungen 
nicht fehlen. Oder, nachdem begrüßenswerterweise 
auf das Fehlerfortpflanzungsgesetz hingewiesen 
wurde, wäre es angebracht, das Gesetz auch in 
seiner allgemeinsten Form anzugeben. Beim Ni- 
vellierinstrument wird auf dessen Prüfung und Be- 


richtigung eingegangen; aber in einem Umfang, der 
einerseits zu viel und andererseits doch wieder zu 
wenig bietet. Bei der Messung von Vertikalwinkeln 
geht es der Bestimmung des Buches entsprechend 
zu weit, wenn in dem beigefügten Zahlenbeispiel 
das Verfahren der Messung in zwei Fernrohrlagen 
mit Ablesung an zwei Nonien angegeben wird. 
Nachdem beim Ökularfadenentfernungsmesser — in 
etwas umständlicher Weise — die Formel für hori- 
zontale Zielung hergeleitet wurde, dürften die bei- 
den Formeln für nicht horizontale Zielung nicht 
fehlen. Im Zusammenhang mit der optischen Ent- 
fernungsmessung wird — mit einer nicht ganz ge- 
lungenen Erklärung des Gerätes — der stereosko- 
pische Entfernungsmesser mit räumlicher Skala von 
C. Pulfrich aufgeführt; im Anschluß daran 
müßte der militärisch wichtigere stereoskopische 
Entfernungsmesser mit räumlich einstellbarer 
Marke erwähnt werden. 

Von einigen Ungenauigkeiten, die mir aufgefallen 
sind, mögen die folgenden angeführt sein. Beim 
Theodolit wird erwähnt, daß man bei vorhandener 
„Nivellierwendelibelle“ — wie eine solche be- 
schaffen ist, wird nicht gesagt — mit dem Theo- 
dolit auch nivellieren könne; dazu braucht man be- 
kanntlich keine Wendelibelle.. Was unter „Okular- 
schraubendistanzmesser“ verstanden wird, ist mir 
unklar. 

Bei dem mit der Photogrammetrie sich befassenden 
Teil tritt die Frage auf, ob in einem Buch wie in 
dem vorliegenden nicht doch auf die linienweise 
Auswertung von Meßbildern mit Hilfe eines der be- 
kannten Geräte eingegangen werden sollte. 

In den Kreisen, für die das Buch gedacht ist, 
wird es sicher viele Leser finden. 


Stuttgart. P. Werkmeister. 138 


Ferner sind bei der Schriftleitung folgende 
Bücher eingegangen (ausführliche Besprechung 
bleibt vorbehalten): 


RICHARD DOERFLING, Mathematik für 
Ingenieure und Techniker. 533 S. m. 
290 Bildern. München und Berlin 1939, Verlag R. 
Oldenbourg. Preis geb. 9,60 M. 


Dr. phil. HANS ERTEL, Elemente der 
OÖperatorenrechnung mit geophysi- 
kalischen Anwendungen. VI+133 S. m. 
8 Abb. Berlin 1940, Verlag Julius Springer. Preis 
brosch. 7.50 M. 


Dr. B. von JUHOS, Erkenntnisformenin 
Natur- und Geisteswissenschaften 
(Pan-Bücherei, Gruppe: Philosophie, Nr. 20). 57 S. 
Leipzig 1940, Rudolf Birnbach, Abt. Pan-Verlag. 
Preis brosch. 2,50 M. 


Dr. phil. habil. Dr. med. SIEGFRIED KOLLER, 
Vorstand der Statistischen Abteilung des Kerck- 
hoff-Forschungs-Institutes Bad Nauheim, Graphi- 
sche Tafeln zur Beurteilung statisti- 
scher Zahlen. VI +73 S. m. 6 Abb. u. 
15 Tafeln. Dresden und Leipzig 1940, Verlag Theo- 
dor Steinkopff. Preis geb. 10 M. 


NACHRICHTEN 


Gesellschaft für angewandte Mathematik und 
Mechanik. 


Die zum Sonnabend, den 4. Mai 1940, nach Braun- 
schweig einberufene wissenschaftliche Tagung der 
Gesellschaft für angewandte Mathematik und 
Mechanik war außerordentlich gut besucht. Über 
hundert Teilnehmer hatten sich zu einzelnen Vor- 
trägen eingefunden. Das Programm war fast zu 
reichhaltig für eine eintägige Tagung. Die Vor- 
träge und Diskussionen des Vormittages waren 


Fragen der Hydro- und Aerodynamik, die des 
Nachmittages in der Hauptsache Fragen der Elasti- 
zitätstheorie gewidmet. Folgende Vorträge wurden 
gehalten: 

A.Busemann, Braunschweig: Neuere Ansätze 
in der Gasdynamik, 

K. Oswatitsch, Göttingen: Kondensations- 
erscheinungen in Überschalldüsen, 

J. Weissinger, Adlershof: Der Tragflügel 
bei unsymmetrischer Anströmung, 
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H. Sehliehting, Braunschweig: Über den Th. Pöschl, Karlsruhe: Mikroaufnahmen des 


Umschlag laminar-turbulent bei beschleunigter und 
verzögerter Strömung, 

Pretsch, Göttingen: Stabilität einer laminaren 
Grenzschicht bei Druckgefälle. 

H. Görtler, Göttingen: Über eine Instabilität 
(der laminaren Grenzschieht an konkaven Wänden, 


H. Reichhardt, Göttingen: Über Wärme- 
übertragung in turbulenter Strömung, 
Voeltz, Braunschweig: Schräggestellter 
Ringspalt. 


C. Weber, Dresden: Über die Minimalsätze 
der Elastizitätstheorie, 
H. Fromm, Danzig: Führungskraft des rollen- 
den Rades und Schiebebewegung des Fahrzeuges, 
H. Neuber, Braunschweig: Schubmittelpunkt u. 
Querschnittsverwölbung  dünnwandiger Träger 


unterhalb der Beulgrenze, 


Zerreißvorganges von Stahl, 

H. Baumann, Hamburg: Über die Allgemein- 
gültigkeit des Schuler-Prinzipes, 

E. Weinel, Jena: Torsionsbeulung eines 
Plattenstreifens. 172 


Persönliches. 

Es verstarben: 

der frühere o. Prof. für Maschinenbau an der 
Techn. Hochschule Berlin, Geh. Regierungsrat 
Dr. E. Josse, 

der o. Prof. für Wehrtechnik, Physik und Balli- 
stik an der Techn. Hochschule und Honorarprof. 
an der Universität Berlin, General der Artillerie 
Dr.-Ing. K. Becker, Chef des Heereswaffen- 
amtes. 


ZUSCHRIFTEN AN DEN HERAUSGEBER 


v.Borb&ly: Über einen Grenzfall der instationären 
räumlichen Tragflügelströmung. Z. angew. Math. 
Mech. Bd. 18 (1938), S. 319 bis 342. 

In einer in dieser Zeitschrift (Bd. 19 [1939], S. 384) 
erschienenen Äußerung hat Herr Küßner die Er- 
weiterung der Wirbelfadentheorie kritisiert, durch 
die der Dorkune und v. Borbely unabhängig von- 
einander eine Näherungslösung für die instationäre 
Bewegung des endlichen Flügels entwickelt haben. 
Es scheint mir, daß die Arbeit von Sears t), auf die 
Herr Küßner seine Bemerkungen stützt und deren 
Grundlage er als einwandfrei betrachtet, nur eine 
Näherungslösung aufstellt, deren Überlegenheit 
nicht so deutlich ist. Im Gegensatz zu Küßners Be- 
hauptungen, glaube ich, daß auch für den Fall des 
unendlichen Tragflügels mit sinusförmiger Abwind- 
verteilung, diese Theorie durchaus nicht eine exakte 
Berechnung der Zirkulation ermöglicht, da sie, 
zwischen induzierter Geschwindigkeit und Zir- 
kulation, dieselbe Beziehung benützt, die für den 
unendlichen Flügel gilt. Die numerischen Ergeb- 
nisse der groben Berechnung von Sears ermöglichen 
vielleicht keinen sinnvollen Vergleich mit denen der 
strengen Rechnung von v. Borbely. Es scheint 
jedoch, daß der Unterschied zwischen den von mir 
und Sears berechneten Werten des Auftriebes nicht 
sehr groß ist, auch nicht für hohe red. Frequenzen. 

Turin. P. Cicala. 


Erwiderung. Ich habe bereits in meiner ersten 
Zuschrift erwähnt, daß auch Herr Sears eine 
Näherung für schlanke Flügel benutzt. Der wenn 
auch kleine Induktionseffekt infolge endlicher 
Spannweite ergibt sich jedoch bei hoher reduzierter 
Frequenz nach Herrn Cicala und Herrn Sears in 
verschiedener Größenordnung, so daß nicht beide 
Theorien zugleich richtig sein können. Die Ursache 
dieser Diskrepanz liegt darin, daß Herr Cicala seine 
Größe K willkürlich gleich der Summe der gebun- 
denen und freien Wirbel längs der Flügeltiefe setzt, 
während K in einer konsequenten Theorie, die ich 
demnächst veröffentlichen werde, ein komplizierter 
Ausdruck ist. 


Göttingen. H. G. Küßner. 165 


1) W. R. Sears: Proe. of the 5. Intern. Congress of 
Applied Mechanies 1938. 


Berichtigung 


zu „Interpolationsverfahren zur Berechnung von 
Flugbahnscharen . . .“ von H. Athen. Z. angew. 
Math. Mech. Bd. 19 (1939), S. 361. 


Herr Dr. Stange, Essen, auf dessen demnächst 


in dieser Zeitschrift erscheinende Ergänzungen zu 
meinem Aufsatz ich wegen ihrer Zweckmäßigkeit 
bei praktischen Rechnungen besonders hinweise, 
machte mich darauf aufmerksam, daß auf S. 363 im 

_ht% . 
Ay für ei 


Anschluß an Formel (5) in n 


Vorzeichenfehler enthalten ist: es muß heißen 
2 


Berlin-Frohnau. H. Athen. 170 


Nachtrag zu der Arbeit: Eigenschwingungen 
gedrückter Kreisplatten. 


Z. angew. Math. Mech. Bd. 20 (1940), S. 37 bis 44. 


Herr Professor Dr. Federhofer, Graz, macht mich 
darauf aufmerksam, daß Eigenschwingungen ge- 
drückter Platten schon in seiner Arbeit „Biegungs- 
schwingungen der in ihrer Mittelebene belasteten 
Kreisplatte“, Ingenieurarchiv Bd. 6 (1935), S. 68 bis 
74, behandelt sind. In dieser Arbeit, die mir leider 
entgangen ist, werden sowohl auf Druck wie auf 
Zug beanspruchte, am Rande eingespannte Kreis- 
platten jedoch ohne federnde Stütze im 
Mittelpunkt untersucht, d. h. in meiner Be- 
zeichnung der Fall o—=0. Mir kam es vor allem 
darauf an, zu zeigen, daß man auch in komplizierten 
Fällen Gleichung und Randbedingung aus dem 
Hamiltonschen Prinzip einfach ableiten kann und zu 
untersuchen, wie sich die Frequenz der Grund- 
schwingung durch die federnde Stützung ändert. 
Herr Professor Federhofer behandelt in der oben 
erwähnten Arbeit für die nicht gefederte Platte 
außer der Grundschwingung die Oberschwingungen 
mit 2 Knotenkreisen (2 — 0;1;2) und keinem oder 
einem Durchmesser als Knotenlinie, ferner solche 
mit zwei Durchmessern als Knotenlinie und keinem 
Knotenkreis. Er weist in seinem Schreiben darauf 
hin, daß man im letzten Fall bei sehr starker Fede- 
rung und starkem Druck noch einen kleinen Be- 
reich erhält, in dem die Frequenz kleiner ist als die 
der Schwingungen ohne Knotenlinie, so daß diese 
also hier erst an dritter Stelle kommen. Schließlich 
bemerkt Herr Federhofer noch, daß die am Schluß 
meiner Arbeit für den dort nur interessierenden 
Fall eines Durchmessers als Knotenlinie ange- 
wandte Knickbetrachtung im allgemeinen Fall zu 
der von Herrn -Nadai 1915 aufgestellten Knick- 
gleichung 

Jp+1(2)=0 
führt. 
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